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1.  INTRODUCTION 


Postflight  processing  of  spread-spectrum  signals  often  requires 
high-rate  data  recording.  One  way  to  reduce  the  recording  rate  require- 
ment is  to  use  single-bit  quantization  (i.e.,  hard  limiting)  of  the  sam- 
pled signal  in  the  recording. 

The  performance  of  an  aided  digital  postflight  receiver  for  car- 
rier and  code  tracking  when  single-bit  quantization  is  utilized  in  the 
data  recording  is  here  analyzed.  The  assumed  signal  structure  is  that 
used  in  the  SATRACK-GPS  program. 

In  the  SATRACK  scheme,  signals  originating  from  N in-view  satel- 
lites (N  £ 7)  ar*1  transponded  by  a missile  to  receivers  on  range  ships. 
The  N satellites  consist  of  N - 1 Global  Positioning  System  (GPS)  satel- 
lites and  one  pseudosatellite  (i.e.,  a range  ship  imitating  a satellite). 
The  signal  received  by  the  range  ship  is  assumed  to  be  single-bit  quan- 
tized before  being  sampled  and  recorded  for  postflight  processing.  Tel- 
emetry and  ground  station  data  are  used  to  aid  the  various  tracking  loops 
in  the  postflight  receiver,  permitting  the  bandwidths  of  these  loops  to 
be  narrowed. 

The  body  of  this  report  combines  the  results  of  Appendixes  D,  E, 
and  F,  where  investigations  are  made  into  the  performance  of  component 
loops  of  the  receiver.  Detailed  diagrams  of  the  component  loops  of  the 
receiver  are  given  in  Appendix  A.  An  example  of  the  performance  of  the 
receiver  for  a given  set  of  input  parameters  and  a given  shipboard  sam- 
pling rate  is  shown  in  Appendix  B.  The  performance  of  an  analog  model 
of  this  receiver  (continuous-time  input,  no  quantization)  is  analyzed  in 
Appendix  C.  Appendixes  C,  D,  E,  and  F report  the  results  of  separate  in- 
vestigations carried  out  at  different  stages  of  the  overall  work.  Thus 
at  times  the  notations,  system  parameters,  and  receiver  diagrams  given 
in  these  appendixes  differ  from  those  given  in  the  text  and  Appendix  A. 

The  analysis  in  Sections  2 through  7 is  for  an  uninterrupted  input 
signal.  In  actuality,  the  SATRACK  missile  transponder  switches  between 
two  pairs  of  receiving  antennas  every  1.152  ms.  The  received  signal  from 
each  antenna  pair  is  tracked  by  its  own  receiver,  hence  the  Input  signal 
to  the  receiver  is  on  for  1.152  ms,  off  for  1.152  ms,  etc.  The  effect 
of  this  antenna  switching  and  shipboard  receiving  antenna  polarization 
are  both  discussed  in  Section  8. 
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2.  DESCRIPTION  OF  THE  SHIPBOARD  RECEIVED  SIGNAL 
IN  THE  SATRACK  PROGRAM 


In  the  SATRACK  scheme,  signals  originating  from  N * 1 in-view 
Global  Positioning  System  satellites  (N  £7),  and  one  pseudosatellite, 
are  transponded  by  a missile  to  receivers  on  range  ships.  The  signal 
received  by  a given  range  ship  has  the  form, 

N ( , 

x(t)  - 53  j ■^Ai(t)PRNi(t  - X1)m1(t  - X^sinto^t  •+  $t(t)]  j+  n(t), 

(1) 

where  Xj  denotes  the  total  time  delay  from  signal  transmission  by  satel- 
lite (i)  to  final  signal  reception  by  the  ship;  the  function  PRN^(t)  is 
a pseudorandom  sequence  (Gold  code)  of  ±1*8  with  a bit  (chip)  length  of 
A * 10"6/1.023  second,  and  a sequence  repetition  rate  (epoch  rate)  of  1 
kHz  (i.e.,  1023  chips  per  epoch);  the  function  m^(t)  is  also  a sequence 
of  ±l*s,  which  among  other  things  carries  satellite  orbit  data;  the  chip 
length  of  m^(t)  is  Am  • 20  ms;  the  phase  function  $i(t)  is  the  result  of 
doppler  shifts  caused  by  satellite,  missile,  and  range  ship  movement;  and 
the  noise  function  n(t)  is  white  noise,  which  is  bandlimited  by  the  ship- 
board receiver,  to  ±1.073  MHz  about  the  missile  transmission  frequency, 
fQ  * 2200  MHz.  A ±1.073  MHz  bandwidth  is  chosen  to  allow  room  for  the 
code  width  of  1.023  MHz,  and  a maximum  doppler  shift  of  ±50  kHz.  The 
power  spectral  density  (PSD)  of  n(t)  is  given  by  N0/2  » io~17.2  (mw/Hz) . 
The  noise  function  may  be  written 

n(t)  - ^?[nI  (t)sln  wot  + ^(tjcos  , 

where  ni (t)  and  n2(t)  are  independent  white  Gaussian  noises,  of  power 
spectral  density  N0/2,  which  are  bandlimited  to  jf|  aC  1.073  MHz.  The 
carrier  power  for  the  ith  signal,  A|,  can  range  between  10~14 .7  atuj 
10-12.5  mw.  The  code  sequences,  PRNi(t),  and  their  epoch  times  (as 
transmitted  by  the  satellites)  are  known  by  the  postflight  receiver,  al- 
though the  delays  (X^)  are  not  known  a priori.  However,  telemetry  and 
ground  station  data  are  available  to  the  receiver,  and  are  converted 
into  aiding  information  for  the  various  receiver  tracking  loops.  The 
data  mj(t)  may  be  partially  known  a priori,  but  will  be  treated  here  as 
totally  unknown. 
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3.  SHIPBOARD  PREPROCESSOR 


The  carrier  signal  received  by  the  ship  is  first  heterodyned  down 
to  100  MHz  and  then  input  to  the  preprocessing  scheme  shown  in  Appendix 
A,  Fig,  A-l.  The  input  signal  to  the  preprocessor  is  heterodyned  to  a 
center  frequency  of  100  kHz.  The  doppler  shift  is  expected  to  be  less 
than  50  kHz  on  either  side  of  this  center  frequency.  The  2.146  MHz  sig- 
nal bandwidth  is  necessary  because  of  the  code  width  (±1.023  MHz)  and  the 
doppler  shift.  A 100  MHz  low-pass  filter  is  used  to  reject  the  double- 
angle signal  and  noise  resulting  from  the  mixing  operation.  The  ship- 
board IF  frequency  was  chosen  at  100  MHz  to  avoid  averaging  the  output 
of  the  mixer  in  the  low-pass  filter  over  any  time  period  larger  than 
that  over  which  the  low  frequency  component  of  this  output  (which  has  a 
bandwidth  of  1.073  MHz)  will  charge  by  0.25  radian  (Ref.  1).  The  I (in- 
phase)  and  Q (quadrature)  channel  signals  recorded  on  the  ship  are  given 
by 

I:  Sign  | “ Xi)cos^i(tn^  + “^n  ~ ^ 

4-  n1(tn)cos(w2tn  - 6)  - n2(tn)sin(u)2tn  - 9)j 
Q:  Sign  j - X^PRN^  - X^sin^^)  4-  w2tn  - 0] 

4-  n^(tn)sin(u2tn  - 9)  4-  n2(tn)cos (w2tn  - 0)  | . 


Ref,  1.  W.  E.  Larimore,  "Design  and  Performance  of  a Second- 
Order  Digital  Phase-Locked  Loop,"  Symposium  on  Computer  Processing  in 
Communications , Polytechnic  Institute  of  Brooklyn,  1969,  pp.  343-356. 
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4.  POSTFLIGHT  DIGITAL  RECEIVER 


The  signal  recorded  on  the  range  ship  Is  processed  at  APL  in  the 
postf light  digital  receiver.  The  object  of  the  receiver  is  to  track  the 
carrier  doppler  shifts  of  the  signals  from  the  various  satellites,  and 
the  different  delays  of  the  PRN  sequences.  Telemetry  and  surface  station 
data  from  the  flight  are  provided  to  the  postflight  receiver  to  aid  in 
this  task. 

The  accuracy  of  the  aiding  signals  determines  the  dynamics  that 
the  receiver  loops  must  track.  Given  the  order  of  a loop,  the  dynamics 
to  be  tracked  essentially  present  a lower  bound  on  the  usable  loop  band- 
width. It  is  assumed  that  the  accuracy  of  the  aiding  signals  is  suffi- 
cient to  ensure  that  the  loop  bandwidths  chosen  in  the  examples  are  above 
these  lower  bounds. 

Each  of  the  N delays  {X^}  and  doppler  shifts  {^(t)}  are  tracked 
with  a separate  receiver.  The  ith  receiver  uses  its  knowledge  of  the 
ith  code  sequence  PRNj[(t)  to  separate  its  signal  component  from  the  other 
N - 1 signal  components  in  x(t)  (see  Eq.  (1)).  The  interference  in  the 
ith  receiver  from  the  other  N - 1 signal  components  in  x(t)  will  affect 
the  performance  of  the  receiver,  and  will  superpose  its  effect  with  that 
caused  by  the  noise  n(t) . The  effects  of  interference  are  deterministic 
in  nature,  and  can  perhaps  be  attenuated  by,  for  example,  processing  the 
strongest  signal  first,  and  then  compensating  the  receivers  for  the 
other  signals  with  the  information  obtained  about  the  strongest  signal, 
etc.  Concentration  here  will  be  on  the  effects  of  input  noise.  Thus 
we  will  deal  with  the  signal  x(t)  as  though  it  were  in  fact  given  by 
xi(t),  a fictitious  signal,  defined  as 

Xi<t)  - ■^2A1(t)PRN1(t  - X^m^t  - X^sinlu^t  + ^(t)]  + n(t)  (2) 

The  ith  postflight  receiver  has  the  form  shown  in  Appendix  A, 

Fig.  A-2.  Note  that  the  computer  aids  obtained  from  telemetry  and  ground 
station  data  are  used  in  both  stages  of  the  receiver.  The  first  stage 
of  the  receiver  estimates  the  function  mi(t  - X^)  and  stores  it  on  tape. 
This  modulation  is  then  stripped  off  the  input  signal  for  use  by  the  sec- 
ond stage.  The  first  stage  also  obtains  information  on  code  epoch  syn- 
chronization and  carrier  phase,  which  may  be  used  in  the  second  stage, 
along  with  (or  combined  with)  the  computer  aids  obtained  from  telemetry 
and  ground  station  data.  The  second  stage  of  the  receiver  uses  the  out- 
puts of  the  first  stage  to  aid  the  precision  estimation  of  $i(t)  and  the 
delay  X^.  If  the  information  on  code  epochs  and  carrier  phase  is  com- 
bined with  the  computer  aids,  off-line  processing  is  required.  For  this 
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reason  we  show  - A^)  being  put  on  tape,  and  the  two  stages  of  the 

receiver  do  not  run  concurrently.  The  various  components  pictured  in 
Fig.  A-2  are  given  in  more  detail  in  other  figures  of  Appendix  A.  As 
to  the  performance  of  the  receiver,  the  relationships  given  in  the  fol- 
lowing discussion  are  valid  for  signal  parameters  within  a few  orders 
of  magnitude  of  the  SATRACK  values  (yielding  quasilinear  performance  in 
the  tracking  loops).  Very  low  signal  levels  cause  threshold  effects 
that  are  not  accounted  for  here.  High  signal  levels  cause  a decrease 
in  the  input-noise  effects,  bringing  them  down  to  the  levels  of  the  ef- 
fects of  noises  generated  internally  by  the  tracking  loops  (e.g.,  see 
Appendix  D) . Signal  power  levels  comparable  to  the  input  noise  power 
(in  a ±1.073  MHz  bandwidth)  cause  a loss  of  the  dither  provided  by  the 
input  noise,  and  a consequent  increase  in  nonlinear  phenomena.  A de- 
tailed example  of  the  type  of  analysis  used  to  obtain  the  following  re- 
sults is  given  in  Appendix  E. 


- 11  - 


THl  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LMJHC.I.  MARYIAND 


5.  PERFORMANCE  OF  THE  FIRST  STAGE  OF  THE  DIGITAL  RECEIVER 


The  main  purpose  of  the  first  stage  of  the  digital  receiver  is  to 
estimate  the  data  modulation,  mf(t  - Af).  The  first  component  of  this 
stage  of  the  digital  receiver  is  the  noncoherent  delay-locked  loop  (DLL) 
(Appendix  A,  Fig.  A-3) . The  purpose  of  this  loop  is  to  lock  onto  the 
FRN  modulation  and  to  remove  it  from  the  input  signal  (see  Eq.  (8)). 

The  input  signal  is  first  heterodyned  down  to  0 frequency  by  an  aided  sig- 
nal designed  to  also  remove  as  much  doppler  shift  as  possible.  The  per- 
formance of  this  loop  is  analyzed  in  Appendix  F.  The  rms  delay  error, 
ei(t)  ■ A ^ (t ) - A^(t),  at  the  loop  output  is  given  by 


where  is  the  one-sided  bandwidth  of  lowpass  filters  preceding  the 
squarers.  This  width  must  be  sufficient  to  pass  both  the  data  (50  Hz) 
and  the  aided  doppler  shift  (say  20  Hz,  although  preliminary  simulations 
show  that  3 Hz  is  probably  sufficient).  Thus  Bw  could  be  chosen  to  be 
70  Hz  (a  1/140  second  averager).  The  factors  of  0.9029  represent  the 
degradation  in  power  caused  by  passing  the  data  through  a filter  that 
passes  only  the  first  major  lobe  of  the  data  power  spectrum.  This  num- 
ber represents  the  value  of  E[m^Lp]2  where  miLp(Lt  - Aj)  is  the  filtered 

data.  The  value  0.9029  is  obtained  by  modeling  the  data  as  having  a tri- 
angular autocorrelation  function  (see  Appendix  F).  If  the  data  were  in 
fact  a 40  ms  square  wave  (the  worst  possible  case),  this  factor  would 
change  to  0.81059.  is  the  two-sided  loop  noise  bandwidth.  This  pa- 

rameter may  be  chosen  to  give  the  desired  error  performance,  with  the 
constraint  that  too  small  a value  of  BnI  will  allow  noises  generated 
within  the  loop  (not  discussed  here)  to  become  significant  factors.  Note 
that,  as  discussed  in  Appendix  D,  local  clock  phase  jitter  is  not  a sig- 
nificant noise  source,  because  of  the  synchronous  control  of  the  local 
clock,  the  loop  sampling  mechanism,  and  the  reading  of  the  input  tape. 

The  factor  of  lO^* 19991  a function  of  the  sampling  rate,  which  is  as- 

sumed here  to  be  2.5  MHz.  In  Fig.  1 the  variation  of  this  factor  with 
sampling  rate  is  plotted.  Note  that  at  2.5  MHz  the  factor  is  1.9991  dB 
below  the  line  in  Fig.  1 for  analog  continuous  time  processing,  account- 
ing for  our  factor  of  10^* 19991  fn  Eq.  (3).  An  ideai  shipboard  bandpass 
filter  (±1.073  MHz)  was  assumed  in  the  analysis.  A comparison  of  Eq.  (3) 
with  the  result  of  Gill  (Ref.  2)  for  analog  continuous  time  processing 

Ref.  2.  W.  J.  Gill,  "A  Comparison  of  Binary  Delay-Lock  Tracking 
Loop  Implementations,"  IEEE  Transactions  on  Aerospace  and  Electronic 
Systems,  Vol.  AES-2,  No.  4,  July  1966,  pp.  415-424.  " 
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Fig.  1 Performance  Degradation  Factors  as  a Function  of  Sampling 
Rate:  (a)  Factor  for  Eq.  3 (i.e.,  noncoherent  delay-locked 
loop);  (b)  Factor  for  Eq.  11  and  Eq.  18  (i.e.,  suppretsed- 
carrier  tracking  loop  or  phaw-locked  loop) 
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shows  that  the  degradation  caused  by  the  shipboard  bandpass  filtering 
and  single-bit  processing  is  embodied  in  the  factors  of  i()0*  1^991  in 
Eq.  (3).  Th«i  multiplication  of  the  band  limited  input  noise  by  the 
reference  code  causes  a 0.4436  dB  decrease  in  noise  power  spectral  den- 
sity. The  effect  of  this  noise  decrease  is  cancelled  by  a 0.4436  dB  de- 
crease in  the  slope  of  the  code  cross-correlation  function,  which,  in 
turn,  is  caused  by  band  limiting  of  the  input  signal  (to  be  discussed). 
Thus  shipboard  band  limiting  produces  no  net  effect.  Consequently  the 
factors  of  io^* 19991  gq,  represent  the  effects  of  single-bit  pro- 
cessing alone. 

We  find  it  useful  later  to  calculate  E[RpR}j^(e^)  ] , where  Rprn^(0 
is  the  ith  code  autocorrelation  function.  Modeling  as  Gaussian, 


where  a(e^)  is  obtained  from  Eq.  (3). 

Because  of  the  single-bit  processing,  the  mean  output  of  the  de- 
lay detector  is  decreased  below  the  output  that  would  be  observed  if 
there  were  no  input  noise.  The  factor  of  reduction  is  inversely  propor- 
tional to  the  aignal/nolae  variance  ratio  (SNRi)  on  the  input  signal, 
for  the  low  SNR  values  that  concern  us.  In  our  case,  with  a -147  dBm 
signal,  SNRj  is  given  by 


(^)<2)(^dJh)  (r)  (2)<2-U6  * 10<i>  (10'17-2H2)(2-1«  « TO6) 


- -41.3266  dB,  (5) 

and  we  find  that  the  scale  factor  with  this  ratio  is  26490.8.  Thus,  if 
the  signal  level  were  in  fact  -135  dBm,  then  SNR^  would  be  -40.1266,  and 
the  scale  factor  would  be 

(26490.8)  io+A*01266  “4*13266  _ 20095t3  (6) 

We  will  discuss  the  internal  loop  signals  as  though  sampling  in 
the  loops  were  being  carried  out  at  the  same  rate  as  the  shipboard  sam- 
pling. In  fact,  the  loop  sampling  rates  are  variable,  and  the  ship- 
board recorded  signal  is  held  at  a constant  level  between  its  samples. 
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With  2.5  MHz  sampling  and  a counter  cycle  of  1/139.997  second, 
the  maximum  count  that  could  be  observed  at  the  output  of  Lowpass^in 
Appendix  A,  Fig.  A-3  is  (2.5  x 1Q6)  (1/139. 997) (2)  - 35  715  counts.  In 
the  no-noise  case,  the  signal  that  would  be  observed  at  this  point  is 


(35715)  n. 


*1? 


(t  - x1)  - xt)  |i-  j -it  + [^(t)  + *|nod  2lr  j(f)  j 


(7) 


where  (t)  ^ ^(t)  - “i^p^  “ *i^  represents  the  filtered 

data,  and  is  the  ith  code  autocorrelation  function.  The  third 

factor  in  Eq.  (7)  is  a cosine  wave  with  triangular  lobes.  With  input 
noise  present,  and  a -147  dBm  signal,  there  is  a scale  factor  of 

“^26490.8  - 162.760  associated  with  this  signal  [at  (^(t)  ■ 0],  So 

we  would  observe 


35715  (^26490.8) 


-1 


m 


(t  - \t)  Rpnj,  U - cos 

Lr  1 


♦ . (t) 
di 


+ na±(t) 


(8a) 


The  transformation  of  the  triangular  wave  to  a cosine  is  discussed  in 
Appendix  D.  The  output  of  Lowpass^^is  given  by 


35715  (V 26490.8)  m (t  - XJ  RpRN  (^  - X±)  sin 

LP  i 


V£> 


nbi(t>’ 


(8b) 


where  noise  terms  nai(t)  and  nfci(t)  are  in  quadrature.  The  signals  in 
Eq.  (8)  are  multilevel  quantized  and  take  new  values  every  (139.997)-! 
second.  They  will  be  referred  to  collectively  as  yn(t).  The  signal/ 
noise  density  ratio  in  Eqa.  (8a)  and  (8b)  is  given  by 


(10 


-0.24427 


) (1/2)  A 


2 

i 


N 

o 

2 


(The  signal/noise  density  ratio  of  these  signals  without  shipboard  hard 
limiting  and  sampling  would  have  been  (1/2)  (A|)  (N0/2)“1).  The  varia- 
tion of  the  factor  10"0* 24427  with  shipboard  sampling  rate  is  shown  in 
Fig.  1.  Thus  we  may  determine  t£e  variance  of  na^(t)  that  would  be  ob- 
served at  the  output  of  Lowpass(l). 

signal  jjiwgr  . (1/2)  [ (35715) ( V 26490. g)"1! 2 (10  °’24427)  (1/2)  (A^) 

noise  deirtity  (139.997)'1  o2In  ,(t»  ' % 

2 
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thus 

o2(nal(t)]  - c2tnbi(t)]  - 37410.6  counts2  (9) 

The  next  phase  of  the  tracking  is  a suppresaed-carrier  tracking  loop. 
The  Intent  is  to  track  the  phase  angle  in  Eq.  (8)  in  order  to  remove  the 
sinusoids  from  the  data  modulation.  (Recall  that  the  PRN  codq  has  been  re- 
moved by  the  noncoherent  DLL.) 


The  suppressed-carrier  tracking  loop  is  shown  in  Appendix  A,  Fig.  A-4. 
The  output  signals  of  the  noncoherent  delay-locked  loop  are  first  heterodyned 
up  to  a 25  Hz  IF  frequency  so  that  the  suppressed-carrier  loop  will  not  have 
to  track  near  zero  frequency.  The  loop  itself  is  a combination  of  a squaring 
loop  and  a Costas  loop.  The  combination  is  needed  because  the  input  signal 
from  the  ship  is  recorded  at  low  enough  frequency  to  have  noise  terms  over- 
lapping zero  frequency.  The  output  of  this  loop  is  two  signals  given  by 


y2i(t) 


{costuijt  - 4»d  (t)  + n3^(t))  I 

sinlw^t  - (t)  4 n3i^t)l j » 


(10) 


where  n3i(t)  is  the  phase  tracking  error.  The  rms  value  of  n31(t)  is  given 
by  (see  Appendix  F) 


H <10°-2“”) 

0 


2A‘<E[RFRM^Ce1)l r (0.9029) 


/ »„<ioo-»m)  v 

’uAjtEtR^^)]}2  (0.9029)1 


(2Bw) 


1/2 


N2  I 


rad 


(ID 


where  Bfl2  is  the  two-sided  loop  bandwidth  of  the  suppressed  carrier  loop,  and 
E(RpRN1(e1)]  was  evaluated  in  Eq.  (4).  The  effect  of  changes  in  sampling 

rate  on  Eq.  (11)  is  given  in  Fig.  1,  analogously  to  the  effect  on  Eq.  (3). 

The  factor  io° -'4427  Eq,  (ll)  represents  0.4436  dB  more  degradation  than 
does  the  factor  100.19991  in  Eq,  (3).  This  difference  is  the  result  of  ship- 
board bandpass  filtering  of  the  input  signal.  Here  the  effects  of  this  oper- 
ation are  analyzed  differently  and  more  accurately  than  they  are  in  Appendix 
F.  The  bandpass  filtering  affects  the  multiplications  of  the  input  signal  by 
reference  codes  in  the  receiver.  The  cross  correlation  of  the  band  limited 
input  code  with  the  reference  code  is  a distorted  version  of  the  code  auto- 
correlation function.  The  degradation  in  the  peak  value  of  the  autocorrela- 
tion function  affects  the  operation  of  any  carrier  tracking  loops  that  follow 
code  removal.  The  degradation  in  the  slope  of  the  autocorrelation  function 
(iu  the  neighborhood  of  40.54  and  -0.54  from  the  peak  value)  affects  the  op- 
eration of  the  delay-locked  loops.  For  an  ideal  bandpass  filter  at  the  chip 
rate,  the  degradation  in  peak  values  is  found  to  be  0.8872  dB.  Thu  degrada- 
tion in  slop*  is  calculated  for  a causal  first-order  bandpass  filter.  (The 
causality  seems  to  be  important  for  this  result.)  The  slope  degradation  is 
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approximately  one-half  of  the  peak  value  degradation,  or  0.4436  dB.  The  re- 
sult given  in  Appendix  F (i.e.,  1.999  dB  degradation)  included  a 0.4436  dB 
degradation  for  both  the  phase-locked  loops  and  the  delay-locked  loops.  This 
is  accurate  for  the  delay-locked  loops,  but  another  0.4436  dB  must  be  sub- 
tracg|^  |or  the  phase-locked  loops.  This  explains  the  extra  factor  of 


10°  • 


in  Eq.  (11),  compared  with  Eq.  (3). 


Equation  (11)  may  be  compared  with  Eq.  (10.52)  in  Ref.  3 where  a sup- 
pressed carrier  loop  with  analog  continuous-time  processing  is  analyzed.  We 
see  then  that  the  factors  of  10®* 24427  represent  the  degradation  caused  by 
shipboard  bandpass  filtering  and  single-bit  quantization.  As  described  above, 
the  effects  of  shipboard  band  limiting  are  a 0.4436  dB  decrease  in  noise 
power  spectral  density  and  a 0.8872  dB  decrease  in  code  cross  correlation 
peak.  Thus  we  may  attribute  a net  degradation  of  0.4436  dB  to  the  shipboard 
bandpass  filtering.  This  leaves  a factor  of  io^»24427  -0.04436  „ 1.9991  dg 
caused  by  the  single-bit  processing  alone,  as  in  Eq.  (3). 


It  is  useful  to  evaluate  E costn^^t)]  later.  Modeling  1131  (t)  as 
Gaussian,  we  obtain 


E cos[n3i(t)]  - exp 


[n_,  (t) ] 

JA 

2 


(12) 


where  o2[n3i(t)]  is  obtained  from  Eq.  (11).  The  outputs  of  the  carrier 
loop  are  mixed  with  the  quantities  in  Eq.  (8)  as  shown  in  Appendix  A, 
Fig.  A-4.  This  removes  the  carrier  from  the  data.  The  result  is  given 
by 


y3i< C)  - (t  - \±)  cos  n3i(t)  RpRN  [ejL(t)]  + n4i(t)  , (13) 

LP  i 

where  n^(t)  is  a lowpass  Gaussian  process  of  70  Hz  bandwidth  and 
100. 24427  (N0/2A^)  power  spectral  density.  There  is  a 180°  phase  ambigu- 
ity in  yn(t) ; thus  the  data  sequence  and  its  inverse  are  mutually  indis- 
tinguishable. This  has  no  effect  on  receiver  performance,  however,  and 
will  be  ignored.  Note  that  the  data  modulation  on  the  signal  transmitted 
by  the  satellite  is  the  source  of  this  180°  phase  ambiguity.  Because  of 
differential  encoding  of  the  data,  message  estimation  is  also  not  af- 
fected. Appendix  C describes  both  the  technique  for  synchronizing  the 
data  clock  (data  transitions  occur  once  every  20  code  epochs,  and  are 
aligned  with  the  code  epochs)  and  the  technique  for  estimating  the  bits 
in  m^(t  - Xi)  given  Eq.  (13). 

Ref.  3.  A.  J.  Viterbi,  Principles  of  Coherent  Communication.  McGravP 
Hill  Book  Co. , New  York,  1966. 
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We  obtain  the  following  bit-estimation  error  probability: 


P<£)  2 P 


tn^(t  - \j)  i »i(t  - 


1/2 


bit  Ik  bit  Ik 


1 ' crf(V^T{EtRPRN1<Ci>1}  jEtc<>8  n3i(t))j) 


. (14) 


where 


Ki  ’ (n  <100'24427)  ) <V 


(15) 


(Recall  that  (A_  ■*  20  me)  and  that  o(ei)  ia  given  by  Eq.  (3).)  The  quanti- 
ties in  Eq.  (14)  are  found  in  Eqa.  (15),  (4),  and  (12),  respectively.  The 
erf (•)  ia  the  tabulated  error  function.  (In  Eq.  (14)  any  intersymbol  inter- 
ference effects  and  the  fact  that  the  signal  in  Eq.  (13)  only  changes  value 
every  1/140  second  have  been  ignored.)  It  is  useful  to  evaluate 
E[mi(t  - X*)  mi(t  - X^)].  First  we  must  reaLize  that,  because  of  the  code 
epoch  estimation  error,  the  probability  that  raj(t  - X^)  + n>i(t  - X^)  at  any 
specific  time  is  not  equal  to  the  bit-estimation  error  probability.  Instead, 


. » /a  - o(c  ) \ 

P(m  4 m)  « Pro'o  [ra^(t  - « I P (f>  1 ^ jj j + 1/2  I -jj J 


(16) 


We  may  then  write 

^(ram)  - - X^)  ra^t  - X^]  -lx  p(m  - m)  - 1 x P(m  m)  (17) 


Equation  (17)  yields  an  important  quantity  for  determining  the  operation  of 
the  second  stage  of  the  digital  receiver.  Incidentally,  this  quantity  can 
be  seriously  degraded  by  the  presence  of  cycle  slips  in  the  suppressed-carrier 
tracking  loops.  Thus  this  phenomenon  must  be  avoided  as  much  as  possible. 

Some  important  simplifications  possible  in  Fig.  A-4  (Appendix  A)  and  the  data 
estimation  are  given  in  Appendix  G. 

An  example  of  the  determination  of  the  performance  of  the  entire  re- 
ceiver, given  a set  of  sample  signal  parameters  and  a sampling  rate,  is  done 
in  Appendix  B.  In  the  next  section,  the  second  stage  of  the  digital  receiver 
is  discussed. 
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6.  PERFORMANCE  OF  THE  SECOND  STAGE  OF  THE  DIGITAL  RECEIVER 


The  phase-locked  loop  and  coherent  delay-locked  loop  that  com- 
prise the  second  stage  of  the  postflight  receiver  for  x^(t)  are  here  dis- 
cussed. This  stage  performs  the  precision  tracking  of  the  doppler  shift 
(4>i(t))  and  the  delay  (A^).  The  operation  of  this  stage  will  be  aided 
by  the  value  of  m^(t  - A^)  determined  by  stage  1,  and  by  the  combination 
of  carrier  and  code  aids  formed  from  both  the  output  of  stage  1 and  the 
telemetry  and  ground  station  data  (see  Appendix  A,  Fig.  A-2).  The 
errors  in  this  stage  can  be  determined  by  iteratively  solving  a pair  of 
simultaneous  equations.  The  digital  phase-locked  loop  is  shown  in  Appen- 
dix A,  Fig.  A-5.  If  given  the  value  of  the  delay  error,  eai(t)  « A-^t) 

- Aia(t)»  of  the  coherent  delay-locked  loop,  then  the  rms  value  of  the 
phase  error,  <t>ej(t)  » ^(t)  - £ia(t)»  of  the  phase-locked  loop  is  given 
bv 


<>Uel(t)] 


N (10°*24427) 
o 


A.  


rad,  (18) 


where  is  the  loop  noise  bandwidth  of  the  phase-locked  loop.  E[mm] 
was  determined  in  Eq.  (17),  and  E{Rpm^  [eai(t)]> may  be  obtained  by  apply- 
ing Eq.  (4)  to  the  rms  value  of  eai(t).  The  factor  of  i()0‘ 24427 
again  a function  of  sampling  rate,  as  shown  in  Fig.  1.  Given  4>ei(t)»  the 
rms  value  of  eaj[(t)  in  the  delay-locked  loop  (Appendix  A,  Fig.  A-6)  may 
be  obtained  from 


aleai(t)] 


N (10°‘ 18067 
o ' 


(Elcos  *el(t))> 


[E(mo) ] 


-)  W 


where  Bj^  is  the  two-sided  loop  noise  bandwidth  of  the  delay-locked  loop, 
E cos  <}>ei(t)  may  be  determined  from  the  rms  value  of  <frei(t)  by  using 
Eq.  (12).  The  factor  1()0‘18067  vand  for  2.5  MHz  sampling.  The 
variation  of  this  factor  with  sampling  rate  is  given  in  Fig.  2.  Note 
that  for  sampling  at  20  MHz  the  degradation  is  0.5043  dB.  This  checks 
well  with  the  result  of  Spilker  (Ref.  4),  who  finds  a 0.414  dB  degrada- 
tion in  binary  delay-locked  loop  performance  with  continuous  time  proc- 
essing following  hard  limiting.  Also  note  from  Fig.  1 that  the  phase- 
locked  loop  degradation  with  20  MHz  sampling  is  1.2103  dB,  As  described 

Ref.  4.  J.  J.  Spilker,  Jr.,  "Delay  Lock  Tracking  of  Binary  Sig- 
nals," IEEE  Trans,  on  Space  Electronics  and  Telemetry,  Vol.  SET-S,  March 
1963,  pp.  1-8. 
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Fifli  2 Parformanca  Degradation  Factor  in  E<*.19,  at  a Function  o! 
Sampling  Rata 
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in  Section  5,  0.4436  dB  of  this  degradation  is  due  to  the  shipboard  band- 
pass filtering,  leaving  0.7667  dB  caused  by  single-bit  quantization. 

This  checks  well  with  the  results  of  Springett  and  Simon  (Ref.  5),  who 
found  a 0.645  dB  degradation  under  somewhat  similar  circumstances  with 
continuous  time  processing.  Equations  (18)  and  (19)  are  simultaneous 
equations  for  a[$ei(t)]  and  a[eai(t)].  These  two  equations  can  be  solved 
by  starting  with  E[RpR^(ei)]  wl,  obtaining  o[6ei(t)l  from  Eq.  (18), 
substituting  this  into  Eq.  (19)  to  get  o[eai(t)],  and  proceeding  itera- 
tively until  o[$ej.(t)}  and  o[eai(t)]  converge. 

Incidentally,  Eq.  (19)  is  derived  by  adding  0.4436  dB  to  the 
degradation  given  in  Appendix  E,  This  is  necessitated  by  the  bandpass 
filtering  of  the  input  code,  not  accounted  for  in  Appendix  E.  Equation 
(18)  is  obtained  from  Appendixes  D and  F.  The  degradation  determined 
in  Appendix  D was  calculated  in  the  absence  of  PRN  modulation,  which  was 
included  in  obtaining  Eq.  (F— 5)  (Appendix  F) , which  we  use  to  get  Eq. 

(18). 


The  scale  factor  degradation  (caused  by  single-bit  processing)  in 
the  output  of  the  delay  detector  in  the  coherent  DLL  is  given  by  the 
square  root  of  the  degradation  found  in  Section  5 for  the  noncoherent 
DLL.  Thus,  if  SNR^  * -41.3266  dB,  then  this  scale  factor  is  ^26490.8'. 
This  factor  is  inversely  porportional  to  yjmq  in  the  coherent  delay- 
locked  loop.  In  the  phase-locked  loop,  the  scale  factor  is  less  than 

this  by  a factor  of  2/tt , so  it  is  given  by  (2/tt)  ^ 2(^496. 81  for  SNRt 

« -41.3266  dB.  Again,  the  factor  is  inversely  proportional  to 
The  physical  cause  of  the  it/2  factor  in  mean  phase  detector  output  is 
that  the  recorded  samples  used  to  deteat  phase  shifts  in  a coherent 
phase  detector  are  the  samples  near  the  peaks  of  the  received  sinusoid, 
which  are  the  most  reliable  samples  in  a noisy  environment  (Appendix  D). 
The  samples  used  to  detect  delay  shifts  in  a delay  detector  are  dis- 
tributed among  reliable  and  unreliable  samples.  An  example  of  the  deter- 
mination of  the  performance  of  the  digital  receiver,  given  a set  .of  sam- 
ple signal  parameters  and  a sampling  rate,  is  done  in  Appendix  B. 


Ref,  5.  J.  C.  Springett  and  M.  K.  Simon,  "An  Analysis  of  the 
Phase  Coherent -Incoherent  Output  of  the  Bandpass  Limiter,"  IEEE  Trans. 
on  Communication  Technology.  Vol.  COM-19,  No.  1,  February  1971,  pp.  42- 
49. 
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I 

7.  LOOP  GAIN  CONTROL 


i 

1 

In  the  noncoherent  delay-locked  loop  and  the  suppressed-carrier 
tracking  loop,  the  scale  factors  affecting  the  detector  outputs  are  In- 
versely  proportional  to  input  SNR^.  In  the  coherent  delav-locked  loop 
and  the  phase-locked  loop,  the  scale  factors  affecting  the  detector  out- 
puts are  inversely  proportional  to  V§Nk^.  Thus,  if  it  were  desired  to  j 

have  constant  bandwidths  in  each  of  these  four  loops  it  would  be  neces*  ' 

sary  to  provide  a factor  proportional  to  1/(SNR^)  in  the  noncoherent 

loops,  and  proportional  to  1/ ^SNR ! in  the  coherent  loops.  A signal  j 

proportional  to  SNRi  is  found  at  node  ® in  Fig.  A-4.  This  signal  j 

may  be  inverted  to  yield  a number  proportional  tol/ (SNR)  for  use  xn  the 
noncoherent  loops.  A signal  proportional  to  V^^l  is  found  at  node  (a)  , 

in  Fig.  A-6.  This  signal  may  be  inverted  to  yield  a number  proportional  y ( 

to  l/y^SRJ  for  use  in  the  coherent  loops.  More  details  of  the  gain  con-  * 

trol  in  our  example  will  be  discussed  in  Appendix  B. 
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8.  EFFECTS  OF  ANTENNA  SWITCHING 


As  mentioned  in  Section  1,  the  missile  transponder  switches  its 
input  between  two  separate  pairs  of  receiving  antennas  every  1.152  ms. 
Synchronous  modulation  on  a pilot  carrier  signal  transmitted  by  the  mis- 
sile indicates  when  switching  occurs.  The  signal  from  each  pair  of  an- 
tennas will  be  tracked  using  its  own  receiver.  The  error  detectors  in 
the  various  component  loops  of  the  receiver  tracking  a specific  pair  of 
antennas  will  be  turned  off  while  the  other  pair  of  antennas  is  in  use. 

The  aiding  signals  will,  however,  continue  to  be  input  to  the  loops  dur- 
ing these  periods  of  no  signal  input. 

It  is  presumed  that  the  difference  between  the  actual  input  phase 
(or  delay)  and  the  aiding  signals  is  low  frequency  in  nature.  The  switch- 
ing rate  is  much  faster  than  any  of  the  receiver  loop  bandwidths,  tending 
to  minimize  the  dynamic  effects  of  the  switching  on  the  receiver  errors. 
Switching  transients  will  have  an  effect,  but  that  effect  will  not  be 
treated  here.  Under  these  circumstances,  the  effect  of  the  antenna 
switching  is  a 3 dB  loss  in  each  tracking  loop  (due  to  the  effective  de- 
crease of  integration  time). 


Thus,  with  2.5  MHz  sampling,  Eq.  (3)  is  replaced  by 


N (lO0,19991)  /N  (10°,19991)\2 

"h + (—2 

2A^(0.9029)  \ A*(0.9029)  / 


Equation  (11)  is  replaced  by 


» (LO0'24427) 
o 


2A2{E[Rp^(Ei)])2(0.9029) 


No(10°-2“27) 


(2B  ) 
w' 


5 L)1'2 

Nil! 

~)\  ' 


!2A^{E[Rpgj,  (Ej) ] }2(0.9029) 


) <2V  (' 


(20) 


)l/2 

(BN2tt  rad  ; 


(21) 
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Equation  (15)  is  replaced  by 


(22) 


Equation  (18)  la  replaced  by 


o [d>el (t> ] - VT" < 


(No(10 


0.24427 


) 


2A, 


{E[RPRN  <ei)n2(E(mm)]2 
i 


and  Eq.  (19)  is  replaced  by 


o(eai(t)] 


. J/v-pW i_ 

l\  2Ai  /\(E(coa[4el(t)]} 


)2  [E(raro) ] 


(24) 


The  results  of  rlght-clrcularly  polarized  shipboard  reception  and  left- 
circularly  polarised  shipboard  reception  of  the  input  signal  are  recorded 
separately  on  the  ship.  These  receptions  may  be  tracked  separately  by  the 
postflight  receiver.  The  result  is  a total  of  four  tracking  receivers  (2  an- 
tenna pairs  * 2 polarizations)  for  each  satellite  signal  transponded  by  a 
given  missile. 
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9.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 


In  this  report  the  performance  of  a single-bit  digital,  receiver 
for  carrier  and  code  processing  has  been  analyzed.  We  have  thus  demon- 
strated the  feasibil-  ':y  of  single-bit  quantization  in  the  recording  of 
spread-spectrum  signals  for  postf light  processing.  Diagrams  pertaining 
to  the  digital  receiver  were  given. 

The  following  facets  of  receiver  performance  were  not  discussed 
here,  and  are  suggested  for  future  investigations. 

1.  Acquisition  behavior:  The  digital  tracking  loops,  when  oper- 

ated in  quasilinear  regions,  are  expected  to  exhibit  acquisi- 
tion behavior  analogous  to  the  behavior  of  analog  loops. 
Nonetheless,  Investigation  here  would  be  useful. 

2.  Nonlinear  phenomena:  Nonlinear  performance  may  be  caused 

both  by  low  signal  levels  (threshold  effects  and  cycle- 
slippage)  and  by  high  signal  levels  (loss  of  dither  from 
input  noise) . These  phenomena  in  the  digital  loops  are  still 
to  be  Investigated. 
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Fig.  A-3  Noncoherent  Delay-Locked  Loop 


Output  (5) 
Output  (6) 


Lowpass  0 
[^Counter] 

Low  paw  (2) 
Counter  |-— 
Lowpats  (3) 

Lowpats  (4) 


Outputs  to  Suppressed- 
Carrier  Tracking  Loop 


1 Outputs  for  Gain  Control 
( Computation  (see  Fig.  A-4) 


Gain 

Control 

Factor 


THfc  JOHNS  HOPKINS  UNlVkHSlTf 

APPLIED  PHYSICS  LABORATORY 

laurel.  Maryland 


I 

I 25  Hz 


y 

i 


I I.dirm  .V:C*WI  V«  4MHfcMK»SW»K*^ 


d 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

uurel.  Maryland 


Appendix  B 

AN  EXAMPLE  OF  DIGITAL  RECEIVER  PERFORMANCE 


In  this  appendix  the  performance  of  the  digital  receiver  is  de- 
termined with  sample  values  of  the  relevant  parameters.  Let  the  signal 
power  be  -130  dBm,  the  noise  power  spectral  density  (PSD)  be 
N0/2  - 10-17.2  mW,  and  the  sampling  rate  be  5 MHz.  The  shipboard  re- 
ceiver bandwidth  is  still  taken  to  be  ±1.073  MHz.  The  following  analysis 
is  done  in  the  absence  of  antenna  switching.  (The  effect  of  antenna 
switching  is  given  in  Section  8 of  the  main  text.)  SNR^  is  given  by 


10 


-13.0 


N , 

(-“)  (2)  (2 .146  x 10b) 


(10“17 *2) (2) (2.146  x IQ6) 


-24.3266  dB  (B-l) 


Thus,  as  demonstrated  in  Eq.  (6),  the  scale  factor  affecting  the  mean 
output  of  the  delay  detector  in  the  noncoherent  delay-locked  loop  is 


(26490. 8)(102* 


43266  - 4. 13266 ^ 


528.561 


(B-2) 


The  scale  factor  in  the  coherent  delay-locked  loop  is  the  square  root 
of  this  quantity,  or  22.990.  The  scale  factor  in  the  phase-locked  loop 
is  (2/tt)(22.990)  - 14.636. 


by 
(71430) 


The  outputs  of  Lowpass  (J)  and  Lowpass  Q)  in  Fig.  A-3  are  given 

\ 1/2 


SNR. 


(l0-2. 43266)  ($28,561) 


\p(c  - 


+ "di(t) 


(71430) 


(SNR  \ 

-2  43266  ) mi  ^ ~ ^i^ ^PRN  ^i  ~ 

(10  c.AJZbb)  (528.56!)  / *LP  1 TRNi  11  d4 


+ n , (t)  , 

di 


(B-3) 
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where  71430  Is  the  maximum  possible  number  of  counts  at  the  outputs  of 
Lowpass  @ or  Lowpass  after  139.997“!  second.  The  variance  of 

ngi(t)  may  be  found  as  in  Eq.  (9) 


Signal  Power 
Noise  Density 


(■j)  (71430) 2 


SNKj 

(1Q-2^3266) 


(528.561)/ 


(139.997)'1c2[nal(t)] 


(10-°*16378)(i)A2 

N 

(f) 


(B-4) 


o2(na1<t)  J " °2  ] 

- 124326.8  counts2  . (B-5) 

Thus  na£(t)  has  a PSD  of  (124326.8) (139. 997)”1  - 888.07  (counts2/Hz) . 

The  variance  of  nai(t),  and  its  PSD,  are  independent  of  the  SNR,  as  can 
be  seen  from  Eq.  (B-4).  They  are  functions  only  of  the  sampling  rate 
and  the  accumulation  time  of  the  low-pass  filter.  We  have  assumed  the 
loop  sampling  rates  to  be  equal  to  the  shipboard  sampling  rate.  If  the 
sampling  rates  internal  to  the  loops  are  faster  than  the  shipboard  sam- 
pling rate,  the  signal  means  will  change  in  predictable  manners.  The 
signal  power/noise  density  ratios  will  remain  constant,  so  that  the 
noise  variances  may  again  be  determined,  using  Eq.  (9),  for  example,  in 
this  particular  case.  The  signal  at  node  (£)  in  Fig.  A-4  is  given  by 


8al(t)  - (71430)  I 2.43266 


(10 


- (71430)1 


— W 2 RJL  (A.  - A ) 

) (528.561 )/  lP  RNi  1 1 

) »4  “tL/*  - si + 


((10-2^3266)  (528.561) 


+ ngli(t)  ‘ 


(B— 6) 


AAA  J A 

For  » Aa,  RpRN^d^  - X1)  » 1,  and  - X±  + A)  w 0.  Also 

m2^  a 0.9029,  and  ngn(t)  has  PSD  given  by 
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8(888. 07)2(139. 997)  + 4(71430)2 


SNR. 


(l(f2*43266)  (528.561)' 


X ^RN  “ V<888-07>  counts2/Hz  • (B-7) 


Equation  (B-6)  may  be  rewritten  as 
2 / SNRi 

■al(t)  * <m30) 


K10 


) (528.561) 


) (0. 


9029)  + n (t)  , 


where  the  PSD  of  n ^^(t)  is  given  by 
8(888. 07)Z(139. 997)  + 4(71430)2  [ 


SNR. 


(10“2,43266)  (528.561), 


1(0.9029) (888.07) . 


Thus,  if  sai(t)  is  averaged  for  Ta  seconds,  its  mean  value  will  be 
given  by 

o / SNR.  \ 


(71430)' 


(10"2,43266)  (528.561) 


(0.9029) 


(B-8) 


(B-9) 


and  the  rms  noise  on  this  observation  will  be 


3(888. 07)2(139. 997)  +4(71430)2' 


SNR. 


ao"2'43266)(528.561), 


(0.9029)  ((*88.07) 


\ 1/2 

(yj  “um" 


(B-10) 

For  SNRj[  in  the  neighborhood  of  -24,3266  dB,  the  rms  error  (with  0.1 
second  averaging)  is  5.6  x 10^  counts.  This  gives  a relative  precision 
of  6%  in  the  estimation  of  SNR^  in  this  neighborhood;  that  is. 


5.6  x 105 


/ (71430)2(0,9Q29)  \ nn-2. 

\(10~2,43266)  (528.561 )/ 


tv  0.06 


-2.43266 


) 


In  place  of  that  part  of  the  gain  control  circuitry  in  Fig.  A-3  and 
Fig.  A-4  that  estimates  2a2[nai(t)],  it  is  possible  instead  to  use  the 
analytically  determined  value  in  Eq.  (B-5) . The  first  noise  term  in 
Eq.  (B-10)  would  then  be  reduced  by  a factor  of  2;  however,  in  this  case, 
if  the  loop  sampling  rate  is  variable,  it  would  then  be  necessary  to 
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average  cj  [nai(t)]  over  the  interval  of  interest,  as  the  sampling  rate 
changes.  Gain  control  proportional  to  (SNRj.)”1  is  necessary  in  the 
noncoherent  DLL,  in  the  suppressed-carrier  loop,  and  at  the  output  of 
the  suppressed-carrier  loop. 

Next,  the  performance  of  the  first  stage  of  the  digital  receiver 
will  be  determined.  Using  Eq.  (3),  we  find  that  (with  5 MHz  sampling, 
and  Bjji  ■ 1 Hz) : 


([(io-17-2)(io0-119A2- 
A I (l(f13,0)(0. 


9029) 


+ 4 P0~1^2^100,119-^-)2(140)  ID./72 

\ (10“1J,U)(0 .9029)  / ) 

o(e.)  , 

— - (6.955  x lo”J) A - 6.799  ns  . 

Then,  from  Eq.  (A) 

E[RPRN1(ei)1  " 1 -(“=:)  (6‘955  X 10~3)  " °* 

Further,  from  Eq.  (11)  we  have  (let  • 1 Hz) 


99AA5 


'll 


(10~17-2)(10°-16378) 
(10”13, °) (0 . 99AA5) 2 (0 . 9029) 

; (io-17«2)(io0-16378) 

i(10~13*°) (0.99A45)2(0.9029) 


j (1A0)  (l) 


t (n~. (t)]  - 0.01022  x 10"J  rad  - 0.5858*  , 

Jl  rms 


(B-ll) 


(B-12) 


(B-13) 


and  obtain  from  Eq.  (12) 


E(cos[n31(t)]}  - exp  q. 


,999948  . 


(B-14) 
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For  the  operation  of  the  bit  estimator,  we  obtain  from  Eq.  (15) 

» 

R'  -( ITfY  — o 16378  V20  * 10"3) 

bi  \ (2  x 10  17 * 2) <10° • 16378>/ 

R^1  - 108.698  ; 
then  from  Eq.  (14) , 

P (E)  w J 1 - erf  ( V 108.^98)  (0. 99AA5)  (0.9999A8)]  1 


P (E)  |l  - erf  (10.3674)] 


P(£)  w 5.6744  x 10 


and  from  Eq.  (16) 


P(.  i i)  ~ (5.6744  x 10-49) 

\ 20  x 10"-5 


Hi 


.799  x 10 
20  x IQ'3 


P(m  * A)  * 1.69975  x 10-'  . (B-17 

Finally,  from  Eq.  (17), 

E(mA)  - (1  - 3.3995  x 10"7)  *»  1 

E(md)  * 1 . (B-l£ 

This  value  of  E(mA)  indicates  that  the  data  estimation  is  essentially 
perfect,  so  far  as  the  second  stage  of  the  receiver  is  concerned.  For 
smaller  signal  power/noise  density  ratios  E(m&)  is  smaller. 

Next,  the  operation  of  the  second  stage  of  the  receiver  is  dis- 
cussed. As  mentioned  in  Section  8 of  the  main  report,  the  gain  control 

signal  for  the  coherent  loops  is  found  at  node  (a)  in  Fig.  A-6.  If  the 

cycle  time  of  the  counter  at  this  node  is  given  by  0.1  second,  then  the 
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maximum  possible  number  of  counts  observable  there  is 
(0.1) (5  x 106) (2)  - 106  . 

The  signal  at  node  @ is  analogous  to  that  at  the  output  of  Lowpass 
(l)  in  Fig.  A-3,  as  described  in  Eq.  (B-3) : 

/ SNR  \1/2 

dO6)  ( TTZiZiT 1 mA  (t  “ x<>®<  (t  " x^rprm  <xi  " X4.>C08 

\a0"2-43266)  (528.561)/  *LP  1 4LP  1 1 iA  ei 


+ ng2i(° 


(B-19) 


The  signal  pover/noise  density  ratio  is  the  same  as  in  Eq.  (B-4) . We  may 
thus  determine  a2[ng2i(t)]. 


Signal  Power  m 
Noise  Density 


1 6 / SNRi  \ 

(?,<10>(a02-43266)(«77.2388)) 


(0.1)o  (ng21(t)) 

(10"°*11942)(i)(A2) 

N 

<f> 


(B-20) 


o2[ng21(t)]  - 1740541  (B-21) 

For  Ai(t)  «A^(t)  and  $ei(t)  *»  0,  we  have  that  RpRNi(Aj  - A^a)  w 1,  and 

cos  $e£(t)  1.  The  value  of  miLp(t  - Ai)miLp(t  - A^)  is  itself  a weak 

function  of  SNRi,  * f*ct  that  we  will  Ignore  and  merely  assign  it  the 
value  1,  as  determined  earlier  in  this  appendix.  For  SNR  in  the  neigh- 
borhood of  -24.3266,  and  with  0.1-sacond  averaging,  the  rms  noise  on  the 

measurement  is  \174054T1  ■ 1319  counts,  and  the  relative  precision  of  the 
estimate  of  (SNRi)*'^  is  3%. 

The  performance  of  the  second  stage  of  the  receiver  is  determined 
next.  We  first  rewrite  Eqs.  (18)  and  (19),  the  simultaneous  equations 
for  the  rms  errors  of  the  second-stage  tracking  loops.  Recall  that  the 
5 MHz  sampling  rate  changes  a factor  in  these  equations,  as  given  in 
Figs.  1 and  2. 
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We  will  solve  these  equations 


°[*el(0]  “ 


o[eal(t)] 


N (10 
o 


0.16378 


No<10 


2Ai 

0.09427 


!)( 


1/2 


(E{RpRN[eai(t)]»2(E[mm])2 


2A‘ 


-)( 


(E{ccn[*el(t)]))2(E(n»])2 


)M 

xw 


rad 

(B-22) 

1/2 

(B-23) 


We  will  solve  these  equations  Iteratively.  Start  by  setting 
ElRpRN^Ctal^))  1 " 1 in  Eq.  (B— 22) . Then,  using  Eq.  (B-18)  and  letting 

Bjq^  * 1 Hz  we  obtain 


oUgl(t)] 


(10~17,2)(100,16378) 


10 


-13.0 


1 1/2 


1 1 


rad  , 


al*el(t)]  - 9.592  x io"3  rad  - 0.5755*  . 

Substituting  Eq.  (B-24)  into  Eq.  (12),  we  obtain 
E(cos  n31(t)]  - 0.999954  . 


(B-24) 


(B— 25) 


Substituting  Eq.  (B-25)  into  Eq.  (B-23),  and  letting  BN3  « 1Hz,  we  have 


o(eai<t)l 


o[eai<t)l 


'(io-17-2)(io°-09427) 


- 6.261  x 10 


10 

-3 


-13.0 


■)  / i )(o.s) 

/ \(0, 999934) (1)/  \ / 


1/2 


(B-26) 


We  then  use  this  quantity  in  Eq.  (4)  to  obtain 
E{RPRNi[eai(t)1}  “ °* 99500  * 

By  substituting  Eq.  (B— 27)  into  Eq.  (B-22),  we  now  obtain 
o[*ei(t)]  - 9.640  x io’"3  - 0.5523°  . 


(B— 27 ) 


(B-28) 
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The  iteration  is  continued  until  o[ea^(t)]  and  c[$ei.(t)]  converge,  which 
has  already  happened  here,  because  substituting  Eq.  (B-28)  into  Eq.  (12) 
yields  E(cos  ^3i(t) ]■  0. 999954,  which  is  the  value  in  Eq.  (B-25).  Thus, 

a[*el(t)]  - 0.5523°  (B-29) 

a(cai(t)]  “ 6.261  x 10“3d  - 6.12  ns  . (B-30) 

These  quantities  characterize  the  performance  of  the  digital  receiver, 
given  the  signal  parameters,  sampling  rate,  and  receiver  parameters 
assumed  in  this  appendix. 

This  performance  characterization  accounts  for  errors  caused  by 
input  noise  only.  As  mentioned  in.  Section  A of  the  main  report,  high 
signal  levels  can  attenuate  the  effects  of  input  noise  to  the  order  of 
the  effects  of  noises  generated  internally  by  the  tracking  loops.  It  is 
shown  in  Appendix  D that  the  noise  generated  internally  by  the  phase- 
locked  loop  is  approximately  0.6°  rras,  independent  of  the  signal  level. 
Thus  the  rms  noise  at  the  output  of  the  phase-locked  loop  is  in  fact 
given  by 

o[$ei(t)]  - V 0.5523^  + 0.6*'-  0.815°  (B-31) 

A similar  result  would  apply  to  the  delay-locked  loop.  The  internal 
noise  levels  in  the  loops  depend  on  the  specific  hardware  implementa- 
tions, as  discussed  in  Appendix  D. 


AA  - 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL.  MARYLAND 


Appendix  C 

SATRACK  CARRIER  AND  PRN  PROCESSING  USING 
AN  ANALOG  RECEIVER  MODEL 


Presented  here  is  a receiver  design  for  SATRACK  carrier  and  code 
processing.  The  discussion  is  in  terns  of  an  analog  input  signal  (not 
sampled,  not  quantized).  In  the  SATRACK  scheme  a sum  of  carrier  signals 
is  transponded  by  the  missile  from  the  satellites  to  the  range  ship. 
These  carriers  are  suppressed  by  pseudorandom  noise  (PRN)  biphase  modu- 
lation (±1  amplitude  modulation,  with  a different  sequence  for  each  sat- 
ellite), and  ±1  data  modulation.  In  the  postflight  receiver  it  is  de- 
sired to  track  the  carrier  doppler  shifts  and  the  different  delays  of 
the  PRN  sequences.  Telemetry  and  surface  station  data  from  the  flight 
is  provided  to  the  postflight  receiver  to  aid  in  this  task.  This  appen- 
dix discusses  a proposed  design  for  a postflight  receiver  and  analyzes 
its  performance. 


A RECEIVER  DESIGN  FOR  SATRACK  CARRIER  AND  CODE  PROCESSING 

The  signal  transponded  by  the  missile  from  the  N in-view  satel- 
lites to  the  range  ship  in  the  SATRACK  scheme  is  received  by  the  ship  in 
the  form: 


x(t) 


PRN^(t  - Xi)m1(t  - X^)cos[wot 


+ VO]  + n(t)  . (C-l) 

Here,  \±  denotes  the  total  time  delay  from  transmission  of  the  signal  by 
satellite  (i)  to  its  final  reception  by  the  ship.  The  function  PRN^(t) 
is  a pseudorandum  sequence  of  ±l’s  with  a bit  (chip)  length  of 
A « 10-6/1,023  second  and  a sequence  repetition  rate  (epoch  rate)  of  1 
kHz  (i.e,,  1023  chips  per  epoch).  The  function  m^(t)  is  also  a sequence 
of  ±l’s,  which  carries  satellite  orbit  data,  among  other  things.  The 
chip  length  of  mj(t)  is  A^,  - 20  ms.  The  phase  function  $i(t)  is  the  re- 
sult of  the  doppler  shifts  caused  by  satellite,  missile,  and  range  ship 
movement.  The  noise  function  n(t)  is  white  noise,  which  is  bandwidth 
limited,  by  the  shipboard  receiver,  to  ±1  MHz  about  the  missile  trans- 
mission frequency  of  f0  ■ 2200  MHz.  The  power  spectral  density  (PSD)  of 
n(t)  is  given  by  (Nc/2)  ■ 10"^7.2  row/Hz.  The  normal  carrier  power  for 
the  ith  signal,  P8^,  is  given  by  10"^  • 7 mw,  although  a wide  variation 
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may  be  observed  in  the  signal  power  received  from  the  different  satel- 
lites. The  nominal  signal  power/noise  density  ratio  is  thus  given  by 
2P„-|/N0  * 25  dB*Hz.  The  code  sequences,  PRNi(t),  and  their  epoch  times 
(as  transmitted  by  the  satellites)  are  known  by  the  postflight  re- 
ceiver, although  the  delays  (X^)  are  not  known  a priori.  The  doppler 
shifts  {$i(t)>  are  also  not  known  apriori.  However,  telemetry  and  ground 
station  data  are  available  to  the  receiver,  and  are  converted  into  aiding 
Information  for  the  various  receiver  tracking  loops.  The  data  m^(t)  may 
be  partially  known  apriori,  but  will  be  treated  here  as  totally  unknown, 
causing  a very  small  deterioration  in  code  and  carrier  tracking  perfor- 
mance. 


The  received  signal  x(t)  is  heterodyned  by  2200  MHz  to  a nominal 
zero  center  frequency  by  the  shipboard  preprocessor  and  then  hard  lim- 
ited, sampled,  and  recorded  for  later  processing  by  the  postf light  re- 
ceiver (Appendix  D),  Here  the  single-bit.  quantization  (hard  limiting) 
and  sampling  of  x(t)  is  Ignored,  and  an  analog  receiver  that  will  per- 
form the  objectives  stated  above  is  discussed.  The  performance  of  a 
specific  design  for  this  receiver  will  be  analyzed,  in  the  expectation 
that  the  performance  of  a digital  receiver  operating  on  the  heterodyned, 
single-bit  quantized,  sampled  version  of  x(t)  will  be  much  like  the  per- 
formance of  the  analog  receiver.  This  expectation  is  not  unfounded,  be- 
cause in  Appendix  D it  is  concluded  that  one  component  of  this  system, 
a single-bit  quantized  digital  phase-locked  loop,  will  perform  in  the 
same  manner  as  an  analog  phase-locked  loop,  suffering  only  a 1.15  dB 
degradation  in  output  phase  variance. 

Each  of  the  N delays  {X^}  and  doppler  shifts  { (t ) > are  tracked 
with  a separate  receiver.  The  ith  receiver  uses  its  knowledge  of  the 
ith  code  sequence  PRNi(t)  to  separate  its  signal  component  from  the 
other  N-l  signal  components  in  x(t)  (see  Eq.  (C-l)).  The  Interference 
in  the  ith  receiver  from  the  other  N-l  signal  components  in  x(t)  will 
affect  the  performance  of  the  receiver,  and  will  superpose  its  effect 
with  that  caused  by  the  noise  n(t) . 

This  appendix  is  limited  to  analyzing  the  effects  of  the  noise 
n(t)  on  the  performance  of  the  ith  receiver.  Thus  the  signal  x(t),  which 
is  input  to  the  ith  receiver,  will  be  dealt  with  as  though  it  were  in 
fact  x^(t),  a fictitious  signal  defined  here  as: 


x^t)  PRNi(t  - A1)m1(t  - X1)cos(o)ot  + ^(t)]  + n(t) 

(C-2) 
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After  x(t)  is  processed  through  the  N receivers  it  might  be  possible  to 
make  a second  pass  through  x(t)  with  each  receiver  having  its  input, 
x(t),  altered  to  remove  the  interfering  signal  components,  using  infor- 
mation supplied  by  the  N-l  other  receivers  from  the  first  pass  through 
x(t). 


The  ith  analog  postf light  receiver  has  the  form  shown  in  Fig.  C-l. 
The  first  stage  of  the  receiver  estimates  the  function  m^(t  ~ X^) , and 
stores  it  on  tape.  The  second  stage  of  the  receiver  uses  this  tape  to 
aid  in  estimating  $i(t)  and  the  delay,  Xj,,  in  the  reception  of  PRN^t). 
The  computer  aids  obtained  from  telemetry  and  ground  station  data  are 
used  in  both  stages  of  the  receiver.  Detailed  diagrams  of  the  various 
components  in  Fig.  C-l  will  be  given  along  with  the  discussion  of  re- 
ceiver performance. 


ESTIMATION  OF  rn^t  - X^ 

The  first  processing  performed  by  the  ith  receiver  on  the  input 
signal  x(t)  is  done  with  a noncoherent  delay-locked  loop  (Ref.  C-l),  as 
pictured  in  Fig,  C-2.  (Recall  that  we  are  working  with  x^(t)  rather 
than  x(t)  in  our  analysis  of  the  effects  of  noise  on  system  performance.) 
This  loop  locks  onto  the  delay  X^  in  the  code  PRNj:(t  - X^) , using  a 
feedback  shift  register  (FSR)  code  generator  that  reproduces  the  pseudo- 
random sequence  PRN^.  The  operation  of  the  loop  is  noncoherent  with  re- 
spect to  both  the  input  carrier  phase  $i(t)  and  the  data  mi(t  - , 

both  of  which  are  unknown  at  this  stage  in  the  processing.  Define  the 
error  in  delay  tracking  as 


A 

where  X*  is  the  estimate  of  delay  X^  generated  by  the  loop.  It  is  shown 
in  Ref.  C-l  that  (see  Fig,  C-2) i 

e^)  - ; |etl  *A/2  . (C-3) 

We  choose  for  the  closed  loop  a transfer  function  that  is  optimum 
in  terms  of  total  mean-squared  delay  error  for  ramp  inputs  of  delay  in 
the  presence  of  noise.  This  transfer  function  is  given  by 


Ref.  "c-l.  W.  J.  Gill,  "A  Comparison  of  Binary  Delay-Lock  Tracking 
Loop  Implementations,"  IEEE  Transactions  on  Aerospace  and  Electronic 
Systems,  Vol.  AES-2,  No.  4,  July  1966,  pp.  415-424. 
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*i(,)  a „ . . (l  + V?b/po) 

xl<*>  ‘ 1 ' II  ♦ V*./,.  + Cs/P0)2] 


(C-4) 


and  the  two-sided  closed  loop  noise  bandwidth  is  given  by 
®N1  " ^*0®  P0  Hz  . 


(C-5) 


The  transfer  function  in  Eq.  (C-4)  can  be  approximated  by  choosing  as 
the  transfer  function  for  the  loop  filter 


(1  + Jz'b/p  ) 

Fi('>  • -<r+t»7»;r 


(C-6) 


where  g is  a normalized  gain  constant  defined  by 


g 


F.l*f8c 


t 


(C-7) 


and  g is  chosen  to  be  large  by  making  gf  large.  It  is  then  shown  in 
Ref.  C— 1 that  the  rma  value  of  the  normalized  error  of  the  loop  delay 
estimate,  e.(t)/A  , is  given  by 


f - <Ei'« 


rms 


(C-8) 


where  B-^  is  the  bandwidth  of  the  bandpass  filters  preceding  the  square 
law  devices  in  Fig.  C-2.  A lower  limit  of  0.1  Hz  set  on  loop-noise 
bandwidth  for  any  of  the  loops  to  be  discussed.  For  lesser  bandwidths 
there  is  danger  of  having  noises  generated  within  the  loops  (which  will 
not  be  accounted  for  in  the  present  analysis)  dominating  the  error  per- 
formance, It  is  obvious  from  Eq.  (C-8)  that  the  smaller  is  made, 
the  better  will  be  the  loop  performance.  In  fact,  considering  a value 
of  oe  » 0.3A  to  be  the  threshold  value  for  proper  loop  operation  (Ref. 
C-2)  we  find  that  it  is  necessary  to  keep  low  in  order  to  guarantee 
proper  performance  and  the  maintenance  of  loop  delay  lock.  The  bandpass 
filters  must  be  wide  enough  to  pass  mi(t  - X^)  relatively  unchanged,  so 
that  mi(t  - X-[)  disappears  in  the  square  law  devices.  Also,  their  pass- 
band  must  be  wide  enough  to  always  encompass  the  output  frequency  of  the 


Ref.  C-2.  J.  J.  Spilker,  Jr.,  '*Delay  Lock  Tracking  of  Binary  Sig- 
nals,"  IEEE  Trans,  on  Space  Electronics  and  Telemetry.  Vol.  SET-9,  March 
1963,  pp.  1-8. 
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mixers,  which  varies  with  the  doppler  shift  of  the  input  signal.  The 
ra^Ct  ~ Xi)  requires  a bandpass  of  Bj.  * 100  Hz.  The  input  doppler  shift 
can  be  as  much  as  ±100  kHz.  Using  Eq.  (8),  given  the  nominal  SATRACK 
signal  and  noise  levels  and  the  lower  limit  BN1  - 0.1  Hz,  the  maximum 
allowable  value  of  for  oe  ■ 0.3A  may  be  determined,  which  is 
Bi may  * 44  kHz . This  bandwidth  will  not  accommodate  the  doppler  shift. 
Therefore  aiding  information  must  be  provided  to  force  the  local  oscil- 
lator to  follow  the  doppler  shift  of  the  input  signal  to  some  degree. 

There  are  a number  of  alternative  methods  to  provide  the  aiding. 
For  the  present  discussion  the  telemetry  and  ground  station  data  that 
are  provided  to  the  postf light  processor  are  utilized.  These  data 
force  the  local  oscillator  in  the  envelope  correlator  to  follow  the  in- 
put doppler  frequency  with  a nominal  offset  of  fixed  frequency  f^»  and 
a phase  error  of 

♦d  <t)  " *i(t)  " ^ref/0  » (C”9) 


where  ♦ref i^t)  is  the  computer  phase  aiding,  which  is  of  sufficient 

accuracy  to  ensure  that  $d^(t)  « 2 it  • 100  rad/s.  Thus,  B^  is  set  equal 

to  100  Hz,  the  value  required  to  pass  m^(t  - X j_) . Then  Bjji  is  arbitrar- 
ily chosen  as  2 Hz,  such  that  from  Eq.  (C-8)  (Je  * 0.08463A,  well  within 
threshold.  Thus  the  code  epoch  estimates  are  accurate  to  an  rms  error 
of  at  « Q.08463A  - 82.73  ns.  The  output  of  the  noncoherent  delay- locked 
loop  is  given,  Ref.  C-l,  by 


yU<t>  ■ ^(t  - xi)RpRNt«:i(t)]costw1t  - ^ (t)]  + n^t)  , 

where 

RpRN(e)  M 1 - ^ 5 M < A , 


(C-10) 


(C-U) 


Ref.  C-3.  W.  C.  Lindsey  and  M.  K.  Simon,  Telecommunication  Systems 
Engineering,  Prentice-Hall,  Inc.,  New  Jersey,  1973. 
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and  n^(t)  Is  white  noise  of  spectral  height  N0/4,  which  is  band-limited 
to  - (Bjl/2)  < | f | & + (B]_/2)  . By  making  the  approximation  that 

ei(t)  is  a Gaussian  random  process,  it  is  found  that 

E[RPRN(ei(t)) 3 “ 0,93247  (C“12) 

so  that  Eq.  (C-10)  may  be  replaced  with 

yii(t)  “ V^si  “i^  " Aj.)  (0.93247)cos[w1t  - $d  (t)  ] 


+ n^t)  , 


(C-13) 


where  Eq.  (C-13)  in  the  following  discussion  is  justified  whenever  the 
noise  bandwidths  of  the  various  loops  are  much  narrower  than  the  fre- 
quency content  of  R_n„[£4 (t) ] . 

The  next  step  in  obtaining  an  estimate  of  m^Ct  - X^)  is  to  re- 
move the  cosine  in  Eq.  (C-13)  from  yj^(t).  This  factor  is  effectively 
a carrier  wave,  which  is  suppressed  by  m±(t  - X*) . Thus  it  is  necessary 
to  employ  a suppressed  carrier  tracking  loop  to  recover  this  cosine. 
Assume  that  a squaring  loop  (which  gives  equivalent  performance  to  a 
Costas  loop)  Is  used.  A decision-feedback  loop  could  be  used  here,  with 
slight  Improvement  in  overall  system  performance,  but  the  loop  acquisi- 
tion problem  would  then  be  complicated  by  the  initial  lack  of  data  clock 
synchronization,  which  is  required  for  decision  feedback.  A diagram  of 
a squaring  loop,  Ref.  C-3,  is  given  in  Fig.  C-3.  The  bandpass  filter, 
D2(s),  must  be  wide  enough  to  pass  both  m^(t  - X*)  and  the  cosine  (l.e., 
100  Hz  as  before).  The  second  bandpass  filter,  03(e),  is  centered  at 
frequency  2fi,  and  must  be  wide  enough  to  pass  the  component  of  y^(t) 
at  2fj_.  A 100  Hz  bandwidth  is  assumed  sufficient  here  also.  The  output 
of  the  squaring  loop  and  frequency  divider  is 


y21(t> 


acostu^t  - $d  (t)  + n3(t)]  , 


(C-14) 


where  03(0  is  the  phase  tracking  error.  Using  the  discussion  of  sup* 
pressed  carrier  tracking  loops  in  this  appendix,  it  is  found  that  the 
rms  value  of  03(f)  becomes 

|l/2 


‘V^rmo 


N 


N 


2P  * (0.93247)' 

V 


2P  • (0.93247)' 


\ 2 

/ 

) 
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where  BN2  is  the  loop  noise  bandwidth  of  the  squaring  loop,  which  is 
assumed  to  be  1 Hz.  We  obtain 


[n3(t)]rros  - 0.07042  rad  - 4.0350°  (C-16) 

As  Indicated  in  Fig.  C-l,  yii(t)  is  multiplied  by  y2i(t)  to  ob- 
tain 

y31(t)  - mi(t-Xj)  [cos  n3(t)]RpRN[e1(t)]  + n^(t)  , (C-17) 


where  n^(t)  is  white  noise  of  spectral  density  N0/2,  bandlimited  to 
|f|  S (Bi/2),  (Bi  - 100  Hz). 

\ 

The  signal  y3i(t)  is  next  put  into  the  data  clock  synchronizer. 
The  data  m^t),  when  modulated  onto  the  signal  x^(t)  before  transmission 
by  the  satellite,  has  its  transitions  at  the  exact  time  of  PRNi(t) 
epochs.  Because  the  chip  length  of  mi(t)  is  exactly  20  times  the  epoch 
length  of  PRNi(t),  which  one  of  20  possible  PRN  epoch 'times  is  aligned 
with  the  m(t)  chip  must  be  determined.  In  order  to  do  this,  the  follow- 
ing test  is  performed.  The  signal  y3i(t)  is  integrated  over  each  PRN 
epoch,  as  determined  in  the  noncoherent  delay-lock  loop,  to  obtain 

(j+1)  x 10~3 

/ y3i(t)dt  ; j - 0,  ...  , (Tp  x io3  - 1)  , 

.1  x 10"3  (C-18) 

where  T is  the  length  in  seconds  of  the  combined  satellite-missile 
pass  (w  300  s).  Thus  there, are  Tp  x 10^  separate  intervals  over  which 
y3i(t)  is  integrated.  The  values  of  (y^|'}  may  Be  stored  on  tape,  to 

avoid  recomputing  them.  The  alignment  of  mi(t)  with  PRNj,(t)  may  occur 
at  any  of  the  first  20  epochs  of  PRNi(t).  If  we  are  testing  for  whether 
the  ith  epoch  is  correct,  the  following  summation  must  be  computed: 


t-i 


V10 

(-"20'— > “ 2 20(k+l)+l-l 


••■I IE  -i.il  * |E  IE -i.1 

J-0  k-0  ,1»7.0k+t 


T «10  -1 
P 


* llE  -5‘  I 


J-Tp*l<r-20+l 


(C-19) 
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Whichever  value  of  l yields  the  highest  value  of  g£  is  chosen  for  the 
alignment  of  m^(t  - A^)  with  the  epochs  of  PRN^(t  - A^) . The  section 
of  this  appendix  on  data  clock  synchronization  shows  that  the  probability 
of  choosing  the  wrong  value  of  i for  the  synchronization  of  the  data  bits 
with  the  code  epochs  is  low  enough  that  computing  time  may  be  saved  by 
using  only  ISO  seconds  of  data  (Instead  of  Tp  w 300  s)  in  the  synchronize 
tion  process,  while  still  having  a probability  of  synch  error  on  the 
order  of  10“°,  which  can  be  ignored.  Of  course,  the  estimates  of  the 
code  epoch  are  in  error  by  an  rms  value  of  82.73  ns,  so  that  the  data 
bits  will  also  be  off  in  alignment  by  this  amount.  The  integration  of 
y3i(t)  over  each  data  chip  interval  determines  whether  m^(t  - Ai)  is 
positive  or  negative  over  that  chip.  These  integrals  are  given  by 

20(k+l)+£-l 

bk  - yJ31  (c-20) 

i«20k+£ 

The  probability  of  incorectly  determining  - A^)  through  integration 

over  a chip  interval  (Am  - 20  ms) , if  the  sta  .'ting  time  of  the  interval 
were  exactly  known,  is  given  from  Eq.  (C-17),  using  Ref.  C-3,  pp.  305- 
319. 


Prob  [m^t  - A^)  i m^t  - A^)  ] 

\ |l  “ erf  J^R^,RpRNte1(t)]co8  t»3(t)J 


- E 


(C-21) 


where  it  is  assumed  that  e.(t)  and  n~(t)  are  relatively  constant  for 
20  ms,  and 


.PA  P x 20  x 10 
„ A 8 m 8 

Eb-  

o o 


(C-22) 


However,  because  there  is  an  rms  offset  of  82.73  ns  in  the  alignment 
time,  R^  effectively  becomes, 


P (20  x 10"3  - 82.73  x lo”9) 
8 

N 


- 3.1622 


(C-23) 
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Next  the  quantity  In  Eq.  (021)  will  be  evaluated.  The  random  variable 
g ie  defined  as 

g - RpRN[c1(t)]coa  n3(t)  (024) 

To  find  the  mean  value  of  g it  is  assumed  that  e^(t)  and  n3(t)  are  inde- 
pendent. Using  Eq.  (11)  and  assuming  that  e(t)  is  Gaussian, 

E*Wei<t)]  “ 0,93247  (C"25) 

has  already  been  obtained.  Assuming  n*(t)  is  Gaussian,  and  using  the 
definition  of  the  characteristic  function  of  a Gaussian  random  variable 
n3  * N(0,o2) , 

♦ <n3>  £ Ee>n3  - «p(-^)  (C-26) 

and 

E cos  n3(t)  - 0.99752  . (027) 

Thus, 

Eg  - ERpRN[ei(t)  ]E  cos  n3(t) 

Eg  - 1.6541  (C-28) 

Because  erf(z)  is  quite  linear  near  z « 1.6541,  and  g is  nearly  evenly 
distributed  in  probability  about  Eg,  the  approximation  is  made  that 

E[erf  (g)  ] - erf  (Eg)  (029) 

Hence,  wc  find  from  Eqs.  (021) , (024) , (C-28),  and  (029)  that 

Pl«i(t  - Xt)  4 mA(t  - Xj)]  - y[l  - erf (1.6541) ] 

- 9.6632  x 10”3  (C-30) 

The  quantity  E[a<(t  - XA)m(t  - X*)],  which  will  be  used  later,  is  eval- 
uated using  Eq.  (030). 
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E[o1(t  - A^m^t  - Ai)]  - 1 x Prob (m  - m)”1  x Prob(ra  4 m) 


1 • f (1  - 9.6632  « IQ’3)  (j^lO-3  - 82.73  . 10'9 


20  x 10 


-3 


1 / 82.73  x 10~9\ 
* \ 20  x 10”3  / 

[<• 


20  x 10 

1.6632  x 10~O 

l-9 


-3,  (20  x 10"3  - 82.73  x 1Q~9^ 


20  x 10 


-3 


1 / 82.73  x lQ-*\ 

2 \ 20  x l(f  3 / 


(t  - A1)rai(t  - A^]  - 0.98067 


(C-31) 


COHERENT  TRACKING  OF  ^(t)  AND  X± 

The  operation  of  the  phase-locked  loop/coherent  delay-locked 
loop  (PLL/DLL),  which  comprises  stage  2 of  the  postf light  receiver  for 
xi(t),  will  be  aided  by  the  value  of  mi(t  - Ai)  determined  by  stage  1, 
the  code  epoch  synchronisation  performed  in  stage  1 (which  will  aid  in 
the  code  acquisition  of  the  delay-locked  loop) , and  the  telemetry  and 
ground  station  data  available  to  the  postflight  receiver.  The  PLL/DLL 
loop  Is  shown  In  Fig.  C-4.  The  errors  in  this  loop  will  be  determined 
Iteratively.  If  given  that  the  phase-locked  loop  tracks  with  error 

<fe  (t)  “ “ ♦j(t)*  then  the  normalized  rms  value  of  the  delay  error 

e#(t)  ■ A^(t)  - A^a(t)  of  the  the  delay-locked  loop  is  given  (Fig.  C!-4 

and  Ref.  C-l)  by 


BN3No 


fT 


4P  [E  cos  4e(t)]  E[m(t  - Aj)m(t  A^) ] 


(C-32) 


where  B^n  is  the  loop  noise  bandwidth  of  the  delay-locked  loop.  We  use 
expectations  in  Eq.  (C-32)  under  the  assumption  that  Bn3  is  narrower  than 
the  frequency  content  of  ♦•(t)  and  m(t  - A^)m(t  - A^),  Also  assume  that 
4e(t)  and  m(t  - A^)  are  independent.  Further,  if  given  a delay  error 


\ 
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ea(t)  in  the  output  of  the  delay- locked  loop,  then  we  obtain  from  phase- 
locked  loop  theory,  Ref.  C-3,  p.  29,  that  the  rms  phase  error  at  the  out- 
put of  the  phase-locked  loop  is  given  approximately  by 


a . r_ BN4^o , . , . 

^ ■ 4/  i pjj  ~ 2 rac*  >(C-33) 

1 2Pg  | ERpRN[efi(t)]  j 1 [Em(t  - X^mCt  - X^P 

where  Bjj^  is  the  loop  noise  bandwidth  of  the  phase-locked  loop.  It  is 
assumed  that  Bjj3  » 0.2  Hz,  B^  ■ 1 Hz,  and  will  now  determine  and 

o^e  iteratively.  First  let  RpRij[ea(t)]  " 1 in  Eq.  (C-33) . Then  obtain 

o^  « 0.057342  rad  and  approximate  $e(t)  as  Gaussian.  Then 

E cos  $e(t)  " 0.99836.  Using  this  and  Eq.  (C— 31)  in  Eq.  (C— 32)  there 
obtains  oGa  ■ 0 . 01&163& . This  quantity  is  then  put  back  into  Eq.  (C-33). 

Approximating  ea(t)  as  Gaussian,  and  using  Eq.  (C-ll) 

EtRpRjj  [ea(t)  ] } “ 0.98551.  Substituting  this  value  into  Eq.  (C-33), 
o(j>e  ■ 0.05816  rad.  The  iteration  between  Eq.  (C— 32)  and  Eq,  (C-33)  is 

reiterated  until  estimates  of  0^  and  oe  converge  (to  five  significant 

digits),  which  happens  on  the  next  iteration,  Then 

o - 0.018164A  - 17.76  ns  (C-34) 

ea 

and 


a - 0.058186  rad  - 3.334°  . (C-35) 

^e 

These  answers  then  are  the  rms  values  of  the  errors  in  the  receiver  esti- 
mates of  4u(t)  and  X^,  Note  that  only  errors  caused  by  input  noise  have 
been  considered,  so  that  both  noise  generated  internally  by  the  various 
loops  and  errors  caused  by  signal  dynamics  have  been  neglected. 


SUPPRESSED  CARRIER  TRACKING  LOOPS 

In  this  section  three  types  of  suppressed  carrier  tracking  loops 
are  discussed:  the  squaring  loop,  the  Costas  loop,  and  the  decision- 

feedback  loop.  The  analysis  is  based  on  linearization  of  the  equations 
given  in  Ref.  C-3.  Comparison  will  be  made  of  the  operation  of  these 
three  loops  when  the  following  input  signal  is  applied  to  them. 
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y(t)  ■ J2?  d(t)  sin(ur  t + 8(t) ] + n (t)  , (036) 

^8  0 W 

where  d(t)  Is  modulation  in  the  form  of  a ±1  pulse  train  with  bit  inter- 
val  I'd,  and  nw(t)  is  white  noise  of  spectral  density  (N0/2). 

A diagram  of  a squaring  loop  is  given  in  Fig.  03.  The  first 
bandpass  filter,  D2 (s)  (bandwidth  B2),  changes  nw(t)  from  white  noise  to 
bandlimlted  white  noise,  nw].(t),  over  frequencies 

fo  - B2/2  * jf | * f0  + B2/2  , 

while  passing  d(t)  relatively  unchanged.  Thus  the  output  of  the  filter 
D2<s)  may  be  written 

z (t)  - ^2P^d(t)  sin[toot  + 6(t)]  + nwl(t)  , (037) 

and  nwl(t)  may  be  written 

«*Ml<t)  - cos  wQt  - nfl(t)  sin  u>ot]  , (038) 


where  n?(t)  and  n8(t)  are  Independent  white  noise  processes  of  power 
spectral  density  N0/2,  bandlimlted  to  |f|  * (B2/2).  Alternatively, 
n^(t)  may  be  written 

nwi(t)  - ^5  } Nc ( t)  cornet  + 8(t)l  - Nfl(t)  sin[wot  + 8(t)]j  , 

(039) 


where 


N (t)  • n (t)  cos  9 (t)  + n (t)  sin  9(t) 
c c • s 

N (t)  - -n  (t)  sin  8(t)  + n (t)  cos  0(t) 
sc  a 


(040) 


and  Nc(t) , Hg(t)  are  independent  white  noise  processes  of  power  spectral 
density  N0/2,  bandlimlted  to  |f|  <82/2.  Let  <|>(t)  ■ 0 ( t)  - 0 ( t) , where 
0 <t)  may  be  obtained  by  frequency  division  of  the  output,  r#q(t),  of  the 
squaring  loop 

r8q(t)  - 2K^  sin[2wot  + 20 (t)]  . (041) 
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The  differential  equation  for  the  operation  of  the  squaring  loop  Is  then 
given  In  Ref.  C-3,  p.  59  by 

' 2$(t)  - 2§(t)  - KF2(s){Pad2(t)  sin  2>p(t)  + v2[t,2iKt)]}  . (C-42) 


The  noise  term,  v2[t,2iji(t)  1 * 1b  characterized  by  (see  Ref.  C-3,  p.  80) 

R <t)  ^ E[v2(t,2ij-)v2<t  + t,2*)]  - 4(PaR^  (t)  + R^  (t)  J , (C-43) 
2 c c 


where  in  the  derivation  of  Eq.  (C-43)  it  is  assumed  that  T^  » rn,  and 
Tn  is  the  correlation  time  of  the  bandlimited  input  noise  nwi(t),  so  that 
m(t)  * m(t  + x)  for  all  values  of  t of  interest,  and  that  ijf(t)  Is  rela- 
tively constant  over  such  values  of  x.  If  the  loop  bandwidth  BN  Is  much 
narrower  than  B2,  then  the  only  portion  of  V£(t)  of  interest  is  the  value 
at  zero  frequency,  which  is 


Sv  (f-0) 
v2 


(C-44) 


In  an  ordinary  phase- locked  loop,  the  equation  of  operation  is  (see  Ref. 
C-3,  p.  28) 

^(t)  - 8(t)  - KF2(a){^P^  sin  *(t)  + N[t,«j;(t)])  , (C-45) 

where  the  PSD  of  N[t,iji(t)]  is  N0/2  at  f ■ 0,  For  such  a loop,  the  rms 
value  of  the  output  phase  error  is  found  to  be,  using  linearization 
(i.e.,  sin  ip  «4>),  (see  Ref.  C-3,  p.  29) 


♦00 


rms 


r,d  • 


(C-46) 


Thus,  by  analogy  between  Eq.  (C-42)  and  Eq.  (C-45),  the  output  phase  error 
of  the  squaring  loop  is 


2*<t) 


rms 


Vfp  N /N  \5 

* ♦ (f) 


B2  r,d  ' 


(C-47) 


so  that 


♦ <t> 


rms 


■VMM 


Bn  rad 


(C-48) 
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As  discussed  In  Ref.  C-3,  p.  63,  the  Costas  loop  (see  Fig.  C-5) , 
is  theoretically  equivalent  to  the  squaring  loop,  so  long  as  the  shaping 
provided  by  G(s)  in  Fig.  C-5  is  equivalent  to  that  produced  by  the  low- 
pass  equivalent  of  D2(b)  In  Fig.  C-3.  Thus  Eq.  (C-48)  also  characterizes 
the  performance  of  the  Costas  loop. 

The  decision-feedback  loop  (Fig,  C-6)  offers  improved  perfor- 
mance over  both  the  squaring  and  the  Costas  loop.  It  makes  use  of  an 
ordinary  phase-locked  loop  with  Incorporated  delay,  Tj.  The  ±1  modula- 
tion d(t)  is  removed  from  the  input  carrier  by  the  output  of  a modulation 
estimator.  The  modulation  estimator  is  a matched  filter  set  up  to  detect 
a ±1  over  a bit  Interval  T<j.  This  matched  filter  is  merely  an  integrate- 
and-dump  circuit  followed  by  a hard  limiter.  The  delay  in  the  phase- 
locked  loop  is  needed  to  allow  the  modulation  estimator  to  make  its  de- 
cision on  the  current  bit  before  any  attempt  Is  made  to  remove  the  bit 
from  the  carrier  in  the  multiplier.  The  differential  equation  of  the  de- 
cision-feedback loop  is  (see  Ref.  C-3,  p.  65) 

£(t)  - 5(t)  - KF2(s)  (Vp>  - 2Pg[iKt)  ] ) sin  ijj(t) 

+ dest<t)NA[t*'Kt)])  ’ <C“49) 


where  the  delay  T^  has  been  ignored,  on  the  presumption  that 
BN  « Td”l.'  N|,[t,ij)(t) ] is  modeled  exactly  as  N[t,<|)(t)]  in  Eq,  (C-45). 
PeWO  J is  the  error  probability  of  the  modulation  estimator,  condi- 
tioned on  the  loop  phase  error.  Pgl<|<(t)]  is  given  by 


where 


and 


?Et*(t))  * |jl  - erf  cos[iKt)]|  > 


*b 


P T 

A ad 
N. 


00 

erf (a)  A / exp(-82)dB  . 

yir  J 


(C-50) 


(C-51) 


(C-52) 


In  the  frequency  region  of  interest  (near  zero  frequency)  N^(t,iji(t)]  is 
white.  Further,  if  the  spectrum  of  de8t;(t)  is  not  wider  than  that  of 
Nj,(t,4*(t)]  then  the  noise  term  dC8t(t)Nfc[t,^(t)]  in  Eq.  (C-49)  is  also 
white  near  zero  frequency,  with  PSD  given  by  the  product  of  the  PSD  of 
N|,[t,Y(t)]  (i.e.*  N0/2)  with  the  power  in  de8*(t)  (i.e.,  unity).  Thus 
this  noise  term  has  spectral  height  N0/2,  as  aid  N|,[t,#(t)]>  We  make 
the  approximation  in  Eq.  (C-49)  that  Pg[$(t)]  can  be  replaced  with  its 
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Quadrature  - Phase  Arm 


In-Phase  Arm 

Fig.  C-5  Costas  Loop 
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expected  value,  which  la  certainly  true  for  small  ^(t).  Tiien  the  rma 
phase  error  at  the  loop  output,  after  linearisation  of  the  loop,  ia  given 
by 


'j'(t) 


rma 


■m 


E erf[^S^  coa  ifi(t)]' 


2 

Bn  rad  , 


(C-53) 


where  BN  ia  the  loop  bandwidth  of  the  phase-locked  loop  component  of  the 
decision-feedback  loop.  If  t{i(t)  ia  small,  which  we  may  verify,  then  we 
may  approximate  coa  i|>  * 1,  ao  that  Eq.  (C-53)  becomes 


(C-54) 


The  answers  obtained  in  Eqa.  (C-48)  and  (C-54)  through  lineari- 
zation agree  quite  well  with  those  displayed  in  Ref.  C-3,  which  were 
computed  using  the  nonlinear  theory. 


PERFORMANCE  OF  DATA  CLOCK  SYNCHRONIZATION 

In  this  aection  we  evaluate  the  probability  that  the  determina- 
tion of  the  synchronization  of  the  data  clock  with  the  code  epochs  (aa 
described  in  the  text)  will  be  in  error.  The  input  to  this  procedure  is 
given  in  Eq.  (C-17)  aa: 

y31(t)  - yPTi  m1(t  - \±)  cos[n3(t)]RpRNte1(t))  + n^(t)  . (C-55) 

Because  of  the  low  frequency  nature  of  the  procedure  used,  Eq.  (C-55) 
may  be  replaced  by 

y3i(t)  “ m£(t  - Xi)E  coa(n3(t)]ERpRN[ei(t)J  + n^(t) . (C-56) 

Then,  using  Eqa,  (C-12)  and  (C-27),  Eq.  (C-56)  may  be  replaced  by 

y31(t)  - (0.93016)mi(t  - X£)  + nA(t)  (C-57) 

If  T seconds  of  data  are  used  (Instead  of  Tp,  aa  in  Eq.  (C-19)),  and  the 
procedure  outlined  earlier  ia  used,  then 

g£  » E(g£)  + (C-58) 

(call  1*  the  value  of  t for  proper  synchronisation).  Than 
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T / 

E(g^)  » J 0.93016  ( -22-5- 

o 


10“3  - 82.73  x 10"9 


20  x 10 


-3 


-)  <lt  , 


(C-59) 


where  the  code  epochs  themselves  have  an  rma  estimation  error  of  82.73  ns. 
If  l is  in  error  by  k out  of  20  epochs  (k  - -10,  -9,  ...  , 0,  1,  ...»  9), 
then 


■ (12ioM)  / 93016  V^if 

o 

E(gt>  ■ (-  20lk|)  a-”01**  1 • 


20  x 10“3  - 82.73  x io"9 


20  x 10 


-3 


dt, 


(C-60) 


The  noise,  n£,  is  given  by 
T 


n£  - j n4(t)hA(t)dt 


(C-61) 


where  h#(t)  is  a sequence  of  ±l’s,  determined  by  the  signs  of  the  {y^. } 
in  Eq,  (C-18) , The  variance  of  n^  is  then 


T T 


W -//  n4(t)n4 (yJh^CtJh^CyJdtdY 


(C-62) 


o o 


Because  of  the  long  integration  time  of  our  process,  the  50  Hz  bandwidth 
of  n^(t)  is  Ignored,  and  it  is  treated  as  though  It  were  white,  with  PSD 
of  N0/2.  Thus  Eq.  (C-62)  is  evaluated, 


(\)2  ’ / /( f)8(t  - Y)dt,i1'  -(f)’ 


(C-63) 


o o 
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Using  signal  space  theory,  Ref.  C-4,  pp.  77-82,  the  probability  of  the 
estimation  of  the  data  clock  synchronization  being  off  by  one  code  epoch 
in  a given  direction  is  the  same  as  the  probability  of  mistaking  two  sig- 
nals space  apart  by  (1/20) (0.930156) ( V P§i)T  in  noise  of  variance  TN0/2, 
which  is 


(i)2[(^H0.930156)  ^ T)2  •(,§-) 


1/2 


(C-64) 


where 

00 

Q(a)  A j ^-jexp(-B2/2)d8  . 
a 


(C-65) 


For  T * 150  sec,  the  quantity  in  Eq.  (C-65)  is  given  by 
Q (5 . 063)  m 3 x 10-7. 

This  quantity  is  the  probability  of  being  off  by  one  code  epoch 
in  a given  direction.  For  both  directions,  the  probability  la  twice 
this,  or  approximately  6 * 10“7.  For  larger  offsets  the  probability  is 
much  smaller.  This  then  is  approximately  the  probability  that  the  syn- 
chronization of  the  data  bits  with  the  code  epoch  times  will  be  incor- 
rect. 


Ref.  C-4.  J.  M.  Wozencraft  and  L.  H.  Jacobs,  Principles  of  Com- 
munication Engineering.  John  Wiley  & Sons,  New  York,  1967. 
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Appendix  D 

A SINGLE-BIT  PROCESSING  SCHEME  FOR 
SATRACK-GPS  CARRIER  TRACKING 


This  appendix  presents  the  results  of  the  first  phase  of  the  in- 
vestigation of  the  digital  receiver.  Three  main  topics  are  treated. 
First,  single-bit-quantized  shipboard  preprocessing  of  the  doppler  car- 
rier signal  is  discussed.  Second,  a specific  design  for  a single-bit- 
quantized  second-order  digital  phase-locked  loop  (one  component  of  the 
digital  receiver)  is  presented.  Third,  models  are  derived  for  the  vari- 
ous noises  present  in  this  tracking  loop,  both  those  noises  inherent  in 
the  input  signal  and  those  due  to  single-bit  quantization  in  the  tracking 
loop. 


SHIPBOARD  PREPROCESSING  AND  THE  DIGITAL  PHASE-LOCKED  LOOP 

In  the  SATRACK  scheme  a sum  of  carrier  signals  is  transponded 
by  the  missile  from  the  satellites  to  the  range  ship.  These  carriers  are 
suppressed  by  pseudorandom  noise  (PRN)  biphase  modulation  (±1  amplitude 
modulation,  with  a different  sequence  for  each  satellite)  and  ±1  data 
modulation.  In  the  postflight  receiver  it  is  desired  to  track  the  car- 
rier doppler  shifts  and  the  different  delays  of  the  PRN  sequences.  Tel- 
emetry and  surface  station  data  from  the  flight  is  provided  to  the  post- 
flight  receiver  to  aid  in  this  task.  The  signals  transponded  by  the 
missiles  and  received  by  the  ships  will  be  multiplied  by  in-phase  and 
quadrature  sinusoids  of  fixed  frequency,  sampled,  quantized,  and  stored 
on  tape  for  postflight  processing  at  APL.  In  the  investigation  to  be 
pursued  in  this  paper,  it  is  assumed  that  single-bit  quantization  is 
used,  and  there  follows  the  analytical  exploration  of  the  receiver  per- 
formance that  would  then  be  obtained.  Specifically,  the  performance  of 
the  digital  phase-locked  loop,  one  component  of  the  digital  receiver,  will 
be  discussed. 

The  shipboard  preprocessing  and  digital  phase-locked  loop  (DPLL) 
have  the  structure  shown  in  Figs.  D-l  and  D-2.  In  this  paper  the  data 
modulation  has  been  ignored  and  it  has  been  assumed  that  the  PRN  used 
for  ranging  has  been  completely  removed  from  the  single-bit-quantized 
doppler  signal  before  processing  in  the  DPLL. 


THE  SHIPBOARD  PREPROCESSOR 

The  carrier  signal  received  by  the  ship  is  first  heterodyned  to 
50  MHz  and  then  input  to  the  preprocessing  scheme  shown  in  Fig.  D-l.  The 
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input  signal  to  the  preprocessor  is  heterodyned  to  a center  frequency  of 
0 Hz  in  order  to  allow  sampling  at  as  slow  a rate  as  possible,  while 
still  somewhat  above  the  Nyquist  rate.  The  1 MHz  bandwidth  on  the  signal 
is  necessary  because  of  the  PRN  present  on  the  carrier  signal  (the  dop- 
pler  shift  itself  is  expected  to  be  less  than  100  kHz  on  either  side  of 
the  center  frequency).  The  50  MHz  low-pass  filter  is  used  to  reject  the 
double  angle  signal  and  noise  that  results  from  the  mixing  operation.  A 
frequency  of  50  MHz  is  chosen  for  the  shipboard  intermediate  frequency 
(IF)  so  as  to  avoid  averaging  the  output  of  the  mixer  in  the  low-pass 
filter  over  any  time  period  larger  than  that  over  which  the  low  frequency 
component  (1  MHz  bandwidth)  will  change  by  0.25  radians  (Ref.  D-l) . 


COMPONENTS  AND  NOISEFREE  OPERATION  OF  THE  DIGITAL  PHASE-LOCKED  LOOP 


The  principles  of  operation  of  the  phase-locked  loop  shown  in 
Fig.  D-2  are  next  discussed  without  considering,  for  the  moment,  any  of 
the  various  noises  present.  Note  that  the  excluaive-OR  (l.e.,  represent- 
ing multiplication  for  sequences  of  ±l's)  gates  1 through  4 are  used  to 
heterodyne  the  input  I and  Q channels,  centered  at  0 frequency,  up  tc  an 
IF  frequency  of  100  kHz.  First,  the  phase  detector,  which  is  composed 
of  gates  5 through  8 followed  by  a counter,  compares  the  phases  of  the 
heterodyned  input  square  wave  and  the  square  wave  generated  by  the  fre- 
quency division  of  the  binary  rate  multiplier  (BRM)  output.  The  latter 
square  wave,  in  a perfect  tracking  situation,  is  at  the  frequency  of  the 
heterodyned  input  waveform  and  90#  offset  from  it. 

To  illustrate  the  above  process  take  two  square  waves,  both  at 
a frequency  of  100  kHz,  which  differ  in  phase  by  72*  (i.e,,  18*  from  per- 
fect track)  and  imagine  them  to  be  put  through  any  of  the  exclusive-OR 
gates  in  Fig.  D-2.  The  resulting  signals  are  shown  in  Fig.  D-3.  The 
counter,  which  samples  50  times  in  every  cycle  of  Divider  1 (D-l)  output, 
sums  the  samples  and  dumps  its  final  answer  (which  is  then  held  as  input 
to  the  digital  loop  filter  until  the  next  dump)  and  resets  to  0 count 
after  exactly  500  cycles  of  the  output  of  Divider  1.  The  total  number 
of  samples  in  every  counter  dump  cycle  is  thus  25  000,  and  we  would'  ob- 
serve a net  count  of  5000  in  this  18*  offset  situation.  Note  that 


18* 


5000 
25  000 


x 90* 


> 


so  the  phase  detector,  operating  in  the  absence  of  noise  (this  is  not  true 
when  noise  is  added),  is  in  fact  linear,  unlike  the  phase  detector  in  the 
analog  phase-locked  loop  (Ref.  D-2)  which  is  sinusoidal  in  its  response. 

Ref.  D-l.  W.  E.  Larimore,  ''Design  and  PSrformance  of  a Second-Order 
Digital  Phase-Locked  Loop,"  Syaposium  on  Computer  Processing  in  Communi- 
cations , Polytechnic  Institute  of  Brooklyn,  1969,  pp.  343-356. 

Ref.  D-2,  A.  J.  Viterbi,  Principles  of  Coherent  Communication, 
McGraw-Hill  Book  Co.,  New  York,1966. 
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Fig.  D-3  Multiplication  in  a Phase  Detector  Gate:  (a)  Time  Axis; 
(b)  100-kHz  Sinusoid;  (c)  100-kHz  Sinusoid  90°  Offset 
from  b;  (d)  100-kHz  Sinusoid  72°  Offset  from  b (18° 
offset  from  perfect  track);  (e)  Product  of  b and  c in 
Phase  Detector  (sum  of  samples  balanced  to  0);  (f) 
Product  of  b and  d in  Phase  Detector  (sum  of  samples 
is  net  positive,  indicating  relative  position  shift  of  b 
compared  with  d) 
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For  a more  mathematical  explanation  of  this  linearity  see  the  section  in 
this  appendix,  "Digital  Phase  Detector  Linearity,"  The  linearity  of  the 
phase  detector  response  only  holds  for  relative  phase  shifts  between 
±90*.  The  overall  response  in  the  absence  of  noise  is  shown  in  Fig.  D-4. 
The  count  also  is  quantized  to  (90/25  000)*  ■ (3.6  x 10-3)*  because  only 
a whole  number  of  counts  can  be  observed. 


Fig.  D-4  Response  (In  Counts)  of  the  Phase  Detector  for  Different 

Relative  Phase  Shifts  of  the  Input  Waves  from  Perfect  Tracking 
(90°  Offset) 


It  is  found,  in  fact,  through  computer  simulation  that  in  the 
absence  of  noise  the  digital  phase  detector  will  respond  to  Inputs  of 
sign  {8in(2*fnTp  + ♦(nTp))>  and  sign  (cos(2rffnTp  4-  $rQ (nTp) ] ) ; 

(n  ■ 1,  ...  , 25  000),  Tp  • 0.2  ps,  with  a response  of 

25  000 

(1/90)  x “ *r  <nTp)1 

n«l  0 

counts  (assuming,  of  course,  |$(nTp)  - $r^(nTp)|  « 90*  for  all  n),  and 

is  most  accurate  when  the  frequency  f is  not  integrally  related  to  the 
sampling  rate,  so  that  the  samples  are  continually  shifting  in  position 
with  respect  to  the  zero  crossings  of  cos  2irfnTp.  It  is  assumed,  inci- 
dentally, that  sampling  is  always  done  at  a greater  rate  than  the  Ny- 
qulst  rate. 
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Two  input  tape  channels  are  processed  for  noise  cancellation, 
which  is  needed  because  of  the  zero  center  frequency  of  the  recorded 
signal.  Four  channels  are  processed  in  the  DPLL  in  order  to  cancel 
double-angle  terms  and  double-angle  noise.  Because  the  loop  operates 
with  a low  IF  compared  to  the  input  noise  bandwidth,  the  double-angle 
noise  terms  would  not  otherwise  be  attenuated  by  loop  smoothing.  When 
noise  is  added,  the  phase  detector  becomes  sinusoidal  in  its  response. 
The  performance  of  the  phase  detector  in  noise  is  discussed  in  the  next 
section. 


The  phase  detector  counter  could  be  designed  to  sample  at  a 
fixed  0.2  us  rate  and  to  dump  every  5 ms,  thus  again  giving  25  000  sam- 
ples per  dump  cycle,  but  the  counting  would  not  then  be  done  over  an  in- 
tegral number  of  cycles  of  Dl  output,  causing  a round-off  error  that  is 
best  avoided.  For  analytical  purposes,  however,  it.  is  henceforth  assumed 
that  sampling  and  dumping  is  done  at  the  above  fixed  rates  and  the  round- 
off effect  is  merely  ignored.  Because  of  the  presence  of  the  IF  fre- 
quency, the  effects  on  predicted  loop  performance  of  using  this  assump- 
tion are  expected  to  be  minimal. 

The  digital  voltage-controlled  oscillator  (VCO)  is  a binary  rate 
multiplier  (BRM)  followed  by  a frequency  divider.  The  frequency  of  the 
BRM  is  controlled  by  the  output  of  the  digital  loop  filter,  which  in 
turn  is  a digital  circuit  chosen  to  make  the  performance  of  the  overall 
loop  be  similar  to  that  of  a second-order  analog  phase-locked  loop.  The 
BRM  operates  at  a frequency  100  times  that  of  the  output  of  Dl,  hence 
Dl  smoothes  over  the  "rough"  output  of  the  BRM,  deleting  pulses  from 
its  uniform  pulse  train  reference  frequency  in  order  to  have  an  overall 
pulse  rate  corresponding  to  the  loop  filter  output.  The  effects  of  the 
nonuniformity  of  the  BRM  output  will  be  discussed  later. 

Before  discussing  the  digital  loop  filter,  some  concepts  of 
phase  locked  analog  loop  design  will  be  reviewed.  A continuous-time 
analog  loop  that  is  optimal  in  its  response  to  a step  change  in  fre- 
quency has  a linearized  closed-loop  frequency  response  (in  terms  of 
phase  angle  outputs  to  phase  angle  inputs)  (Ref.  D-3)  of 

H(s)  - ....... IP -ti , (D-l) 

t sz  . . , 

— 2 — + TS  4-  1 


Ref.  D-3.  R.  Jaffee  and  E.  Rechtin,  "Design  and  Performance  of 
Phase  Lock  Circuits  Capable  of  Near-Optimum  Performance  Over  a Wide 
Range  of  Input  Signal  and  NoiBe  Levels,"  IRE  Trans,  on  Information 
Theory.  Vol.  IT-1,  March  1955,  pp.  66-72. 
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i 


and  an  open-loop  frequency  response  of 


G(s) 


2(ts  + 1) 


2 2 
T 8 


> 


(D-2) 


where  (Ref.  D-4) 

T - , (D-3) 

2BN 

and  B^  Is  the  two-sided  closed  loop  noise  bandwidth  (He).  Let  be 
2 Hz  here. 

In  the  design  of  the  digital  loop  filter , the  linearized  open- 
loop  frequency  response  of  the  analog  second-order  loop  given  by  Eq. 
(D-2)  can  be  expressed  as  two  factors.  One  is  a factor  resulting  from 
the  analog  loop  filter  and  the  other  factor  results  from  the  analog  VCO. 

The  factor  resulting  from  the  analog  loop  filter  is 


2K(ts  + 1) 

TS 


(D-4) 


The  step  response  of  this  filter  is  given  in  Fig.  D-5.  A combination  of 
BRM's  and  counters  has  been  formulated  that  responds  to  an  input  step 
function  from  the  phase  detector  counter  (i.e.,  a sequence  of  held 
counter  dumps  of  equal  magnitude)  with  a response  approximately  that  in 
Fig.  D-5b  (aside  from  quantization,  which  will  bo  discussed  later). 


Fig.  D-B  Response  of  Analog  Loop  Filter  to  a Step  Function: 
(a)  Analog  Step  Function;  (b)  Step  Response  of  Filter 


Relf.  D-4.  R.Sydnor,  J.  j.  Caldwell  and  B.  E.  Rose,  "Frequency 
Stability  Requirements  for  Space  Communications  and  Tracking  Systems, " 
Proc.  IEEE.  Vol.  54,  No.  2,  February  1966,  pp.  231-236. 
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The  factor  in  Eq.  (D-2)  caused  by  the  analog  VCO  is 

1/Kts  . (D-5) 

This  response  relates  the  VCO  input  to  the  phase  of  the  VCO  output.  The 
phase  of  the  step  response  of  an  analog  VCO  (the  step  response  is  a con- 
stant frequency  starting  at  time  0)  is  given  in  Fig.  D-6. 

In  response  to  a digital  step  function  from  the  loop  filter,  the 
phase  of  the  digital  VCO  output  behaves  in  a manner  quite  similar  to  that 
of  the  analog  VCO  in  Fig.  D-6,  but  with  some  quantization  (to  be  dis- 
cussed later)  imposed. 

Thus  aside  from  quantization  and  other  noise  sources  the  digi- 
tal phase-locked  loop  may  be  modeled  as  an  analog  loop  (Fig.  D-7).  The 
digital  nature  of  the  loop  must  be  accounted  for  when  both  the  causes 
and  effects  of  the  various  noises  that  occur  are  analyzed.  A linearized 
model  (in  terms  of  phase  angle)  for  the  DPLL  is  given  in  Fig.  D-8 
(T  ■ 0.005  s).  The  phase  signal  input  to  the  loop  in  Fig.  D-7  is  assumed 
itself  to  be  of  low  enough  frequency  that  the  sampling  operation  does  not 
cause  any  aliasing,  and  merely  represents  a (1/T)  factor  in  frequency 
response.  Note  that  the  input  phase  to  the  model  in  Fig.  D-8,  <p (t ) , is 
defined  as  the  difference  between  the  input  phase  of  the  actual  loop  in 
Fig.  D-2  and  the  loop  phase  aiding  information,  $(t)  ■ 4>i(t)  - 4>ref(t), 
which  is  expected  to  be  low  frequency  in  nature.  Note  also  that  the 
output  phase  in  Fig.  D-8,  4>r(t),  is  the  difference  between  the  output 
phase  ({^(t),  of  Fig.  D-2  and  4>ref(t).  The  only  places  where  the  sig- 
nals in  Fig.  D-8  do  not  correspond  to  those  in  Fig.  D-2  are  at  the  loop 
input  and  at  the  VCO  input  and  output. 
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Fig.  0 7 Model  of  Noi»a*Free  Digital  Phase-Locked  Loop 


Fig.  D-8  Linearized  Model  of  Noioe-Free  Digital  Phaee-Locked 
Loop  (r  - 0.76) 
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Evaluation  of  the  frequency  response  of  the  loop  in  Fig.  D-8 
(t  ■ 0.75)  indicates  that  its  magnitude  is  within  1%  of  the  magnitude 
of  the  response  of  the  second-order  loop  in  Eq.  (D-l)  for  all  frequencies 
below  11  Hz,  at  which  frequency  the  magnitude  of  the  loop  response  is  an 
attenuation  of  greater  than  28  dB  from  the  magnitude  of  the  response 
at  zero  frequency.  The  phase  of  the  frequency  response  of  the  loop  in 
Fig.  D-8  is  within  5%  of  the  phase  of  the  response  of  the  second-order 
loop  in  Eq.  (D-l)  for  all  frequencies  below  2 Hz,  which  is  more  than 
twice  the  half-power  frequency  of  the  second-order  loop.  At  about  2.6 
Hz  the  phase  response  of  the  loop  in  Fig.  D-8  crosses  through  -90°  and 
continues  to  fall  off  at  1.8#/Hz  over  all  positive  frequencies,  while 
the  phase  response  of  the  second-order  loop  is  asymptotic  to  -90°  as 
the  frequency  increases.  It  is  found  in  this  manner  that  if  the  spectral 
composition  of  the  signal  $(t)  - <J>ref(t)  is  concentrated  below  2 Hz  the 
loop  in  Fig.  D-8  may  be  modeled  as  being  identical  to  the  second-order 
loop  of  Eq.  (D-l)  (in  the  noiseless  case).  If  the  spectral  composition 
has  other  significant  components  the  more  exact  model  (Fig.  D-8)  must  be 
used  instead. 


NOISE  SOURCES  IN  THE  DIGITAL  PHASE  PROCESSING  SCHEME 

In  this  section  the  major  noise  sources  present  in  the  digital 
processing  scheme  are  listed  and  briefly  described. 

1.  Input  Noise:  The  input  noise  is  defined  as  the  stochastic 

process  n(t)  which  additively  corrupts  the  carrier  signal  aboard  the 
range  ship  (see  Fig.  D-l).  This  noise  in  turn  is  caused  primarily  by  the 
noisy  reception  of  the  satellite  signal  aboard  the  missile,  which  tran- 
sponds the  signal  to  the  range  ship.  In  finding  the  effect  of  this  noise, 
the  much  smaller  errors  over  one  5 ms  measurement  interval  (Including 
phase  quantization)  are  also  accounted  for  because  the  ship  records  a 
sampled  signal  instead  of  a continuous  time  signal. 

2.  Shipboard  Oscillator  Instability:  The  phase  noise  caused 

by  jitter  of  the  oscillator  controlling  the  I and  Q channel  processing 
and  sampling  of  the  carrier  Input  signal  on  the  range  ship. 

3.  Instability  in  the  APL  Frequency  Standard;  The  noise 
caused  by  stochastic  Irregularity  in  the  cesium  standard  governing  the 
frequency  of  the  outputs  of  the  IF  mixer  and  BRM  in  Fig.  D-2.  Its  ef- 
fects will  be  shown  to  be  nonexistent  and  are  replaced  by  synchronous 
sampling  noise  (discussed  next). 
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4.  Synchronous  Sampling  Noises  The  errors  that  result  because 
the  sampling  in  the  phase  detector  counter  is  in  synchronism  with  the 
frequency  generation  by  the  BRM. 

5.  Phase  Texture  Noise  of  Signal  Recorded  on  Ship;  The  noise 
occurring  because  the  waveform  recorded  on  the  ship  denotes  Incremental 
phase  displacements  by  changing  the  values  of  samples  that  are  coarsely 
distributed  in  time  over  a number  of  5-ms  measurement  intervals. 


6.  VCO  Phase  Texture  Noise;  The  errors  generated  because  the 
output  phase  of  the  VCO  is  ’’rough"  in  nature.  This  happens  because  the 
BRM  output  is  formed  by  deleting  positive  pulses  nonuniformly  from  a 
uniform  pulse  train,  at  a high  frequency. 

7.  VCO  Frequency  Quantization  Noise:  The  errors  that  result 

because  the  BRM  output  frequency  is  restricted  to  quantized  levels. 

8.  Loop  Filter  Quantization  Noise;  The  errors  generated  be- 
cause operations  in  the  loop  filter  can  yield  only  quantized  results. 

A more  detailed  discussion  of  these  noise  sources  follows. 


DETAILED  TREATMENT  OF  THE  NOISE  SOURCES 
Input  Noise 

2 

In  SATRACK-GPS,  the  carrier  power,  A , received  at  the  ship  can 
range  from  -125  to  -147  dBm.  We  shall  choose  the  weakest  level  (-147 
dBm)  as  the  nominal  level  for  our  discussion.  The  input  noise  is  a 
stationary  random  process  with  a two-sided  PSD  of  Nq/2  * 10-17 *2  mw/Hz, 
and  la  considered  to  be  2 MHz  wide  and  centered  at  50  MHz  in  Fig.  D-l, 
after  having  been  previously  heterodyned  down,  along  with  the  carrier, 
from  the  nominal  transmission  frequency  of  2200  MHz. 

The  immediate  effect  of  input  noise  is  to  change  some  of  the 
samples  recorded  on  the  ship  from  positive  to  negative  (or  vice  versa) . 
Thus,  given  a sinusoidal  carrier  wave  as  in  Fig.  D-9a,  thaf  would  gener- 
ate the  sequence  of  1-blt  quantized  samples  shown  in  Fig.  D-9b,  the 
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Fig.  D-9  Effect  of  Input  Noil#  on  Sampled  Data:  (a)  Sinusoidal  Carriar 
Wave;  (b)  1-Bit  Quantized  Samples  of  Sinusoid;  (c)  Sinusoid 
with  Additive  Input  Noise;  (d)  1-Bit  Quantized  Samples  of 
Sinusoid  with  Noise 
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effect  of  the  additive  input  noise  might  be  to  change  the  sinusoid  to 
the  signal  shown  in  Fig.  D-9c,  which  would,  in  turn,  generate  the  se- 
quence of  samples  shown  in  Fig.  D-9d.  Note  that  because  the  input  noise 
n(t)  is  stationary  and  independent  of  the  carrier  wave,  the  likelihood 
of  a given  sample  being  changed  by  the  noise  varies  with  the  height  of 
the  noiseless  carrier  signal  at  the  precise  time  of  the  sample.  In 
effect,  Instead  of  representing  a pulse  of  a certain  polarity,  each  sam- 
ple now  becomes  a random  variable  with  possible  values  of  -fl  or  -1, 
each  of  these  having  certain  probabilities.  The  sequence  of  random 
variables  generated  over  a certain  5 ms  measurement  interval  would  have 
a joint  probability  distribution.  The  count  observed  at  the  phase  de- 
tector counter  output  in  Fig.  D-2  at  the  end  of  the  5 ms  would,  in  turn, 
also  be  a random  variable  with  mean  and  variance  depending  on  the  phase 
relation  of  the  VCO  output  to  the  IF~mixer  input  waveform  to  the  DPLL 
(Fig.  D-2).  Investigation  of  the  mean  and  variance  of  this  random  vari- 
able for  different  phase  relations  between  the  VCO  output  and  DPLL  input 
square  wave  will  be  pursued,  but  first  heuristic  discussion  will  be 
given. 


If  the  VCO  output  is  shifted  6 degrees  from  the  DPLL  input  in 
the  noiseless  case,  the  total  count  at  the  output  of  the  phase  detector 
counter  after  5 ms  (with  the  nominal  0.2  us  sampling  rate)  would  be 
100  000  x (6/90),  one-fourth  of  this  count  coming  from  each  of  the  four 
gates  (5,  6,  7,  8).  If  the  nth  sample  at  the  DPLL  input  would  be  a +1 
in  the  noiseless  case,  it  will  now  be  a random  variable  with  a mean  of 
p(n)  - [1  - p(n)],  where  p(n)  is  the  probability  that  the  sample  is 
still  a +1  and  1 - p(n)  is  the  probability  that  the  sample  is  now  a -1. 
The  mean  of  this  random  variable  is  2p(n)  - 1,  which  is  always  *1  (it 
is  equal  to  1 when  p(n)  ■ 1,  the  noiseless  case).  In  the  presence  of 
input  noise,  when  all  the  samples  occurring  during  the  5 ms  Interval  are 
accumulated  in  the  phase  detector  counter,  the  mean  of  the  sum  of  all 
samples  is  equal  to  the  sum  of  the  means  of  each  sample,  hence  the  mean 
is  scaled  down  by  some  factor  from  the  number  of  counts  that  would  have 
been  observed  without  noise.  It  is  possible  to  heurlstlcally  observe 
how  this  factor  would  change  as  the  angle  of  offset  of  the  VCO  output  and 
DPLL  input  from  perfect  lock  (90*  difference)  changes. 

In  Figs,  D-lOa,  b,  c,  and  d are  plotted,  respectively,  a sinu- 
soidal wave  with  no  additive  noise;  a 1-bit-quantized  version  of  this 
sinusoid,  to  ba  regarded  as  an  Input  to  the  DPLL;  a square  wave  signal 
tracking  the  input  square  wave  (90s  offset);  and  a square  wave  signal 
slightly  shifted  (6s)  from  perfect  tracking  of  the  input  square  wave. 

In  the  phase  comparator,  a product  would  be  formed  between  the  input 
square  wave  and  the  tracking  square  wave.  The  phase  comparator  results 
are  shown  in  Fig.  D-lOe  for  the  perfect  tracking  square  wave,  and  in  Fig. 
D-lOf  for  the  tracking  sinusoid  with  slight  offset.  Figure  D-lOg  shows 
the  original  sinusoid  from  Fig.  D-lOa,  with  emphasis  on  the  portions  of 
the  sinusoid  which  after  1-blt  quantization  and  multiplication  by  the 
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Fig.  D-10  Effect  of  Input  Noise  with  Phase  Comparator:  (a)  Noise- 
less Input  Sinusoidal  Wave;  (b)  1-Bit  Quantized  Version 
of  Sinusoid;  (c)  Squarewave  Tracking  Input  Sinusoid;  (d) 
Squarewave  Slightly  Offset  in  Phase  from  Perfect  Track 
of  Input  Sinusoid;  (e)  Output  of  Phase  Comparator  for 
Inputs  b and  c;  (f ) Output  of  Phasa  Comparator  for 
Inputs  b and  d;  (g)  Sinusoid  from  a,  but  Shaded  Where 
Results  of  e is  Positive;  (h)  Sinusoid  from  a,  but  Shaded 
Where  Result  of  f is  Positive 
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VCO  output  in  Fig.  D-lOc  in  the  phase  comparator  gave  positive  output  in 
Fig.  D-lOe.  Figure  D-lOh  again  shows  the  original  sinusoid  from  Fig.  D- 
10a,  but  now  with  emphasis  on  those  portions  of  the  sinusoid  which  after 
1-bit  quantization  and  multiplication  by  the  VCO  output  in  Fig.  D-lOd 
have  given  positive  output  in  Fig.  D-iOf.  The  new  portion  of  the  sinu- 
soid giving  positive  output  in  Fig.  D-lOh  but  not  giving  positive  output 
in  Fig.  D-lOg  is  the  portion  of  the  input  to  the  DPLL  that  detects  the 
phase  shift  of  the  VCO  output.  Note  that  the  new  portion  of  the  sinu- 
soid giving  positive  output  in  Fig.  D-lOh  but  not  giving  positive  output 
in  Fig.  D-lOg  is  near  the  peak  of  the  wave,  so  that  when  noise  is  in- 
troduced only  the  most  reliable  samples  are  used  to  indicate  the  phase 
offset  (an  efficient  technique).  As  the  phase  offset  increases,  the  re- 
liability of  the  samples  used  to  indicate  the  additional  offset  decreases. 
Thus  the  net  output  of  the  phase  detector  counter  every  5 ms  is  not  a 
linear  function  of  the  phase  offset  between  the  DPLL  input  and  the  VCO 
output,  but  exhibits  a decreasing  slope  as  the  absolute  value  of  the 
offset  Increases.  In  fact  this  output  is  found  to  be  sinusoidal  for  the 
input  carrier  level  and  noise  in  the  SATRACK-GPS  problem  (as  will  be 
shown) . 


The  mean  and  variance  of  the  phase  detector  output  count  after 
a 5 ms  measurement  interval  for  various  phase  shifts  between  DPLL  input 
and  VCO  output  will  now  be  evaluated.  The  narrowband  input  noise  n(t) 
in  Fig.  D-l  may  be  written  (Ref.  D-2) 

n(t)  -^^[^(t)  cos  w^t  + ^(t)  sin  wot]  , (D-6) 

where  the  PSD’s  of  n(t)  and  of  nj(t)  are  shown  in  Fig.  D-ll  (ni(t)  and 
n2(t)  have  the  same  PSD's,  and  are  independent).  Considering  the  lowpass 
filter  in  Fig.  D-l  to  fully  attenuate  the  double-angle  terms  resulting 
from  the  50  MHz  mixing  operations,  and  considering  the  sampling  follow- 
ing the  low-pase  filtering  to  be  ideal,  the  inphase  (1)  and  quadrature 
(Q)  channels  recorded  on  tape  have  the  signals  shown  in  Fig.  D-2. 

I channel : 

sign}A(tn)cos(4ii(tn)  - 0]  + n1(tn)cos  e + n2(tn)sin  ej  (D-7a) 
Q channel: 

sign  | A(t  ) sin  |>,(t  ) - 8]  - n,(t  )sin  0 + n,(t  )cos  ej  (D-7b) 
( n x n x n z n > 
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Fig.  D-1 1 Comparison  of  Power  Spectral  Densities:  (a)  Power 

Spectral  Density  of  n(t>;  (b)  Power  Spectral  Density  off 
ni  (t),  n2(t) 


The  angle  6 is  a random  phase  angle  introduced  by  the  mixing  operation. 
It  is  assumed  that  0 la  uniformly  distributed  in  probability  between  0 
and  2ir  radians.  The  results  of  the  IF  mixing  are  next  listed  as-the 
outputs  of  gates  (1,  2,  3,  4)  (Fig.  D-2) : 


gate  1: 


-A(tn)cos[4i(tn)  - * - 6 - «tn]  + A(tn)cos[^(tn)  - 6 + t|)  + wtJ'J 
sign  J -n^t  >cos<9  + <|/  + wtn>  - n2(tn)sin(0  + ♦ + wtn> 

^(tnJcoetwtn  + ♦ - 0)  - n2(tn)sin(wtn  + ♦ - 0) 


<D-8a) 


gate  2: 


sign  \ 


A(tn)ein[41(t||)  - 6 - ♦ - «tnl  + A(tn)sin[01(tn)  + ♦ - 0 + wtj 

+n,(t  )sin(wt  + ♦ - 0)  + n_(t  )cos(wt  + <|>  - 8) 
in  n a u u 

+n1(tn)ain(-wtn  - ♦ - 0)  + n2(tn)coB<utn  + i|»  4-  8) 


(D-8b) 
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gate  3: 


A(tn)cos[  ♦^(tn)  + - 0 + mtn)  + AU^coel^Ct^)  - ^ — 0 — wtn)j 

sign  / +n1(tn)cos(<Dtn  + \p  + 6)  + n2(tn>sin(wtn  + <j»  + $) 


•Hi.  (t  )cos(-wt  + 6 - <p)  + n_ (t  )sin(~wt  + 6 - i>) 
x n n a n n 


<D-8c 


gate  4: 


A(tn)sin[wtn  + t|<  + * A(tn)«ln[«tn  + * + 6 - ♦1<tn)] 

sign  ^ +n1(tn)sin(wtn  + <|i  - 6)  + n2(tn)cos(wtn  + ♦ - 8) 

■Hi,  (t  )sln(ut  -f  ii  + 6)  - n,(t  )cos(u»t  + t + 8) 
in  n a xx  n 

(D-8d) 

The  angle  ^ Is  another  uniformly  distributed  random  phase  angle,  Intro- 
duced by  the  IF  mixing  operation.  Finally,  these  DPLL  inputs  are  mixed 
with  the  VCO  outputs  in  gates  5,  6,  7,  and  8 to  yield  (only  the  gate  5 
answer  will  be  given  below) : 


gate  5: 


»<t  (t  ) - .r  <«.H  - ♦,(.„)  - 29  + 2,  + 2.«„] 

o o 

+A(tn)ai«[4r  (tn)  ♦ ♦1(t||)  - 201  + A(t„>Bin^ro<tn>  “ W + 2»  + 2utnl 

■Hi. (t  )alnj4  (t  ) - 20)  n,(t  )coa{#  (t  ) - 20] 
in  r0  n in  *o  n 

4n.(t  )eln[*  (t  ) + 2*  ♦ 2ut  1 - n,(t  )coaU  (t  ) + 2ip  + 2ut  ] 

In  r#  n nzn  t0  n n 

~n.(t  )aln[*  <t  ) ♦ 20  - 20  ♦ 2wt  ] - n,(t  )coe(4  (t  ) + 2*  - 26  + 2wt  1| 
X n n n i n n n | 

-n,(t  )aln[*  (t  )J  + n,(t  )coaU  (t  )J 
l n rrt  n * n rfl  n 


(D-9) 
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Each  of  the  three  other  gate  outputs  also  contains  A(tn)sin[$i(tn) 

- 4>ro(tn))  as  one  term  in  their  sign  (•)  function,  and  each  contains  the 

noise  term  -ni(tn)slnt^r£)(tn) ] + n2(tn)cos[<j»ro(tn)] . All  other  terms 

cancel  between  the  four  gate  outputs  (not  exactly,  of  course,  because  of 
the  presence  of  the  sign  (•)*  and  the  fact  that  sign  (a)  4 sign  (b)  4 1 
sign  (a  + b)  in  general) . It  is  necessary  to  use  all  four  gate  outputs 
because  the  IF  frequency  is  low  compared  to  the  bandwidth  of  the  input 
noise,  causing  the  double-angle  noise  terms  to  have  frequency  compo- 
nents overlapping  zero  frequency,  that  would  not  be  removed  by  the  loop, 
even  though  the  loop  is  low-pass  in  nature.  (For  example,  the  PSD  of 
the  noise  term 

ni(tn)sin[*ro(tn)  4 2<j>  4 2utn]  - n2(tn)coa[*ro(tn)  4 2t|i  4 2wtn] 

is  plotted  in  Fig.  D-12  for  an  IF  frequency  of  100  kHz.)  The  IF  frequency 
cannot  be  raised  high  enough  to  completely  prevent  this  overlap  phenome- 
non if  the  sampling  rate  (S  MHz)  is  to  be  kept  substantially  above  the 
Nyqulst  rate  for  the  total  signal  (which  would  then  be  at  a high  fre- 
quency) } also  the  PSDcs  of  the  various  noises  are  not  really  square  in 
shape  (as  was  assumed  for  simplicity)  but  decay  over  wider  frequency  bands. 


Fig.  D-1 2 PSD  of  Double  Frequency  (200  MHz)  Noise  Term  in  the 
Output  of  Gete  5.  (This  noise  is  caused  by  the  overlap  of 
e box  2-MHz  wide  centered  at  200  kHz  and  a box  2*MHz 
wide  centered  at  -*200  kHz.) 


Using  the  distributions  of  the  random  phase  angles,  and  assuming 
that  the  phase  offset  le  constant  over  a 5 m?»  time  interval  (l.e., 
♦ro^n)  “ ♦i(tn)  * d#(l  * n * 25  000)),  the  autocorrelation  function  of 
the  outputs  of  gates  5 through  8,  and  the  cross  correlation  functions 
of  these  outputs  can  be  evaluated.  In  this  computation  it  is  assumed 
that  the  terms  interior  to  the  sign  (*)  functions  are  Gaussian  (the 
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noise  terms  are,  in  fact,  Gaussian,  and  do  dominate  the  signal  terms). 
Using  the  fact  that  for  two  jointly  Gaussian  random  variables  x and  y 
(Ref.  D-5) : 


E 


(D-10) 


This  information  is  then  used  to  evaluate  the  mean-squared  output  of  the 
phase  detector  counter  (which  is  normalised  to  degrees  through  multipli- 
cation by  90/100  000  because  there  are  100  000  possible  net  counts  in 
the  5 ms  Interval  (with  one  count  occurring  every  0.2  ps  in  each  of  the 
four  channels) , and  a total  count  of  100  000  denotes  a 90°  phase  shift 
(from  a perfect  tracking  situation)  in  the  noisefree  case,  as  previously 
discussed) . 

Let  y(6)  denote  the  output  (in  degrees)  of  the  phase  detector, 
with  input  phase  offset  6*.  With  the  aid  of  a computer,  it  is  found 
that  to  four  digit  accuracy  as  4 is  varied  from  perfect  track, 


E[y(6)]2  - 0.16829  + [(0.45199)  sin  6)2  . (D~ll) 

Then,  y(6,i),  the  output  of  the  phase  detector  for  two  inputs  offset 
from  perfect  track  by  6"  over  the  1th  5 nis  measurement  interval,  can  be 
modeled  as 


y(6,i)  « np(i)  + m(«)  , (D-12) 

where  the  variance  of  the  effective  noise  at  the  phase  detector  output 
is  constant  as  6 is  varied,  and  the  mean  response  of  the  phase  detector 
to  a {*  phase  shift  is 

m(6)  • 0.45199  sin  « . (D-13) 

By  Eq.  (D-13),  near  sero  phase  offset  the  factor  attenuating  the  mean 
response  is  given  by 

0. 45199 (ir/180)  - 7.8887  * 10“3  (D-14) 

(Note  that  (sin  6*)/6*  t/180  for  8 + 0) 

and  in  order  to  have  the  remainder  of  the  loop  respond  properly  to 
a 8*  phase  shift,  a multiplication  by  the  inverse  of  the  quantity 


Ref.  D-5.  A.  Popoulis,  Probability,  Random  Variables,  and  Sto- 
chastic Processes,  McGraw-Hill  Book  Co.,  New  York,  1965,  pp.  483-484. 
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in  Eq.  (D-14) , the  factor  126.76  must  be  introduced.  This  factor  is  de- 
pendent, of  course,  on  the  relative  amplitude  of  the  carrier  A(t)  and 
the  noise  n(t).  This  dependency  is  the  reason  for  the  gain  control  in 
the  loop;  i.e.,  to  adjust  the  relative  amplitude  when  the  input  signal/ 
noise  ratio  changes  (see  section  on  gain  control  design).  Thus  the  ef- 
fective noise  variance  at  the  phase  detector  output  becomes 

nd2(i)  »(l26.762)(o.l6829)-  2704.1  degrees2  , (D-15) 

which  has  an  rms  value  of  52.001°.  Henceforth  this  noise  effectively 
added  at  the  phase  detector  output  shall  be  referred  to  as  the  "effec- 
tive input  noise". 

Experimentation  with  different  averaging  Interval  lengths,  noise 
PSD  levels,  noise  bandwidths,  and  sampling  rates  (within  a few  orders 
of  magnitude  of  the  nominal  values  of  these  parameters)  Indicates  that 
the  factor  by  which  the  mean  has  been  attenuated  by  the  noise  is  propor- 
tional to  the  square  root  of  input-noise-variance/carrier-power  ratio, 
and  is  Independent  of  averaging  interval  length  or  sampling  rate.  Thus 
the  variable  factor  that  is  to  be  inserted  in  the  loop  must  be  propor- 
tional to  1/^SNR  (see  section  on  gain  control  design) . 

With  a fixed  sampling  rate,  the  variance  of  the  effective  noise 
at  the  phase  detector  output  (before  or  after  multiplication  to  correct 
the  mean)  is  Inversely  proportional  to  the  number  of  samples  in  the 
averaging  interval  (equivalently  the  number  of  input  noise  "correlation 
lengths"  in  the  averaging  interval) , provided  that  the  averaging  inter- 
val is  at  least  ten  correlation  lengths  of  the  input  noise  (0.5  ps  per 
correlation  length,  with  input  noise  being  2 MHz  wide  at  its  center  fre- 
quency). Thus  the  effective  input  noise  may  be  considered  white  every 
5 ms,  because  additive  white  noise  would  display  this  dependence  on  the 
averaging  Interval  length  (number  of  5 ms  samples  in  this  case).  After 
the  holding  operation  following  the  phase  detector  counter  dumping,  the 
autocorrelation  (R(t)]  of  this  noise  is  triangular  in  form,  having  a 
maximum  at  t ■ 0,  falling  linearly  to  0 at  t ■ ±T  (T  * averaging  inter- 
val), and  remaining  at  0 for  |t|  > T.  Thus  the  power  spectral  density 
of  this  noise,  which  then  has  a [(sin  f)/f]2  form,  has  a level  at  zero 
frequency  Independent  of  the  interval  length.  Because  of  the  low  fre- 
quency nature  of  the  loop  [relative  to  1/ (averaging  interval)],  only 
the  zero-frequency  level  is  Important. 

If  the  bandwidth  of  the  input  noise  is  decreased,  and  the  averag- 
ing interval  length  held  constant,  less  correlation  lengths  of  the  Input 
noise  are  considered,  hence  the  variance  of  the  effective  noise  at  the 
phase  detector  output  before  multiplication  by  the  scale  factor  rises  in 
proportion  to  the  factor  of  noise  bandwidth  decrease.  The  scale  factor, 
however,  also  decreases  (because  the  total  input  noise  variance  decreases) 
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in  proportion  to  the  square  root  of  the  input  noise  bandwidth  decrease. 
Thus  the  effective  variance  of  the  input  noise  after  multiplication  by 
the  scale  factor  remains  the  same. 

It  may  be  concluded  that  the  PSD  of  the  effective  input  noise  at 
the  output  of  the  hold  circuit  is  dependent  for  its  level  near  0 fre- 
quency on  the  signal-power/ lnput-noise-spectral-density  ratio  (+25  dB 
for  SATRACK-GPS)  precisely  as  in  the  analog  processing  case,  where  Ref. 
D-2  shows  that  the  output  phase  variance  of  an  analog  phase-locked  loop, 
of  noise  bandwidth  B^,  is  given  by 

E(N2)  “ | BN  radian2  . (D-16) 

The  variance  in  Eq.  (D-15)  results  after  the  delay  detector,  which  it- 
self has  a two-sided  noise  bandwidth  of  200  Hz.  When  there  is  neither 
hard  limiting  or  sampling  of  the  input,  the  variance  to  be  expected 
after  filtering  by  a 200  Hz  noise -bandwidth  filter  is  given  by  Eq.  (D-16) 
as 


2 / 10“16-9  \ 2 

E[Nj(l) ] ~ TT^I  200  * 0.6325  radian  , 

d \2  x 10 “14*7/ 

E[N2(i)J  - 2076.38  degrees2  . 


(D-17) 

(D-18) 


The  variance  found  by  application  of  Eq.  (D-16)  ia  1.15  dB  less  than  the 
variance  observed  in  Eq.  (D-15).  Thus  the  effect  of  1-bit  processing 
relative  to  input  noise  is  a deterioration  of  1.15  dB. 

Using  Eq.  (D-16)  for  comparison,  we  now  may  state  the  dependence 
of  the  VCO  output  phase  variance  on  input  noise  in  the  1-blt  case 


E[N2(i)) 


(1.3034)  radian2 


1 (D-19) 


This  relation,  plotted  in  Fig.  D-13,  applies  as  long  as  the  loop  is 
above  threshold,  and  hence  behaves  linearly.  Threshold  for  the  DPLL 
will  be  discussed  in  later  reports. 

As  tho  sampling  rate  is  varied,  keeping  the  averaging  Interval 
length  at  5 ms,  the  effective  noise  at  the  phase  detector  output  is 
asymptotically  near  51*  rms  as  the  sampling  rates  become  high  compared 
to  the  Hyqulst  frequency.  Analysis  performed  subsequent  to  tho  work  de- 
scribed in  this  appendix  shows  that  the  deterioration  decreases  as  the 
sampling  rate  increases,  asymptotically  approaching  0.77  dB.  Thus  the 
rms  phase  error  would  approach  47.6*  (see  Section  6).  This  effective 
noise  remains  relatively  constant  as  the  sampling  rate  is  decreased 
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A2/(No/2)  (decibel*) 

Fig.  D-13  Relation  Between  the  Phaae  Uncertainty  (par  hertz  of 
loop  bandwidth)  at  the  VCO  Output  Caused  by  the 
Input  and  Input  None  (carrier  power/noiae  spectral 
density)  Ratio,  A2/ (Nq/2).  (log-log  scale) 
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(52.001®  rtns  for  5 MHz  sampling  rate  until  the  Nyquist  frequency  (2  MHz) 
is  reached.  With  sampling  rates  lower  than  the  Nyquist  frequency,  con- 
secutive samples  are  almost  uncorrelated,  and  the  effective  noise  vari- 
ance begins  to  rise  in  inverse  proportion  to  the  number  of  samples  in 
the  5 ms  averaging  interval,  which  is  decreasing  as  the  sampling  rate  de- 
creases. The  exact  dependence  of  output  phase  variance  on  sampling  rate 
(for  a 2 Hz  loop)  is  shown  in  Fig,  D-14. 

If  any  of  the  outputs  of  gates  5 through  8 are  deleted  from  the 
sum  in  the  phase  detector  counter,  the  effective  noise  is  increased  be- 
cause the  double-angle  noises  overlap  zero  frequency  and  are  not  averaged 

by  the  loop.  By  deleting  the  signals  from  gates  5 and  7 the  effective 

input  noise  variance  doubles.  The  factor  attenuating  the  mean  remains 
the  same,  because  the  noise  and  signal  levels  in  the  gate  6 and  8 outputs 
are  not  changed.  The  number  of  samples  processed , however,  is  smaller, 
so  that  less  averaging  takes  plate,  causing  an  increase  in  the  additive 
white  noise  variance  at  the  output  of  the  phase  detector  counter  before 
multiplication  by  the  scale  factor. 

Experimentation  with  different  loop  IF  frequencies  Indicates  that 
an  IF  frequency  higher  than  the  current  100  kHz  will  in  fact  result  in 
slightly  less  effective  input  noise.  This  would  not  be  true  in  analog 

processing,  because  the  double-angle  noise  terms  in  the  gate  5 through  8 

outputs  cancel  each  other.  For  one-bit  processing,  however,  the  sign 
(*)  function  prevents  exact  cancellation,  so  that  as  the  frequency  of 
the  double-angle  noise  is  raised  the  more  they  are  attenuated  both  by  the 
averaging  in  the  phase  detector  counter  and  by  the  loop  filter.  At  any 
rate,  the  attenuation  in  effective  Input  noise  caused  by  changing  the  IF 
frequency  to  1 MHz  Is  only  0.071  dB  at  the  phase  detector  counter  output. 

In  summary,  the  effect  of  input  noise  is  to  add  a noise  source  at 
the  output  of  the  phase  detector  and  to  attenuate  the  mean  output  of  the 
phase  detector,  making  its  response  sinusoidal.  Since  in  the  noise-free 
case  the  response  of  the  phase  detector  is  linear,  it  is  postulated  that 
the  mean  phase  detector  response  to  phase  offset  behaves  somewhat  as  in 
Fig.  D-15,  as  the  noise  power  grows  in  relation  to  the  carrier  power. 

It  was  also  found  that  the  effective  input  noise  at  the  phase 
detector  output  is  52.001®  rms,  with  the  nominal  signal  and  noise  levels 
in  the  SATRACK-GPS  system.  The  gain  control  necessary  to  keep  the  mean 
output  of  the  phase  detector  at  the  proper  level  (i.e.,  in  degrees,  for 
example)  is  developed  in  the  section  on  gain  control  design.  In  terms 
of  phase  noise  at  the  VCO  output,  the  input  noise  causes  1.15  dB  worse 
performance  in  the  1-blt  case  than  in  the  analog  processing  case,  and 
has  variance  at  the  output  given  by 
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Fig.  D-15  Postulated  Mean  Response  of  Phaaa  Da  ltd  or  to  Phase 
Offsets,  with  y (-  Input  SNR)  at  a Paramatar 


Shipboard  Oscillator  Inatabllity 

Tha  50  MBs  oscillator  on  tha  ship  la  subjact  to  phase  jitter  af- 
fecting the  aignal  recorded  for  processing  at  APL.  As  discussed  In  Ref, 
D-4.  this  litter  noise,  as  an  additive  phase  angle,  has  a two-sided  PSD 
of  (c/|t>P)  deg*/HX'  where  the  constant,  e,  depends  on  the  quality  of 
the  shipboard  oscillator.  Thus  the  signals  generated  by  the  shipboard 
oscillator  In  Fig.  D-l  ere  in  fact 

\/2  sin  [w  t + 8 + n (t)l 

v o osc 

yfi  cos  l«Qt  + « + no#c(t)l  . <D-20) 

"Shlptlne"  tg,  a function  of  real  tine  t,  is  defined  as 

tg(t)  8 t + 8 q(t)  (D-21) 


. jJ 
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so  that  the  signals  generated  by  the  shipboard  oscillator  may  be  written 
sin  (w  t +9) 

ypTcoa  (“QtB  + 0)  . (D-22) 

Furthermore,  the  samples  are  not  taken  at  real  times  t«  « n x 0.2  ps, 
but  at  the  real  times  when  t8  ■ n x 0.2  us;  i.e.,  q”Mts  ■ n x 0.2  us). 
Thus  the  I and  Q channel  inpuBs  to  the  IF  mixer  in  Fig.  D^l  are  actually 


I channel: 


'Afa-^t,  )1  cos  ♦1[q"1(tJ  )]  - 6 - nogc[q"1(t8  ) ]| 


n 


n 


n 


sign  < + n-Iq"1^  )]  cos  jo  + u^lq^Ct  )] 

n ' n 


Q channel: 


<+  n2[q_1(tg  )]  sin  0 + noac[q~1(ts  )] 
n n 


Alq'^t  )]  ein  )]  - O - n.-Jq^Ct  )jf 

n n "n 


sign 


.-1 


-1. 


- nx[q  *(t  )]  sin  0 + nogc[q  *(t  )] 

n ’ *n 


+ n2lq”1(t,  )]  cos  je  -f  nogc[q“1(t#  )] 
n n 


(D-23) 


Note  that  in  addition  to  addltlvely  corrupting  the  phase  angle  signal, 
♦ilq-Mt,^  J,  the  phase  jitter  also  has  an  effect  on  the  input*  noise. 

It  is  assumed  for  the  present  that  the  latter  effect  is  not  significant 
because  all  noise  terms  involving  no*c[q“l(t^)]  cancel  in  the  outputs 

of  gates  5 through  8,  presuming  that  the  loop  tracks  the  low  frequency 
**oecl<r1(tsn)]  component  of  its  Input  phase  signal  (i.e.,  the  VCO  output 

signals  at  time  are 
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sign  (cos  jwq~1(t  ) + ♦ [q  X(ta  ) ] + t|»  - 9 - nQBc[q  1(ta  )]j) 
n o n n 

and 

sign  (-sin  [u>q-1(tg  ) + <f>r  tq“1Ct  )]  + ♦ - 6 - n08CU"1(ts  )]|)  , 

non  n 

where  ty^(t)  is  the  "APL  tine,"  and  t Is  real  tine).  This  assumption 
will  be  justified  in  the  discussion  below. 


Instability  in  the  APL  Frequency  Standard 

The  frequency  standard  at  APL  is  subject  to  the  sane  type  of 
instabilities  as  the  shipboard  oscillator  discussed  above.  The  reading 
of  the  input  tape,  which  has  recorded  1-bit  entries  every  0.2  ps  accord- 
ing to  shipboard  oscillator  tine,  the  oscillator  in  the  loop  IF  mixer, 
the  input  reference  frequency  to  the  BRM,  and  the  sampling  occurring  in 
the  phase  detector  are  all  synchronously  controlled  by  the  same  frequency 
standard o (For  example,  the  frequency  division  following  the  BRM  has  a 
square  wave  transition  for  every  50  "down"  transitions  of  the  BRM  output, 
and  the  sampling  in  the  phase  detector  counter  occurs  for  every  second 
positive  transition  in  the  BRM  output.)  This  fact  has  a number  of  con- 
sequences. First,  because  the  input  tape  reading  and  the  IF  mixer  oscil- 
lator are  synchronous,  the  output  sequence  of  the  IF  mixer  is  the  same 
sequence  that  would  have  been  observed  had  the  APL  frequency  standard 
been  synchronous  with  the  shipboard  oscillator.  Next,  if  the  frequency 
standard  has  a phase  jitter  this  might  cause  a VCO  output  square  wave 
transition  to  be  delayed  in  time.  If  the  sampling  of  the  VCO  output  in 
the  phase  detector  were  asynchronous  with  the  VCO  output  itself,  a sample 
value  of  this  output  in  the  phase  detector  might  be  changed.  But  the 
sampling  and  the  output  ere  synchronous,  so  the  sampling  is  delayed  along 
with  the  VCO  square  wave  output  transition,  hence  the  sequence  of  0.2  ys 
samples  of  the  VCO  output  is  exactly  the  sequence  that  would  have  been 
generated  with  a perfect  frequency  standard.  Furthermore,  when  the  out- 
put of  the  IF  mixer  is  read  into  the  phase  detector  to  be  compared  with 
the  VCO  output,  the  reading  of  nominal  0.2  ys  samples  is  synchronous  with 
the  VCO  frequency  generation,  hence  there  are  no  phase  errors  Introduced 
in  the  sequence  of  0,2  ys  samples  of  the  comparison.  In  conclusion,  the 
output  of  the  DPLL  is  the  same  sequence  that  would  be  generated  if  the 
frequency  standard  at  APL  were  synchronous  with  the  shipboard  oscillator. 
In  other  words,  once  the  I and  Q channel  inputs  have  been  determined  on 
the  ship,  the  remainder  of  the  processing  operates  as  though  the  APL  fre- 
quency standard  and  the  shipboard  oscillator  were  in  fact  synchronous. 

All  frequency  instabilities  can  thus  be  modeled  by  considering  only  the 
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differences  between  shipboard  oscillator  time  and  real  time  (neglecting 
APL  frequency  standard  time) , and  the  output  of  the  VCO  at  APL  time 
t^(t)  Is  estimating  phases  that  occurred  at  ship  time  t8  (t)  [i.e., 
real  time  q-3(ts  )],  as  claimed  in  the  case  of  shipboard  Oscillator  In- 
stability. n 

Synchronous  Sampling  Noise 

As  mentioned  in  the  previous  paragraph,  the  effects  of  APL  fre- 
quency standard  instabilities  are  not  felt  in  the  DPLL  under  considera- 
tion, because  of  the  synchronous  nature  of  the  reading  of  the  input  tape, 
the  IF  frequency  generation,  the  reference  frequency  input  to  the  BRM, 
and  the  sampling  in  the  phase  detector.  Analogously  to  the  first  law  of 
thermodynamics,  this  loss  of  a noise  source  does  not  transpire  without 
the  introduction  of  another  noise  source.  The  synchronous  sampling  of 
the  VCO  output  introduces  a noise,  which  will  be  discussed  next. 

Because  of  the  nature  of  the  frequency  division  following  the 
BRM  (which  causes  a Bquare  wave  transition  for  every  50  down  transitions 
of  the  BRM  output,  and  the  sampling  in  the  phase  detector,  which  occurs 
for  every  second  positive  transition  of  the  BRM)  the  sequence  of  0.2  ys 
sample  pulses  would  be  unchanged  if  the  continuous  time  output  of  the 
divider  were  shifted  any  amount  less  than  0.025  ys  in  either  direction 
(the  time  interval  between  transitions  of  the  reference  input  to  the  BRM). 
This  0.025  ys  shift  is  0.9*  at  the  nominal  100  kHz  output  of  the  di- 
vider (the  IF  frequency  with  no  offset).  The  probability  distribution 
of  this  error  is  modeled  as  uniform  between  -0.9*  and  0.9*  for  a given 
5 ms  measurement  interval,  and  independent  in  successive  measurement  in- 
tervals. This  probability  distribution  yields  a variance  of  0.27  degree^, 
hence  an  rms  value  of  0.52*.  This  noise  source  enters  the  loop  at  the 
frequency  divider  output. 

Phase  Texture  Noise  of  Signal  Recorded  on  Ship 

It  was  mentioned  earlier  that  the  input  noise  computation  also 
accounted  for  those  errors  over  one  5 ms  measurement  interval  because 
the  ship  records  a sampled  signal  instead  of  a continuous  time  signal. 
However,  there  is  also  a lower  frequency  noise  source  present  in  the 
ship's  output  signal,  which  will  now  be  described. 

Imagine  the  input  signal  to  the  ship  to  be  at  a constant  fre- 
quency, but  suddenly  shifted  6*  in  phase.  In  the  absence  of  input  noise 
the  output  1-bit  quantized  pulse  train  from  the  ship  will  denote  this 
phase  shift  by  changing  the  sign  of  a fractional  part  (6/180)  of  its 
samples  to  the  nearest  whole  number  of  pulses.  So,  because  there  are 
25  000  samples  in  5 ms,  a total  of  25  000  x (6/180)  samples  will  be 
changed.  Thus  one  pulse  every  5 ms  will  change  for  a phase  shift  of 
7.2  x 10-3  degrees.  If  the  pulse  shift  were  Instead 
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■ 3.6  x 10  ^ degree  , 

then  one  pulse  in  every  two  hundred  5-ms  intervals  would  be  changed 
(i.e.,  one  sample  would  change  in  every  second  of  record).  The  5 ws  in- 
terval containing  the  pulse  that  changed  would  (if  considered  alone) 
communicate  a 7.2  * 10“^  degree  phase  shift.  However,  taken  together 
with  the  other  199  intervals  in  which  no  pulses  changed,  the  phase  shift 
communicated  ia  3.6  x 10”5  degree.  Thus  there  is  a degree  of  time  non- 
uniformity  in  the  output  from  the  ship  in  its  communication  of  phase 
shift  information.  Calculation  of  this  noise  source  as  it  would  effec- 
tively appear  at  the  output  of  the  phase  detector  counter  shows  it  to  be 
more  than  two  orders  of  magnitude  smaller  than  the  input  noise.  In  fact, 
in  the  presence  of  the  input  noise  this  noise  is  attenuated  even  further, 
as  will  be  shown  next. 

The  mechanism  (Fig.  D-10)  accounting  for  the  phase  detector  out- 
put in  the  presence  of  input  noise  is  the  weighting  of  various  portions 
of  the  input  sinusoid  as  positive  or  negative  after  multiplication  in 
the  phase  detector.  When  the  input  is  slightly  shifted  positively  in 
phase  (the  sampling  times  remaining  constant)  the  shift  might  be  too 
small  to  cause  any  samples  to  be  changed  in  a 5 ms  Interval  in  the  noise- 
free  case,  but  in  the  presence  of  input  noise  the  samples  that  after  mul- 
tiplication in  the  phase  detector  would  give  positive  output  become  more 
reliable  (are  Increased  in  magnitude  on  the  input  sinusoid) , and  the  sam- 
ples that  would  give  negative  output  after  multiplication  become  less 
reliable.  Thus  the  mean  phase  detector  response  shows  phase  change  even 
though  in  the  noise-free  case  the  change  would  not  be  seen.  Further,  if 
the  shift  were  in  fact  graat  enough  to  cause  a sample  to  change  sign  in 
the  noiseless  case  this  sample  would  be  a very  unreliable  one  when  input 
noise  is  added,  coming  near  a zero  crossing  of  the  input  sinusoid.  This 
apparent  dither  phenomenon,  caused  by  the  input  noise,  attenuates  the 
phase  texture  noise  of  the  signal  recorded  on  the  ship. 

VCO  Phase  Texture  Koise 

The  output  of  the  BRM  is  obtained  by  selective  deletion  of  posi- 
tive pulses  from  a 20  MHz  reference  square  wave  (50%  duty  cycle)  supplied 
by  the  frequency  standard.  The  selection  of  pulses  to  be  deleted  is 
controlled  by  the  BRM  gate,  which  is  set  at  the  output  count  of  the  loop 
filter  to  control  the  BRM  output  frequency.  In  the  present  design  the 
gate  has  24  bits  of  precision,  which  means  that  the  output  frequency  of 
the  BRM  may  be  set  to  quantized  frequencies  with  spacing  of  1.2  Hz  be- 
tween them  (20  x 10®  H*)/2*^  « 1.2  Hz. 

To  illustrate  the  operation  of  the  BRM,  if  the  gate  were  6 bits 
in  length,  and  were  set  to  101011  ■ 43,  then  the  positive  pulses  corres- 
ponding to  the  I'u  In  the  binary  integer  would  be  seen  in  the  following 


7.2  x IQ 
200 


-3 
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sense.  First,  consider  the  binary  Integer  above  to  be  reversed  to  ob- 
tain 110101.  Next  take  the  set  of  binary  numbers  from  0 through  63. 

In  going  from  one  number  in  the  sequence  to  the  next  if  any  of  the  binary 
digits  corresponding  to  a (1)  in  110101  is  converted  from  a (0)  to  a (1) 
then  we  see  the  corresponding  positive  pulse.  So,  in  changing  from 
000000  to  000001  we  see  the  positive  pulse  for  #1,  while  in  changing  from 
111111  to  000000  we  see  no  pulse  for  #0.  A portion  of  the  input  square 
wave  reference  to  the  BRM  and  the  output  for  this  gate  setting  are  shown 
in  Fig.  D-16.  After  every  64  positive  input  pulses  the  pattern  repeats, 
unless  the  number  in  the  gate  is  changed.  Thus  the  phase  of  the  BRM  out- 
put is  nonuniform,  and  differs  from  that  of  a uniform  50-percent  duty 
cycle  square  wave  at  a frequency  of  (43/64) f0,  where  f0  is  the  reference 
frequency  of  the  BRM.  Phase  nonuniformities  of  this  type  are  analyzed 
by  Levinge  (Ref.  D-6) . 


» miuuiMJMJiMJuiniuimu1  - 


(0.05 /is  = 1 division) 


Fig.  D-16  Relation  of  the  Input  Squarewave  Reference  to  the  BRM 
and  the  Output:  (a)  Input  Squarewave  Reference  to  BRM; 
(b)  Output  of  6-Bit  BRM  with  Gate  Set  to  43 


It  is  found  that  for  a 4-bit  gate  the  highest  peak-to~peak  phase 
error  in  a pattern  period  (16  cycles  of  the  input  reference  waveform, 
for  a 4-bit  gate)  for  any  possible  gating  number  (0  to  15)  is  2.3ir  radi- 
ans. Extrapolating  a linear  relation  (between  number  of  gating  bits  and 
highest  peak-to-peak  phase  error  in  a pattern  period  for  any  possible 
gating  number)  given  by  Levinge,  we  find  that  the  maximum  peak-to-peak 
error  with  24-bit  gating  is  15.3ir  radians.  Because  of  the  frequency 
division  of  the  BRM  output  by  100,  this  peak-to-peak  error  becomes 
15.3tt/100  rad,  which  is  27.54°  at  the  output  of  the  frequency  divider. 


Ref.  D-6.  R.  W.  Levinge,  "Frequency  Control  Using  Binary  Rate 

Multipliers  for  Automatic  Testing  on  CW  Radar  Systems,"  Proc  IEE. 

I Vol.  121,  No.  10,  October  1974,  pp.  1059-1066. 
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The  actual  value  of  the  phase  error  at  any  given  time  for  the  particular 
gating  number  causing  the  maximum  peak-to-peak  value  moves  between  the 
peak  values  of  13.77°. 

From  Levinge  (Ref.  D-6) , given  the  number  of  gating  bits,  the 
peak-to-peak  phase  error  in  a pattern  period  averaged  over  the  possible 
gating  numbers  is  about  one-half  of  the  maximum,  for  small  numbers  of 
gating  bits.  We  assume  that  this  relation  holds  for  larger  numbers  of 
gating  bits,  specifically  the  24-bit  case  of  interest,  hence  an  approxi- 
mate model  is  used,  which  says  that  the  peak-to-peak  phase  error  for  a 
particular  gating  number  is  proportional  to  the  gating  number,  with  a 
high  value  of  27.54°  for  gating  number  22Zf  - 1. 

We  imagine,  to  find  the  largest  possible  effect  of  this  noise, 
that  the  output  of  the  BRM  is  given  a new  frequency  by  changing  of  the 
gating  number  only  once  every  5 ms,  and  thus  starts  at  a random  place  on 
the  phase  error  pattern  for  that  frequency  (beginning  with  zero  phase 
error  at  that  place  and  following  the  time  evolution  of  the  phase  error 
curve  for  that  gating  number  in  order  to  incrementally  generate  phase 
error).  The  5 ms  interval  is  0.00596  part  of  one  pattern  period,  which 
is  2 2^  x 0.05  us  » 0.8389  s.  From  the  phase-error  time  response  given 
by  Levinge  (scaled  in  time  axis  and  also  in  peak-to-peak  error  to  fit 
each  of  our  gating  number  peak-to-peak  values)  we  estimate  that  the 
phase  error  (i.e.,  the  difference  between  the  phase  generated  and  the 
phase  that  would  be  generated  by  a uniform  square  wave  at  the  proper 
frequency)  accumulated  in  5 ms  for  gating  #n  will  be  [1 .0719n/ (2^  - 1)]® 
in  either  the  positive  or  negative  direction.  This  error  averaged  over 
the  5 ms  interval  is  [0.53595n/ (22^  - 1)]°  in  either  the  positive  or  nega- 
tive direction,  each  with  a probability  of  0.5.  The  rma  value  of  this 
two-point  distribution  Is  also  [0.53595n/(22*  - 1)]°.  If  each  gating  num- 
ber were  equally  probable,  the  overall  rms  value  would  be  0.3094°,  which 
would  be  a noise  source  uncorrelated  every  5 ms  and  added  at  the  fre- 
quency divider  output.  The  actual  noise  is  less  because  the  gating  num- 
ber is  changed  more  often. 

VCO  Frequency  Quantization  Noise 

The  allowable  output  frequencies  of  the  VCO  are  governed  by  the 
fact  that  the  BRM  is  gated  by  24  bits.  The  input  frequency  to  the  BRM 
from  the  standard  ie  20  MHz,  thus  the  quantization  is  (20  * 10^/224 
« 1.2  Hz.  There  is  a frequency  division  by  100  following  BRM,  so  the 
output  of  the  VCO  is  quantized  to  1,2/100  - 0.012  Hz.  In  one  5 ms  aver- 
aging interval  the  average  phase  error  that  would  be  introduced  by  a 
frequency  error  of  this  size  is  0.012  Hz  x (0.005/2)  s x 360°  «*  0.0108°. 
The  probability  density  distribution  of  this  error  is  taken  as  uniform 
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between  -0.0054°  and  -K).0054°.  The  variance  caused  by  this  noise  source 
is  9.72  x 10"6  degrees2,  or  an  rms  value  of  3.118  x 10“3  degrees  at  the 
output  of  the  digital  VCO.  This  noise  is  modeled  as  uncorrelated  every 
5 ms. 

Loop  Filter  Quantization  Noiae 

The  loop  filter  quantization  noise  described  briefly  earlier  is 
not  treated  in  this  discussion.  It  is  not  expected  to  be  at  all  signifi 
cant  in  comparison  with,  for  example,  the  input  noise. 


I 

I 

I 


SUMMARY  OF  THE  NOISE  SOURCES 

The  various  noise  sources  (as  calculated  using  the  nominal 
SATRACK-GPS  parameters  given  in  the  previous  section)  and  their  places 
of  entry  into  the  DPLL  can  be  summarized  as  follows: 

1.  Input  Noise:  Effective  input  noise  is  introduced  at  the 

output  of  the  phase  detector  counter.  The  input  noise  neces- 
sitates the  introduction  of  a multiplication  by  126.76  at 
the  phase  detector  counter  output  because  of  attenuation 
that  it  causes  in  the  mean  output  of  the  phase  detector. 

The  effective  input  noise  has  an  rms  value  of  52.001°  fol- 
lowing this  multiplication,  and  is  uncorrelated  every  5 ms. 
The  effective  input  noise  with  1-bit  processing  is  1.15  dB 
worse  than  in  the  analog  case.  This  same  number  holds  for 
the  increase  in  the  1-bit  case  of  the  phase  noise  variance 
at  the  VCO  output  caused  by  input  noise. 

2.  Shipboard  Oscillator  Instability:  This  noise  source  is  in- 

troduced into  the  phase  of  the  input  signal,  and  has  a PSD 
of  c/|wP  (in  (phase  angle) 2 /Hz].  Also  it  causes  the  sam- 
ples taken  on  the  ship  to  be  recorded  at  times  offset  from 
the  Intended  times. 

3.  Instability  in  the  APL  Frequency  Standard:  In  the  digital 

phase-locked  loop  design  no  errors  are  caused  by  this  insta- 
bility. The  loop  functions  as  though  the  APL  standard  were 
synchronous  with  the  shipboard  oscillator. 

4.  Synchronous  Sampling  Noise:  This  noise  has  an  rms  value  of 

0.52*  and  is  introduced  at  the  output  of  the  frequency  di- 
vider. It  is  uncorrelated  every  5 ms. 
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5.  Phase  Texture  Noise  of  Signal  Recorded  on  Ship:  Negligible 

because  of  dither  caused  by  the  input  noise. 

6.  VCO  Phase  Texture  Noise:  The  VCO  phase  texture  noise  has  as 
rms  value  of  0. 3094 e , and  is  introduced  at  the  output  of  the 
frequency  divider.  It  is  uncorrelated  every  5 ms. 

7.  VCO  Frequency  Quantization  Noise:  This  noise  has  an  rms 

value  of  3.118  * 10“'^  degrees,  and  is  introduced  at  the  out- 
put of  the  frequency  divider.  It  is  uncorrelated  every  5 ms. 

8.  Loop  Filter  Quantization  Noise:  The  quantization  noise  from 

the  loop  filter  is  not  treated  here  because  it  is  not  ex- 
pected to  be  significant  in  comparison  to  input  noise. 


DIGITAL  PHASE  DETECTOR  LINEARITY 

In  this  section  a more  methematiral  explanation  is  given  for  the 
fact  that  the  output  of  the  digital  phase  ..*tector  is  in  fact  linear  in 
its  response  to  the  phase  shift  between  the  two  input  square  waves,  while 
in  the  analog  case  the  response  is  sinusoidal. 

In  the  analog  case,  if  two  sinusoids  <^2  sin  (tat  + 8)  and 
•y2  cos  (cot  + are  multiplied,  we  obtain 

y(t)  * 2 sln(u>t  + 8)  cos(ut  + $)  ■ sin  (8  - $)  + sin(2wt  + 0 + $) 

(D-24) 

This  output  is  plotted  in  Pig.  D-17,  If  this  output  is  applied  to  a 
low-pass  filter  the  result  is, 

yLP(t)  “ Bin  <8  - ♦)  , (D-25) 

which  is  sinusoidally  related  to  the  phase  difference  between  the  two 
input  sinusoids.  However,  if  the  input  sinusoids  are  1-bit  quantized  be- 
fore multiplication  the  result  is 

y^(t)  • sign[^T  sin(uit  + 0)1  x sign[  ^2  cos(uit  + <j>)] 

■ sign[2  sln(u>t  + 0)  cos(wt  + $)] 

y^Ct)  - sign[sin(8  - $)  + sln(2wt  + 0 + <p) ] , (D-26) 

which  is  a square  wave.  This  output  is  also  shown  in  Fig.  D-17. 
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Fig.  0-17  Output  of  Multiplier  for  Sinusoidal  Inputs  [y (t)  ] , and  for  1-Bit 
Quantized  Sinusoidal  Inputs  [yq(t)] . 


The  probability  density  function  for  the  amount  of  time  that  a 
sinusoidal  wave  spends  at  various  levels  is  given  by 


0 ; otherwlsej 


(D-27) 


Equation  D-27  may  be  deduced  from  the  following  argument.  Given  a func- 
tion f(t)  that  is  strictly  monotonic  over  a segment  [c,d],  and  for  which 
f(c)  « a,  f(d)  ■ b,  then  the  amount  of  time  the  function  spends  between 
a and  b ia  given  by  (d  - c),  which  may  be  written  as  f"^(b)  - f“^(a). 

If,  over  the  segment  [c,d],  f(t)  has  a continuously  differentiable  in- 
verse function  f~^(*)»  then  we  know  that 


b 

f f_1(y)dy  " f_1(b)  - f_1(a)  » d - c 
a 


(D-28) 


So  if  [e,g]  c [a,b]  (or  [e,g]  c [b,a]  if  f(t)  is  decreasing),  the  prob- 
ability measure  for  the  amount  of  time  spent  between  e and  g is  given 
by: 
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if  f(’)  is  increasing  \ 


if  f(')  is  decreasing 


. 


which  may  be  abbreviated 

g 

P[e,g]  * TV^Idy  , 

e 


(D-29) 


(D-30) 


where  the  integration  is  always  in  the  positive  direction*  Thus  the. 
probability  density  function  for  the  amount  of  time  that  f(*)  spends  at 
various  levels  is  proportional  to 

where  the  constant  of  proportionality  is  determined  by  the  requirement 
that  the  total,  probability  over  all  levels  be  unity.  This  argument  is 
easily  applied  to  the  case  f(t)  ■ sln(t),  where  the  function  ia  strictly 
monotonic  over  wide  regions. 


By  Eq.  (D-27),  the  portion  of  time  during  y„(t)  in  Fig.  D-17  is 
below  zero  is  given  by 


-sin(8  - $) 


— {arcsin{-sin(8  - 4>)  I * arcsin(-l)  } 


- i (♦  - 9 + i] 

4 - e . 1 

- V +2  • 


(D-31) 


The  portion  of  time  during  which  yq(t)  is  above  zero  is  thus  given  by 


(D-32) 
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Finally,  when  averaging  (low-pass  filtering) 
of  Yq(t),  the  average  value  is 


is  done  over  many  cycles 


(D-33) 


So  that  the  output  of  the  phase  detector  counter  is  a linear  function 
of  the  phase  angle  difference  of  the  two  input  square  waves. 


DESIGN  FOR  THE  GAIN  CONTROL  FUNCTION 

Here  a design  is  described  for  the  gain  control  function  that 
will  maintain  the  scale  factor  following  the  phase  detector  in  the  digi- 
tal phase-locked  loop.  The  scale  factor  is  at  a level  that  ensures 
that  the  mean  output  of  the  phase  detector  (after  multiplication  by 
90/100  000)  is  the  phase  offset  (in  degrees)  of  the  VC0  output  from  the 
DPLL  input. 

The  output,  y(6,i)  of  the  phase  detector  (discussed  earlier  in 
this  appendix),  given  two  inputs  offset  by  6°  from  perfect  track  during 
the  ith  measurement  interval,  is  given  by  (after  multiplication  by 
90/100  000): 

y(6,i)  - m(6)  + np(i)  , (D-34) 

where 

m(6)  - 0.45199  sin  6 (D-35) 

and  tip(i)  is  a random  variable  with  variance  equal  to  (0.16829)  for  the 
nominal  S/N  ratio,  sampling  rate,  and  phase  detector  counter  averaging 
interval  of  5 ms  in  the  SATRACK-GPS  system.  The  changes  in  the  S/N 
ratio  must  be  detected  so  compensation  can  be  made  for  the  resultant 
changes  in  phase  detector  mean  output  attenuation.  If  the  VC0  output 
square  waves  are  mixed  with  the  DPLL  input  waves  with  0°  nominal  offset 
phase  Instead  of  90°  offset  phase  as  in  the  four  gates  (5  through  8)  of 
Fig.  D-2,  the  count  that  should  be  received  after  multiplication  by 
90/100  000  is 

yQ(6,i)  « 0.45199  cos  6 + n (i)  , 

po 


(D-36) 
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where,  for  the  ranges  of  S/N  of  interest  the  coefficient  in  front  of 
cos  6 is  proportional  to  V S/N,  and  the  variance  of  n (i),  (0.16829),  is 
independent  of  S/N.  0 


If  it  is  assumed  that  the  loop  is  close  to  a perfect  tracking 
condition,  then  6 w 0,  and 


y (i)  w 0.45199 
o 


JW”)  actual 
1 <S/N>nominal 


+ n (i)  , 


(D-37) 


where  (S/N)  Q ^ ^ nominal  SATRACK-GPS  signal/noise  ratio  (-41  dB 

in  a 2 MHz  bandwidth),  and  (S/N)actua^  the  actual  signal/noise  ratio. 


If  multiplication  by  90/100  000  is  now  omitted,  a net  number  of 
counts  should  be  observed  every  5 ms,  given  by 


^count^ 


502.21 


(S/N) 

(S/N) 


actual 

nominal 


+ n (i)  , 

*count 


where  the  variance  of  npcount(i)  is  207  765.4  counts2.  If  the  S/N  ratio 

is  assumed  to  change  slowly  enough  to  be  considered  constant  over,  say, 
l-aecond  intervals,  then  ycount(i)  «*y  he  averaged  over  two  hundred  5~ms 
intervals  to  get,  for  the  jth  second: 


y ( j ) 

J count, avv  J 


502.21 


J (S/l,).CtU.i 


+ n 


count ,av. 


(J)  , 


where,  because  npcQunt^)  i®  white  every  5 ms,  the  variance  of 

n . (j)  is  207  765.4/200  - 1038.827  counts2,  an  rma  of  32.23 

count, av.  J * 

counts,  not  enough  to  significantly  disturb  the  measurement  of  the  mean 

value  of  502.21  counts  that  would  be  observed  if  the  S/N  ratio  were  at 

the  nominal  SATRACK-GPS  level. 


As  was  mentioned  earlier  in  this  appendix,  we  wish  to  introduce 
a multiplier  into  the  DPLL,  having  a nominal  value  of  126.76  that  is 
Inversely  proportional  to 


J(SW.ctu.J 
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Thus  we  use 

6.366  x 10A 

y ( 1 ) 

3 count, av.  J 


as  the  multiplier,  which  would  then  be  changed  once  per  second,  to  pro- 
vide the  loop  gain  control.  Any  externally  supplied  gain  information 
that  may  be  available,  of  course,  can  be  used  to  aid  in  the  determina- 
tion of  the  loop  gain  control  factor. 


A NOTE  ON  EXPERIMENTAL  VERIFICATION 

Relative  to  the  phase  noise  present  at  the  VCO  output,  caused  by 
input  noise,  the  performance  predicted  for  the  processing  scheme  dis- 
cussed here  (in  terms  of  the  sfgnal-power/noise-density  ratio)  is  also 
predicted  for  a scheme  in  which  the  input  carrier  is  not  beat  to  0 fre- 
quency, but  directly  sampled  and  1-bit  quantized.  For  this  latter 
scheme,  processing  is  carried  out  using  only  one  channel  of  Information 
provided  the  sampling  rate  is  at  least  as  high  as  the  Nyquist  rate. 

Such  implementation  would  require  sampling  that  is  too  fast  for  our  prob- 
lem, but  that  can  be  used  conveniently  with  lower  frequencies  to  experi- 
mentally check  the  theoretical  results. 
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Appendix  E 

PERFORMANCE  ANALYSIS  OF  THE  COHERENT  BINARY  DELAY-LOCKED 
LOOP  FOR  THE  SATRACK  DIGITAL  RECEIVER 


This  Appendix  continues  the  Investigation  of  the  digital  receiver, 
begun  in  Appendix  D.  The  coherent  binary  delay-locked  loop  is  analyzed. 

A receiver  for  the  SATRACK  suppressed  carrier  signals  was  discussed  in 
Appendix  C.  It  was  an  analog  receiver  that  would  perform  the  combined 
code  and  carrier  tracking.  The  received  signal  x(t)  is  in  fact  hetero- 
dyned by  2200  MHz  to  a nominal  zero  center  frequency  (bandwidth  ±1  MHz) 
by  the  shipboard  preprocessor,  hard  limited,  sampled  at  5 MHz  (2.5  times 
the  Nyquist  frequency) , and  recorded  for  lAter  processing  by  the  post- 
flight receiver  (Appendix  D) . The  postflight  receiver,  then,  is  a digi- 
tal receiver,  with  somewhat  different  performance  than  the  analog  re- 
ceiver that  has  been  discussed,  although  the  digital  receiver  is  de- 
signed to  perform  in  a manner  similar  to  the  analog  receiver.  Our  con- 
cern in  this  appendix  is  the  second  stage  of  the  digital  receiver  (Ap- 
pendix C,  Fig.  C-4),  the  phase-locked  loop/delay-locked  loop.  Specifi- 
cally, having  already  discussed  the  digital  phase-locked  loop  (Appen- 
dix D)  the  performance  of  the  coherent  binary  delay-locked  loop  follow- 
ing hard  limiting  of  the  input  signal  is  now  discussed. 


THE  BINARY  DELAY-LOCKED  LOOP  PRECEDED  BY  SHIPBOARD  HARD  LIMITING  AND 
HETERODYNING 

The  binary  delay-locked  loop  and  associated  SATRACK  receiver  pro- 
cessing is  illustrated  in  Fig.  E-l.  As  in  Appendix  C,  in  discussing  the 
receiver  for  the  signal  from  the  ith  satellite,  the  effects  of  interfer- 
ence from  signals  transmitted  by  the  other  N-l  in-view  satellites  is 
ignored  for  the  present.  The  signal  received  by  the  ship  is  thus  assumed 
to  be  given  by 

*i(t)  ■V2PsiPRNi<t  ' V rai(t  " V Sin  lV  + *i(t)1  + n(t)  ’ (E'1) 

where  A*  denotes  the  total  time  delay  from  transmission  of  the  signal 
from  satellite  (f)  to  its  final  reception  by  the  ship.  The  function 
PRN^(t)  is  a pseudorandom  sequence  of  il's  with  a chip  length  of 
A * (10“6/1.023)  second  and  an  epoch  rate  of  1 kHz  (i.e.,  1023  chips 
per  epoch).  The  function  m^(t)  is  also  a sequence  of  ±l's  that  carries 
data  from  the  ith  satellite.  The  chip  length  of  m^(t)  is  Ajj  ■ 20  ms. 

The  phase  function  i^(t)  is  due  to  the  doppler  shifts  caused  by  satel- 
lite, missile,  and  range  ship  movement.  The  noise  function  n(t)  is 
Gaussian  white  noise  that  has  been  bandlimlted  by  the  shipboard  re- 
ceiver to  ±1  MHz  about  the  missile  transmission  frequency  of  fQ  ■ 2200 
MHz.  The  PSD  of  n(t)  is  given  by  N0/2  ■ 10-17.2  mW/Hz.  The  noise  func- 
tion n(t)  may  also  be  written  as 
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n(t)  sin  u t + n„(t)  cos  « t]  , (E-2) 

V i O / O 

where  ni(t)  and  are  independent  Gaussian  white  noise  processes, 

bandlimited  to  | f J SI  MHz,  with  PSD'e  of  N.,/2.  The  nominal  carrier 
power  for  the  ith  signal,  Ps^,  is  given  as  10-14.7  although  a wide 

variation  may  be  observed  in  the  signal  power  received  from  the  differ- 
ent satellites.  The  nominal  signal/noise  density  ratio  is  thus  given 
by  2Pg^/N0  ■ 25  dB-Hz.  The  code  sequences,  PRNj^(t) , and  their  epoch 
times,  as  transmitted  by  the  satellites,  are  known  by  the  postflight 
receiver,  although  the  delays  (A^)  are  not  known  aprlori.  In  the  GPS 
program  a set  of  37  Gold  codes  have  been  chosen  as  the  PRN  codes  to  be 
used  by  the  various  satellites.  The  doppler  shifts  (♦i(t)}  are  also 
not  known  apriori;  however,  telemetry  and  ground  station  data  are 
available  to  the  receiver,  and  are  converted  into  aiding  information 
for  the  various  receiver  tracking  loops. 

In  this  appendix  only  the  received  signal  emanating  from  one  sat- 
ellite is  dealt  with,  therefore  the  subscript  (1)  in  Eq.  (E-l)  will  be 
dropped.  Also,  the  small  effect  on  the  delay-locked  loop  of  errors  in 
the  estimation  of  mi(t  - Aj)  in  the  first  stage  of  the  receiver  will  b» 
ignored  here.  Thus  we  take  for  the  shipboard  received  signal: 


x(t)  «-J2P]  PRN(t  - X)  sin  [oi  t + «t> (t)  ] + n(t)  (E-3) 

The  shipboard  processing  is  of  the  form  given  in  Appendix  D,  Fig. 
D-l,  The  in-phase  and  quadrature  (I  and  Q)  channel  signals  received  on 
tape  from  the  ship  are  heterodyned  up  to  an  IF  frequency  of  100  kHz  be- 
fore processing  (gates  1,  2,  3,  4).  The  exclusive-OR  gates  represent 
straightforward  multiplication  for  sequences  of  ±1’».  The  principle  of 
operation  of  the  delay-locked  loop  is  discussed  in  Ref.  E-l.  Th®  delay 
detector  counter  dumps  its  output  and  resets  at  a nominal  rate  of  1 kHz, 
after  adding  together  40  000  pulses  from  gates  9 through  16  (the  pulses 
are  taken  at  0.2  us  intervals).  When  the  loop  input  is  not  hard  limited, 
Gill  (Ref.  E-2)  shows  that  the  mean  delay  detector  output  is  linear  in 
the  clock  delay  error  (e) , e(t)  A A(t)  - A(t),  for  |ej  * A/2,  with  a half 
maximum  reading  for  e * A/2.  That  is 

»(e)  - 2^r(f)  ; |t | *f  (E-4) 


Ref.  E-l.  J.  J.  Spiiker,  Jr,,  "Delay  Lock  Tracking  of  Binary  Sig- 
nals,1' IEEE  Trans,  on  Space  Electronics  and  Telemetry,  Vol.  SET-9,  March 
1963,  pp.  1-8. 

Ref.  E-2.  W.  J.  Gill,  "A  Comparison  of  Binary  Delay-Lock  Tracking 
Loop  Implementations,  IEEE  Transactions  on  Aerospace  and  Electronic  Sys- 
tems, Vol.  AES-2,  No.  4,  July  1966,  pp.  415-424. 


- 110  - 


lMt  JOHNS  MOPHlNS  UNIVERSITY 

APPUEO  PHYSICS  LABORATORY 

LAuftfi. Maryland 


Sign{  A(rn)  PRN  (tn  - M tin  <tn)  --  0)-ni<tn)  sin  0 + nj  (tn)  co*  ()} 


.4 


d 


{t_^|  Gain  Control 
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With  hard  limiting  preceding  the  loop,  a half  maximum  reading  would  be 
20  000  counts  at  the  delay  detector  counter  output  after  1 ms.  The 
digital  loop  filter  is  designed  so  that  the  closed  loop  response  is 
given  approximately  by 


H(s) 


TS  + 1 


2 2 
r s 


+ TB  + 1 


(E-5) 


This  response  is  optimal  in  terms  of  the  total  mean-squared  delay  error 
for  ramp  inputs  of  delay  in  the  presence  of  noise.  The  two-sided  noise 
bandwidth  of  the  loop,  is  then  given  (Ref.  E-3)  by 


The  necessity  for  the  abundance  of  signals  in  our  loop  is  (a)  the 
input  signal  has  been  heterodyned  to  a low  frequency  on  the  ship,  caus- 
ing input  noise  terms  to  overlap  in  the  power  spectrum,  and  (b)  the  IF 
frequency  (100  kHz)  used  in  the  postprocessing  is  not  high  enough  to  en- 
sure that  double-angle  noise  terms  found  at  the  outputs  of  gates  9 through 
16  do  not  overlap  zero  frequency.  For  these  reasons,  extra  signal  terms 
are  carried  to  allow  some  of  the  noise  terms  to  cancel  each  other. 


Gill  (Ref.  E-2)  shows  that  the  normalized  output  delay  error  of 
the  loop,  when  the  input  is  not  hard  limited,  is  given  by 


where  only  the  errors  caused  by  input  noise  have  been  aseessed.  Atten- 
tion here  will  be  restricted  to  the  errors  caused  by  input  in  the  antici- 
pation that  the  loop  bandwidth  used  will  be  sufficiently  wide  to  atten- 
uate those  noises  generated  within  the  loop.  In  this  case,  the  corre- 
sponding result  to  Eq.  (E-7)  is  to  be  determined  where  hard  limiting 
has  been  done  on  the  ship. 

From  Fig.  E-l  it  is  easy  to  calculate  the  outputs  of  gates  9 
through  16.  For  example,  the  output  of  gate  9 is  given  by 


Ref.  E-3.  R.  Syndor,  J.  J.  Caldwell,  and  B.  E.  Rose,  "Frequency 
Stability  Requirements  for  Space  Communications  and  Tracking  Systems," 
Proc.  IEEE,  Vol.  54,  No.  2,  February  1966,  pp.  231-236. 
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*,<t)  * 


(P!LN(Tn  - A + |>]  | (VvV)  + co*  1*<tn>  + i(tn>  * 2t  + 2*  + ^n1 

- co»(9<t  ) - i(t  ) - 2*  - 2«  t 1 - co*(9(t  ) + i(t  ) - 26)  PKN  (t  -») 

nn  n n n 

- n,  (t  ) co* 1 2u.t  2*  + i<t  ))  - n,(t  > »ln(2ut  ♦ 2*  + $(t  )) 

in  n n i n n n 

- 0,  (t  ) co* ( 6(t  ) - 26)  + n,(t  ) *in(i(t  ) - 26) 

In  n 2 n n 

+ n.  (t  ) co.(2nt  ♦ 2*  ♦ *(t  ) - 26]  - n,<t  ) *ln[2ut  + 2*  + i(t  > - 29) 

In  n n an  r\  n 

♦ nj(tn)  eo*(i(tn))  + 


(E -8) 


/ 


A * 

It  is  assumed  that  the  delay  error  e(t)  - X(t)  - X(t)  and  the  phase- 
locked  loop  error  $e(t)  * $(t)  - $(t)  are  constant  over  a 1 ms  counting 
interval.  The  output  of  the  delay  detector  counter  at  the  end  of  this 
interval  (after  being  divided  by  40  000  for  normalization)  is  given  by 


y 


1 

40  000 


16  5000 


£ £ 

J-9  i-1 


z. (t  + i x t ) ; 
J P 


(Tp  - 0.2  ub) 


(E-9) 


(see  Fig.  E-l) . We  evaluate  the  mean-squared  value  of  y: 

16  5000 


E(yz) 


IM*  ® g 


+ i X Tp) 


E(y2) 


/ l \2 

" \40  000/ 


16  16  4999 

50“  E + E E 2 ( 5000  - 1)  (t  - Tp) 

k,j-9  i-1 


J-9  J J 


Tk 


J (E-10) 


where  R^,  z (t)  represents  a cross-correlation  function  defined  by 

J * 


R (t)  - E[z , (t)  z,  (t  - t)}  . 
Zj*k  J * 


(E-ll) 


In  the  evaluation  of  R (t)  the  following  facts  and  approximations  are 
used : * j k 

1.  P8(t),  e(t),  and  $#(t)  are  approximated  as  constant  over  the 
delay  (t)  of  interest. 

2.  8 and  i are  uniformly  distributed  random  phase  angles. 
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3.  The  terms  interior  to  the  sign  [»]  functions  are  approxi- 
mated as  Gaussian  (the  noise  terms  are,  in  fact,  Gaussian 
and  do  dominate  the  signal  terms) . 

4.  PRN(t)  is  approximated  as  Gaussian,  in  the  sense  that  the 
Gaussian  Moment  Factoring  Theorem  (Ref.  E-4)  is  applied  to 
products  involving  this  function  (see  example  below). 

5.  The  "(3,7)  Gold  code"*  is  used  for  PRN(t),  and  its  autocorre- 
lation function  [RpRN(3  7)(t)]**  is  used  in  the  evaluation 

of  E[y2] . 

6.  n^(t)  and  n2(t)  are  zero  mean  independent  processes,  with 
identical  correlation  functions. 

Assumption  3 allows  use  of  the  fact  that  for  two  jointly  Gaussian  random 
variables  x and  y (Ref.  E-5) ; 

E[sign[x]  * sign[y]]  ■ ^ Arcs  in  (E-12) 

x y 


As  an  example  of  the  evaluation  of  Rz.z.  (*)»  we  find  (see  last  section 
of  this  Appendix)  that 


*9*9 


(T) 


P 


♦ W'  + 2 + ,) 


x|co«2U#)  + co«(2ut)  + (y) 


+ R (t) 
nl 


4 co*2(wt)  RpujjCt)  . (E-13) 


where  Rn^(r)  is  the  autocorrelation  function  of  n^(t),  given  by 

A <T)  - ^ » 2 X 106  * f 6-  ,T>  , 

ni  2 2 X 106  TTT 


(E-14) 


*ITT  recommendations 

**Private  communication  with  G.  Jamison,  APL. 

Ref.  E-4.  J.  M.  Wozencraft  and  L.  M.  Jacobs,  Principles  of  Com- 
munication Engineering.  John  Wiley  & Sons,  New  York,  1967, 

Ref.  E-5.  A.  Popoulis,  Probability,  Random  Variables,  end  Sto- 
chastic Processes,  McGraw-Hill  Book  Co.,  New  York,  1965,  pp.  483-484. 
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and  the  Gaussian  Moment  Factoring  Theorem  has  been  used  to  make  the  ap- 
proximation that 

E[PRN(t  - X +4)  x PRN<t  **  A)  * PRN(t  - A + | - T>  x PRN(t  - t - A))  * 

+ Z5  + + RP*NC'  + f + t>  * W'  + 2-')  • <E‘15) 

Given  a fixed  value  of  e,  thia  approximation  (Eq.  (E-15) ) can  be  seen 
to  be  inaccurate  for  certain  specified  values  of  t.  Nonetheless,  aver- 
aging over  the  full  range  of  5000  values  of  t in  an  averaging  period, 
the  PRN  terms  Involved  heurlstically  should  exhibit  Gaussian  behavior, 
justifying  the  approximation.  Another  possible  approximation  in 
Eq.  (E-15)  used  in  a different  context  by  Spilker  (Ref.  E-l) , is  to  take 
analytically  a low  frequency  average  of  the  signal  components  of  the 
gate  outputs  in  anticipation  of  the  averaging  to  be  performed  by  the 
counter  and  subsequent  division  by  AO  000.  Thus,  for  example,  in 
Eq.  (E-8)  we  would  make  the  approximation 

PRN (t  - A - |)  PRN(t  - A)  * + f)  (E~16) 


This  method  leads  to  results  within  2X  of  those  obtained  using  Eq.  (E-15), 
with  the  delay  detector  counting  interval  of  1 ms. 


Using  a computer  program  to  evaluate  the  quantity  in  Eq.  (E-10) 
(actually  a scaled-down  version  of  this  quantity) , it  is  found  that  as 
e is  varied  (letting  the  phase-locked  loop  error  ♦ H 0,  for  the  present 
analysis) : 


E[y2(e)]  - 4,5225  x 10"5  + 


5.0496  x 10“3(|-) 


(E-17) 


Thus  y(e,  j),  the  output  of  the  delay  detector  for  a clock  offset  of  e 
over  the  jth  l-ms  counter  interval,  may  be  modeled  as: 

y(e,  J)  " (j)  + m(e)  , (E“18) 


where  the  variance  of  the  effective  noise  at  the  delay  detector  output 
is  constant  as  g is  varied. 


n,  (j) 

da 


2 - 4,5225  x lo"5 


(E-19) 
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and  the  mean  response  of  the  delay  detector  to  a delay  shift  c in  the 
input  signal  is  given  by 

m(e)  - (5.0496  x lO-3)  <£)  ; |e|  (E-20). 

In  order  to  reobtain  the  half-maximum  reading  at  e - 4/2,  which  occurs 
at  the  delay  detector  output  when  no  hard  limiting  is  done,  gain  con- 
trol multiplication  by  198.04  * (5.0496  x 10~3)-1  must  be  introduced. 

This  factor  is  dependent  on  the  relative  power  of  the  carrier  P8(t)  and 
the  noise  n(t),  which  is  the  reason  for  putting  the  gain  control  in  the 
loop  [the  gain  control  adjusts  this  constant  when  the  input  signal/noise 
ratio  changes  (see  Appendix  D)].  Note  that  the  scale  factor  of  198.04 
is  approximately  w/2  times  the  scale  factor  126.76  which  was  used  in 
the  digital  phase- locked  loop  (Appendix  D) . This  difference  Is  caused  by 
the  fact  that  phase  information  in  the  hard-limited  signal  is  contained 
in  the  relatively  reliable  samples  recorded  near  carrier  sinusoid  peaks 
(see  Appendix  D) , while  delay  information  is  contained  in  samples 
spread  uniformly  over  the  carrier  sinusoid.  Note  that  each  heterodyning 
operation  performed  on  the  signal  doubles  the  number  of  noise  terms 
associated  with  each  signal  [there  are  8 in  Eq.  (E-8)] , and  also  doubles 
the  number  of  channels  that  are  carried  (one  signal  received  by  the  ship, 
two  channels  recorded  on  the  ship  after  heterodyning  to  zero  frequency, 
four  channels  carried  after  the  post-processing  IF  heterodyning  opera- 
tion) . Each  of  these  operations  has  the  canceling  effects  of  raising 
the  scale  factor  by  VT,~  and  doubling  the  total  number  of  available  sam- 
ples. From  Eq.  (E-19) , the  effective  poise  variance  after  gain  control 
multiplication  is 

E[nd(i)]2  - (198. 04)2  (4.5225  * 10~5)  - 1.7737,  (E-21) 

This  variance  results  after  the  delay  detector,  which  itself  has  a two- 
sided  noise  bandwidth  of  1 kHz.  In  the  case  without  hard  limiting,  the 
variance  to  be  expected  following  filtering  by  a 1 kHz  noise  bandwidth 
filter  is  given  by  Eq.  (E-7)  as 

, / N w . n3  v $ 

E <f>  " 1-5811  • (*-22) 


Thus  the  hard  limiting  has  resulted  in  a variance  degradation  of 


1.7737 

1.5811 


1.1218 


0.4991  dB  . 


(E-23) 
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Thus  for  the  delay-locked  loop  under  consideration,  with  shipboard  hard 
limiting,  we  find  that  the  output  delay  error  caused  by  input  noise  has 
variance  given  bys 

E (f)2  «(^)(~)  (1.1218)  . (E-24) 

The  above  answer  was  obtained  using  a shipboard  sampling  rate  of 
5 MHz  (2.5  times  the  Nyquist  frequency).  A study  was  made  of  the  ef- 
fects of  changing  this  sampling  rate.  Graphs  of  the  normalized  delay 
variance  at  the  loop  output  (for  a 0.5  Hz  loop)  are  given  on  different 
scales  in  Figs.  E-2  and  E-3.  Note  that  as  the  sampling  rate  is  in- 
creased the  variance  approaches  a lower  limit  (*» 8 x 10~*0 , as  expected, 
because  unquantlzed  sampling  at  the  Nyquist  rate  would  be  sufficiently 
fast  to  reconstruct  the  entire  received  signal.  The  lower  limit  vari- 
ance is  very  close  to  the  variance  that  would  be  predicted  using  Eq. 

(E— 7 ) for  a loop  with  no  prior  hard  limiting  or  sampling.  Yet  Spilker 
(Ref.  E-l)  predicts  a 0.4  dB  degradation  in  the  case  where  band-limited 
white  noise  is  directly  added  to  a binary  PRN  sequence  that  is  then  hard 
limited  (not  sampled)  before  being  input  to  a binary  delay-locked  loop. 
The  fact  that  the  result  is  more  optimistic  is  likely  because  of  the  1 
assumption  by  Gill  (in  the  derivation  of  Eq.  (E-7))  that  the  bandwidth 
of  the  loop  input  is  substantially  wider  than  the  PRN  bandwidth,  caus- 
ing the  PRN  x input  noise  terms  to  have  higher  spectral  density  than 
in  this  case,  where  the  input  to  the  loop  is  band  limited  so  as  only  to 
include  the  main  lobe  of  the  PRN  spectrum  (±1  MHz).  Also,  for  rates 
lower  than  the  Nyquist  rate  the  variance  is  inversely  proportional  to 
the  sampling  rate,  because  the  noise  components  of  consecutive  samples 
are  then  almost  uncorrelated.  It  may  be  noted  from  Fig.  E-2  that  if 

the  shipboard  sampling  rate  were  lowered  from  5 MHz  to  2 MHz  (the 

Nyquist  rate),  the  deterioration  in  coherent  delay-locked  loop  perfor- 
mance would  be  approximately  1.33  dB. 

Incidentally,  use  of  a maximal  length  feedback  shift  register 
sequence  (Ref.  E-l)  for  the  PRN  function  instead  of  the  Gold  code  in- 
troduces less  than  0.2Z  change  in  the  answers.  The  difference  is  ex- 

pected to  be  more  significant  in  terms  of  intersignal  interference, 
which  has  not  been  discussed  here,  since  the  Gold  codes  have  been 
chosen  for  the  GPS  program  due  in  part  to  their  low  cross-correlation 
properties. 


EVALUATION  OF  THE  AUTOCORRELATION  FUNCTION  RZ9z9(t) 

R_  _ (t)  - E[z„(t)  zQ(t  - t)]  . (E-25) 

Z9Z9  9 9 


- 118  - 


I 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


g 0.001 

(0 


« (1/f) 


0.0001*- 

0.1 


Nyquist  Proposed 
Frequency  Sampling  Rate 


Sampling  Rate  (MHz) 


Fig.  E-2  Normalized  Delay  Variance  at  Output  of  0.5-Hz  Loop  as  a Function 
of  Shipboard  Sampling  Plata,  (log-log  scale) 
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Ki 


Froa  Eq.  (E-8)  of  the  text,  Eq.  (E-25)  may  be  written  as  (use  assump- 
tions 1,  2,  and  6 


R,  * <T> 

*9  9 


[PWKc  - X + |)  PRM(t  - x - X + |)]  {P#[co.2  ♦,  + | co.  2ut 

+ ~ co«  2ut  + (PRM<t  - X)  PM(t  - X - t>] 

+ * <t)  co.2(i  - 28)  + t R (t)  20.  2ut  +t»  (t)  co.  2ut 

*•  «2 

+ R (t)  .In2  (9  - 26)  + \ R (t)  co.  2wt  tjl  (t)  co.  2ut 

+ R (i)  co.*  ♦ + R (t)  .in2  ♦) 

"l  n2 


(E-26) 


Using  the  fact  that  n^t)  and  n2(t)  have  identical  correlation  functions, 
and  the  approximation  in  Eq.  (E-15),  Eq.  (E-13)  is  quickly  obtained. 
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Appendix  F 

PERFORMANCE  ANALYSIS  OF  THE  NONCOHERENT  CARRIER  AND  CODE 
TRACKING  LOOPS  FOR  THE  SATRACK  DIGITAL  RECEIVER 


INTRODUCTION 

This  appendix  continues  the  analysis  done  in  Appendixes  C,  D,  and 
E.  The  performance  of  the  first  stage  of  the  digital  receiver  is  inves- 
tigated. In  Appendix  C a receiver  for  the  SATRACK  suppressed  carrier 
signals  was  discussed.  The  discussion  Included  an  analog  receiver  that 
would  perform  the  combined  code  and  carrier  tracking.  The  received  sig- 
nal x(t)  is  in  fact  heterodyned  by  2200  MHz  to  a nominal  zero  center 
frequency  (bandwidth  11.073  MHz)  by  the  shipboard  preprocessor.  It  is 
then  hard  limited,  sampled  at  5 MHz  (2.33  times  the  Nyqulst  frequency), 
and  recorded  for  later  processing  by  the  postflight  receiver  (Appendix  D) . 
A ±1  MHz  bandwidth  was  considered  in  Appendix  D,  but  ±1.073  MHz  is  used 
henceforth,  as  this  is  the  actual  code  width,  plus  possible  doppler  shift. 
The  postflight  receiver,  then,  is  a digital  receiver,  with  somewhat  dif- 
ferent performance  than  the  analog  receiver  that  has  been  discussed,  al- 
though the  digital  receiver  is  designed  to  perform  in  a manner  similar 
to  the  analog  receiver.  In  Appendix  D and  Appendix  E the  performance  to 
be  expected  from  the  component  loops  of  the  second  stage  of  the  digital 
receiver  was  treated.  This  appendix  discusses  the  component  loops  of 
the  first  stage  of  the  digital  receiver,  the  noncoherent  delay-locked 
loop,  and  the  suppressed-carrler  tracking  loop,  and  then  evaluates  the 
overall  performance  to  be  expected  of  this  stage  of  the  receiver. 

The  Noncoherent  Delay-Locked  Loop  Preceded  by  Shipboard  Hard  Limiting 
and  Heterodyning 

A diagram  of  the  noncoherent  delay-locked  loop  for  the  SATRACK 
digital  receiver  is  given  in  Fig.  F-l.  The  shipboard  preprocessing  is 
of  the  form  given  in  Appendix  D,  Fig.  D-l.  The  excluslve-OR  gates  rep- 
resent straightforward  multiplication  for  sequences  of  tl’s.  The  input 
signal  to  the  ship  is  assumed  to  be  of  the  form  given  in  Appendix  E, 

Eq.  (E-3) . 

xt(t)  -^2\(t)  PRN1(t  - At)  m1(t  - At)  sin[wot  + ♦1(t)]  + n(t),  (F-l) 

where  interfering  signals  have  been  ignored,  but  the  data,  mj.(t  - A^), 
are  included.  Discussion  of  the  factors  and  parameters  in  this  signal 
is  given  in  Appendix  E.  Note  that  the  input  I and  Q channels  are 
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heterodyned  down  to  a nominal  0 Hz  in  the  loop  multiplication!  and  that 
the  mixing  frequency  contains  aiding  information  obtained  from  telemetry 
and  ground  station  data.  The  resulting  frequency  is  intended  to  be 

0 Hz  ±20  Hz,  where  a 20  Hz  error  in  the  siding  frequency  is  allowed  here. 
A shipboard  recording  IF  frequency  of  100  kHz  is  used,  instead  of  the 
zero  IF  frequency  mentioned  in  the  introduction.  Such  a shipboard  IF 
frequency  (u>2  “ 2tt  * 100  kHz)  simplifies  the  second-stage  processing. 

The  I and  Q recorded  channels  thus  change  from  those  given  in  Appendix  D, 
Eq.  (D-7)  to  the  following: 

1 channel : 

•ig“  | W P,!N(tn  ' Xi>  mlUn  - V co,[W  + “2*»  ' 91 

+ n^t^)  cc>u(w2tn  - 0)  - n2(tn)  sin(us2tn  - 8)| 

Q channel : ' 

sign  jAi(tn)  PRN(tn  - X£)  m1(tn  - X£)  sinl^t^  + - 9] 

+ n^t^  sin(u2tn  - 0)  + n2<tn)  cos(o>2tn  - 0)| 

The  necessity  for  the  multiple  signals  in  Fig.  F-l  is  caused  by  the  fact 
that  we  have  heterodyned  the  Input  signal  on  the  ship  to  a low  frequency, 
causing  input  noise  terms  to  overlap  in  the  power  spectrum. 

In  the  noncoherent  delay-locked  loop  shown  in  Appendix  C,  Fig. 

02,  there  is  only  one  input  signal,  and  the  bandvldths  of  FRNi(t)  and 
n(t)  are  assumed  narrow  in  comparison  with  the  IF  frequency.  Reference 
F-l  discusses  the  principal  of  operation  of  the  noncoherent  delay-locked 
loop  given  In  Appendix  C,  Fig.  02.  The  digital  loop  In  Fig.  F-l  works 
on  the  same  principles  as  this  loop,  but  suffers  some  degradation  in 
performance  because  of  the  1-bit  quant izet ion  and  sampling  prior  to  the 
shipboard  recording.  The  level  of  this  degradation  will  be  determined 
next. 


The  signal  et  the  output  of  the  center  bandpass  filter  in  Appen- 
dix C,  Fig.  02  is  given  by  (see  Rtf.  F-l) 


A 

AiRFRNjL^Xi  " Ai*  co»l*di<t)]  + »r(t)  . 


(F-2) 


Ref.  F-l.  W.  J.  Gill,  "A  Comparison  of  Binary  Delay-Lock  Tracking 
Loop  Implementations,"  IEEE  Transaction!  on  Aaroapaee  and  Electronic 
Systems , Vol.  AES-2,  No.  A,  July  1966,  pp.  415-424. 
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where  RpRN^*)  is  the  autocorrelation  function  of  PRNj^t),  and  $dl<t)  is 
the  error  in  the  phase-aiding  reference  signal  in  Appendix  C 

♦d  (t)  - MO  - ♦ f (t)  . 
i i 

The  term  mj,(t  - Xy  is  changed  to  m^Lp(t  - X^)  to  account  for  the  effect 

low-pass  filter  on  the  data.  nr(t)  is  white  noise  of  spectral  height 
N0/A,  which  is  bandliaited  to  fi  - Bj^/2  £ |f|  < fi  + Bi/2.  The  signal 
corresponding  to  this  one  in  the  digital  case  is  found  at  the  output  of 
Lowpass  (y  in  Fig.  F-l.  Analysis  is  done  similar  to  that  in  Appendix  E 
in  the  case  where  mi(t  - X*)  is  ignored  and  <o»x  £ 0,  4>d*(t)  - $Q).  This 

signal,  after  sunning  for  1/139.997  sscond  (a  total  of  71  430  possible 
counts),  is  given  by 

461.8626  Rp^  (X^  - X^)  cos  *f  n#(t)  , (F-3) 


where  n„(t)  has  a variance  of  124326.4  counts?.  A summation  Interval 
of  1/40  second  (i.e,,  a 70  Hz  low-pass  filter)  was  picked  here  because 
a low-pass  filter  was  needed  that  passes  the  data  m^(t  - X*),  whose 
spectrum  is  significant  to  50  Hs  past  the  highest  frequency  (20  Hz)  that 
the  aided  carrier  can  assume.  A short  explanation  of  the  result  in 
Eq.  (F-3)  is  given  in  the  last  section  of  this  appendix.  Because  d±(t) 

will  be  a frequency  within  the  passbsnd  of  the  low-pass  filter  following 
the  multiplication,  we  infer  that  in  the  digital  case  the  signal  corre- 
sponding to  Eq.  (F-2)  will  be 

461,8626  m (t  - X ) IL,  (X.  - X ) cos[$.  (t) ] + n (t) , <F-4a) 

iLp  i riwi  i i dt  r 


where  nr.(t)  has  a variance  of  124326.4  counts?.  This  signal  contains 
noise  terms  overlapping  0 frequency.  Thus  we  keep  a quadrature  version 
of  Eq.  (F-4a)  (output  of  Lowpass  0 in  Fig.  F-l). 


461.8626  mt  (t  - X^  (X1  - \±)  sin[*d  (t)  ] + nr  (t) 
LP  1 12 


(F~4b) 


where  nr2(t)  is  in  quadrature  with  nri(t).  The  signal/noise  density 
ratio  of  each  of  the  quadrature  signals  of  Eq.  (F-4)  is  given  by 


(1/2)  (461. 8626) 2 
(124326.4)  (1/139.997) 


120.102 


20.7955  dB»Hz  . 


(F-5) 
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This  ratio  la  1.1942  dB  worse  than  the  signal/noise  density  ratio  of 
aJ/No  " 21.9897  dB  • Hz  that  would  be  obtained  if  the  processing  were 
done  In  continuous  tine  without  quantization  (using  a 0 IF  frequency , 
as  we  have,  and  keeping  I and  Q channels).  The  signals  in  Eq.  (F-4) 
correspond  to  the  DLL  (delay-locked  loop)  estimate  of  the  PRN  epoch 
times.  These  signals  are  sent  to  the  next  loop  for  further  processing. 
The  corresponding  signals  with  early  and  late  code  multiplication  are 
processed  in  the  DLL.  Using  analysis  similar  to  that  of  Gill  (Ref.  F-l) , 
it  is  found  that  the  loop  rms  delay  error  in  our  case  is  given  by 


a 

A 


■w 


(10°*11942)  [N  (100,11942)]2  (2B 
° ^ + _£ 


rms 


2Afc<\P>' 


1/2 


(F-6) 


where  By  is  the  bandwidth  of  the  low-pass  filters  (By  - 70  Hz)  following 
the  multiplication.  The  factors  of  (lO^* 11942)  in  Eq,  (F-6)  represent 
the  1.1942  dB  loss  due  to  digital  processing.  Otherwise,  the  result  is 
identical  to  Gill's  (see  Appendix  C,  Eq.  (C-8)),  with  2By  replacing  the 
IF  filter  bandwidth  in  his  result.  The  derivation,  however,  is  changed 
because  the  I and  Q signals  are  present.  An  approximation  of  E(ra£p)  may 
be  obtained  by  approximating  miTP(t  - ) as  the  result  of  passing 

mi(t  - A*)  through  a 50  Hz  low-pass  filter.  The  worst  performance  can  be 

square  wave  with  a 40  ms  cycle  time.  In  that 
case  E(m£p)  - 0.81059.  If  the  data  are  more  random  in  nature,  the  cor- 
responding autocorrelation  function  might  be  a triangle,  extending  to 

0.9029.  The  latter  will  be  used  for  the  pres- 
former  is  a better  approximation,  then  By 
could  be  widened  to  pass  the  data  more  effectively,  hence  increasing 
E(«2p). 


expected  if  the  data  are  a 
■S)  - 0.8105! 

Let ion 

T ■ ±20  ms.  Then  E(m? „)  ■ 
ent  approximation.  It  the 


The  discussion  above  has  been  based  on  a sampling  rate  of  5 MHz. 
Lowering  the  sampling  rate  results  in  additional  degradation.  In  Fig. 

F-2  the  factor  increase  in  the  signal/noise  density  ratio  of  Eq.  (F-4) 
is  plotted  as  a function  of  the  sampling  rate.  Observe  that  the  line  in 
Fig.  F-2  that  refers  to  the  +21.9897  dB  • Hz  signal/noise  density  ratio 
obtained  by  continuous  time  processing  without  any  quantization;  also, 
note  that  with  5 MHz  sampling  the  result  falls  1.1942  dB  below  this  line, 
as  derived  above.  This  is  the  factor  100*11942  in  gq,  (F-6),  with  2.5 
MHz  sampling  (NL.16  * Nyquist  rate).  Fig.  F-2  shows  that  the  degradation 
is  1.9991  dB.  Much  storage  could  be  saved  by  sampling  at  this  lower 
rate,  Eauatlon  (F-6)  then  becomes 


MS 


rms 


fN  (100,19991)  (N  (100,19991)]2  (2B 

O , O . W 

1 2 r _.  2 v •»  T 99/—%  v'e  2 


2AI(E(mLP)] 


VE(mLP)] 


M" 


(F-7) 
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Fig.  F-2  Degradation  of  Signal/Nolta  Density  Ratio  in  Eg.  F4  at  a 
Function  of  Sampling  Rata 
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Next  the  form  that  Eq.  (F-4a)  would  assume  will  be  determined  for  a 2.5 
MHz  sampling  rate.  Recall  that  summing  is  over  1/139.997  second.  Pre- 
viously (with  5 MHz  sampling)  there  was  a maximum  possible  net  count 
(from  two  channels)  of  (5  x 10^) (1/139. 997) (2)  - 71  430  counts.  Ob- 
serve the  I and  Q answers  given  in  Eq.  (F-4) : In  the  no-noise  case  the 

mean  of  the  signal  in  Eq.  (F-4a)  would  be  given  by 
(71  430)  m^Lp(t  - X^)  RpRNi(Xi  - X ) for  4^(0  = 0.  Thus  the  scale 

factor  is  71  430/461.8626  * 154.6564.  As  mentioned  in  Appendix  D,  this 
scale  factor  is  a function  of  the  signal/noise  variance  ratio,  not  of 
sampling  rate.  Thus,  for  sampling  at  a 2.5  MHz  rate,  with  a total  pos- 
sible net  count  of  35  715  in  1/139.997  second,  the  mean  value  of  the 
count  is  given  by 

(iMJMa)  miTD(t  " Xi>  ^RN  <Xi  “ Xi)  cos 

' Lr  l i 


- (230.9313)  mt(t  - \±)  RpRN^(XJL  - XJ  cos  ^ (t)  (F-8) 

We  have  observed  that  the  signal/noise  density  ratio  would  be  degraded 
by  1.9991  dB  from  the  continuous  time  processing  case  with  no  quantiza- 
tion (the  latter  case  has  A^/N0  « 21.9897  dB  • Hz).  Thus  we  may  solve 
for  the  variance  of  the  noise,  nr„(t),  which  would  be  added  to  the  mean 
count  in  Eq.  (F-8). 

var.  (n^t))  - ‘ 74821-14  c°u"ts2  <F‘9> 

Thus,  with  2.5  MHz  sampling,  the  signal  corresponding  to  Eq.  (F-4a)  is 
given  by 

(230.9313)  mt  (t  - X±)  RpRN  (Xf  - A ) cos  ^ (t)  + nr  (t)  (F-ICa) 
LP  i 13 

and 

(230.9313)  m (t  - X ) IL,  (X,  - X.)  sin  * (t)  4-  n^  (t),  (F-lOb) 

1LP  1 11  ai  r4 

where  nr^(t)  is  characterized  by  Eq.  (F-9).  The  performance  of  the  non- 
coherent delay-locked  loop  with  2.5  MHz  sampling  was  given  in  Eq.  (F-7). 
Let  Bv  ■ 70  Hz  (50  Hz  for  the  data,  20  Hz  for  the  aiding  uncertainty), 
and  let  Bjji  ■ 1 Hz  in  Eq.  (F-7).  We  then  obtain  that 

o - (e.)  « 0.10676  A - 104.36  ns  (F-ll) 

e i rms  v ’ 
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For  convenience,  scale  the  results  in  Eq.  (F-10)  by  a factor  to 


obtain  (let  ^ Xj  “ Xj) : 

*1  "lLp(t  - V “pRNj^l5  cos  \M  + "r5(t) 


(F-12a) 


and 


Ai  mi  (t  “ V npRNj  ^ei^  8in  *d/fc*  + nr,(t)  * 

Lr  1 10 


(F-12b) 


where  nr,.(t)  and  nr^(t)  are  low  frequency  quadrature  noise  terms  of  one- 
sided bandwidth  70  Hz,  and  of  spectral  height  1.9991  dB  greater  than 
N0/2  [as  was  discussed;  i.e.,  (1.5846)  N0/2].  These  noise  terms  could 
be  written  as 


n - n (t)  cos  $ (t)  + n (t)  sin  $ (t) 

r5  a d±  b d± 

and 

n - n (t)  sin  A.  (t)  - n.  (t)  cos  Aa  (t)  (F-13) 

r,  a a,  D a. 


where  na(t)  and  nfc(t)  are  Independent  Gaussian  low  frequency  noise  terms, 
of  bandwidth  70  Hz,  and  spectral  density  (1.5846)  (N0/2).  The  next  phase 
of  the  tracking  is  a noncoherent  suppressed-carrler  loop,  to  be  dis- 
cussed next.  The  intent  is  to  track  the  phase  angle  in  Eq.  (F-12)  in 
order  to  remove  the  sinusoids  from  the  data  modulation. 


The  Suppressed-Carrier  Tracking  Loop 

The  suppressed-carrler  tracking  loop  is  shown  in  Fig.  F-3.  The 
output  signals  of  the  noncoherent  delay-locked  loop  are  first  hetero- 
dyned up  to  an  IF  frequency,  so  the  suppressed-carrler  loop  will  not 
have  to  track  through  zero  frequency.  The  loop  Itself  is  a combination 
of  a squaring  loop  and  a Costas  loop.  The  combination  is  used  because 
the  Input  noises  on  the  recorded  signal  overlap  zero  frequency.  Analysis 
of  the  combined  loop  indicates  that  its  performance  is  analogous  to  the 
performance  of  the  squaring  loop  (see  under  "Suppressed  Carrier  Track- 
ing Loops"  in  Appendix  C).  The  loop  output  is  two  signals  given  by 


y21(t)  - coslwjt  - + n3<t)1 

and 

y2i,q^  “ “ ♦di^)  + *3(0]  (F-14) 
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where  n3(t)  Is  the  phase  tracking  error.  Analysis  indicates  that  the 
rms  value  of  nj(t)  is  given  by  (for  2.5  MHz  sampling) 

f No(100,19991) 

( 2Ai(EtR1>RN  E(\p) 

(F-15) 


N (100*19991) 
0 


(2Bw) 


1/2 


N2 


The  effect  of  different  sampling  rates  on  Eq.  (F-15)  is  analogous  to 
the  effects  on  Eq.  (F— 7) . 


To  evaluate  Eq.  (F-15),  E[RpRNl(ei)]  must  be  evaluated  first.  The 

is  models 
approximately  by 


term  e-i  is  modeled  as  Gaussian,  and  uses  the  fact  that  Rprn-i  00  is  given 


RPRN1(x)  * 


(1  - JxJ-) ; 1*1 


0 ; othr-vise 

We  then  find  that  (using  Eq.  (F-ll)) 


1 - I— — 1 (0.10676)  - 0.91482 

W*T/ 


(F-16) 


(F-17) 


(F-18) 


Substituting  Eq.  (F-18)  into  Eq.  (F-15),  and  letting  loop  bandwidth 
BN2  “ 1 Hz,  we  obtain 


[n,(t) ] - 0.11308  rad  - 6.479*  . 

j rms 


(F-19) 


The  outputs  of  the  carrier  loop  are  mixed  with  the  quanti*  U-is  In  Eq. 
(F-12)  as  shown  in  Fig.  F-3,  This  mixing  attempts  to  remove  the  sinu- 
soids from  the  data.  The  result  is  given  by 


y3i<t)  - 


- *t)  cos{n3(t)J  *pgn  ♦ *»4<t)t  (F-20) 


where  n$(t)  is  a lov-nass  Gaussian  process  of  bandwidth  70  Hz  and  power 
spectral  density  io0, 19991  (Hq/2).  Appendix  C describes  the  technique 
for  estimating  the  bits  In  ni(t  - A1)  given  Eq.  (F-20),  and  evaluates 
the  probability  of  data  estimation  error,  given  Eq.  (F-20).  Because  the 
low-pass  filter  effects  mj(t  - A*),  we  must  compensate  for  the  fact  that 


i 
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E(m£p)  - 0.9029.  The  probability  of  estimating  the  bit  in  the  data 
sequence  incorrectly  is  given  by 


where 


* V*  ■ xi> 
tk 


#k 


<•  1/2^1  - E cos  n3(t)j  j . (P-21) 


\lm 


[(0.9029)  A2]  (20  x 10‘3  - 104.36  x l(f9) 

N <100’19991) 
o ' 


- 1.8019  . 

If  n-j(t)  is  modeled  as  Gaussian,  we  obtain 


(F-22) 

(F-23) 


E{cos[n3(t)] } - exp 


- (l/2)o2[n3(t)]J  - 0.99363  (F-24) 


If  Eqs.  (F-18),  (F-23),  and  (F-24)  are  substituted  into  Eq.  (F-21)  the 
value  obtained  is 


*[£4(t  - 


+ m4(t  - X4) 


fk 


] "(1/2)  [1  - erf (1.2202)]  « 0.042208  . (F-25) 


EXPLANATION  OF  EQUATION  (F-3)  RESULTS 

A short  explanation  of  the  derivation  of  the  result  given  in 
Eq.  (F-3)  is  in  order  here.  In  Appendix  E the  output  of  a coherent  de- 
lay detector  was  analyzed.  In  that  detector,  a difference  between  early 
and  late  channel  outputs  was  made.  In  the  present  situation,  only  one 
channel  la  considered;  one  channel  is  deleted  from  the  earlier  analysis. 
Also,  In  that  analysis  each  channel  is  a sum  of  four  signals  Instead  of 
two  signals  (because  of  the  extra  heterodyning  operation  that  was  made 
in  that  analysis  and  use  of  the  approximation  given  in  Appendix  E, 

Eq.  (E-16)).  If  the  four  channels  are  summed  for  1/139.997  second  (a 
total  of  142  860  possible  counts)  and  divided  by  142  860,  the  output  is 
found  to  be 

(4.811685  x 10"3)  Rpw,i(A1  - cos(^)  + nr#(t)  , (F-26). 

where  nr4(t)  has  a variance  of  1.218352  * IQ-3.  In  the  no-noise  case, 
the  mean  response  would  have  been 
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RpRN  (^  - x±)  ^1  - | - It  + [<j>d  (t>  + ir]  mod  2ir 

(the  secondAfactor  la  a cosine  with  triangular  lobes),  which  la  just 
when  ^^(t)  “ °*  So  scale  factor  is 


= - 207.8274  . 

4.811685  x 10““} 

The  signal/noise  density  ratio  Is 

. 266.036  . 24.2494  dB.Hz . 

1.218352  x 10° 


(P-27) 


(F-28) 


With  one  less  heterodyning  operation,  the  scale  factor  Is  multiplied  by 
1/VT.  This  Is  because  the  signal/noise  variance  ratio  in  each  compo- 
nent signal  is  3.01  dB  higher  (Appendix  D).  Thus,  with  a possible  total 
of  71  430  counts,  and  division  by  71  430,  the  output  is  given  by 

(6.80475  x io'3)  (At  - cos  + n^t)  , (F-29) 


where  the  variance  of  nrb(t)  may  be  determined  by  the  fact  that  the 
slgnal/nolse  density  ratio  is  unchanged  by  the  deletion  of  a heterodyn- 
ing operation. 


,-3x2 


v.r  tnrb(t))  - t1*-80^5  ■ ^.«670  . 10~5  (F-30) 


If  Eq.  (F-29)  is  translated  into  counts  through  multiplication  by  71  430, 
we  obtain 


(486.0633)  Ry^  (Xt  - k±)  cos  + ng(t) 


(F-31) 


2 

where  n#(t)  has  a variance  of  124326.4  counts  . 

A signal  degradation  of  0.4436  dB  caused  by  the  input  signal  to 
the  ship  being  bandpass  filtered  (±1.073  MHz)  must  also  be  added  (per- 
fect bandpass  is  assumed  for  analytical  ease).  The  degradation  may  be 
derived  approximately  by  multiplying  a sequence  of  ±1' s (with  a trian- 
gular autocorrelation  function)  by  another  sequence  that  has  been  low- 
pass  filtered  at  the  chip  rate  (here,  1.023  MHz),  the  product  has  a 
mean-squared  value  0.4436  dB  down  from  unity,  (The  code  autocorrelation 
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is  approximated  as  triangular  even  though  Gold  codes  are  actually  used 
here.)  This  effect  was  not  accounted  for  in  the  analytical  model  used 
to  obtain  Eq.  (F-26) ; however,  the  effect  of  the  bandlioiting  on  the  in- 
put noise  was  taken  into  account.  Thus  in  practice  Eq.  (F-31)  would 
become 

(461.8626)  RpRN^(A1  - \±)  cos  ♦q  + ng(t)  , (F-32) 

where  n8(t)  has  a variance  of  124326.4  counts^.  Incidentally,  the  ef- 
fects of  bandlimitlng  of  signal  and  noise  on  system  performance  approxi- 
mately cancel  each  other.  Thus  Eq.  (F-3)  of  the  text  has  now  been  de- 
rived. 
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Appendix  F 
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Appendix  G 

SOME  PRACTICAL  SIMPLIFICATIONS 


In  this  appendix,  we  present  some  practical  simplifications  and 
their  effects  based  on  Fig.  A- A. 

First,  note  that  the  operations  of  Lowpass  and  Lowpass  (f) 
in  Fig.  A-3  may  be  performed  after  the  heterodyning  operations  in  Fig. 
A-A.  This  will  permit  the  redundant  heterodyning  operations  on  Output 
(D  and  Output  (b)  to  be  combined  with  the  other  heterodyning  operations 
in  Fig.  A-A. 

Second,  note  that  the  outputs  of  the  suppressed-carrier  tracking 
loop  are  low  frequency  (<  A5  Hz) . The  outputs  of  the  heterodyning  opera- 
tions on  Output  (5)  and  Output  (f)  are  also  low  frequency.  Thus,  we 
could  reduce  the  internal  sampling  rate  in  the  suppressed-carrier  track- 
ing loop,  and  accumulate  samples  of  Output  ($)  and  Output  ® before 
mixing  these  outputs  with  the  suppressed-carrier  loop  outputs.  The  ac- 
cumulation time,  Ts,  would  have  to  satisfy  Ts  <<  0.01  s in  order  not  to 
cause  extensive  lntersymbol  interference.  The  signal  y3i(t),  which  is 
sent  to  the  data  estimator,  would  no  longer  change  value  at  the  ship- 
board sampling  rate,  Tp,  but  every  T8  second  instead. 

We  not  determine  the  intersymbol  effects  of  the  above  procedure. 

We  assume  that  the  number  of  accumulated  samples  per  data  bit  Interval 

A Ara 

is  a fixed  whole  number  (Ng  “ )•  The  proper  choice  of  which  Ng  sam-- 

s 

pies  belong  to  each  data  bit  is  the  function  of  the  data  clock  synchro- 
nizer. The  synchronizer  uses  the  estimated  code  epochs  to  determine 
which  of  the  20  code  epoches  per  data  interval  is  aligned  with  the  data 
bit  epoch.  Accumulated  samples  that  fall  less  than  (Ts/2)  second  after 
the  estimated  data  epoch  are  ascribed  to  the  previous  bit  in  the  bit- 
value  estimation.  The  synchronizer  uses  the  estimated  code  epoches  to 
determine  which  of  the  20  code  epoches  per  data  interval  is  aligned  with 
the  data  epoch.  The  absolute  value  of  the  offset,  ke,  between  the  exact 
data-blt  boundaries  and  the  Interval  over  which  the  Ng  samples  integrate 
lies  approximately  between  0 and  (T8/2)  + o(c^).  The  data  bits  are  es- 
timated by  adding  together  the  Ns  samples  over  a given  data  bit  interval 
to  obtain  a positive  or  negative  answer.  The  effective  signal  power 
degradation  factor,  ds,  to  be  put  into  Eq.  15  is  a function  of  ke.  If 
the  bit  stream  is  a pair-wise  uncorrelated  sequence,  this  factor  if 
given  by: 
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0.  * <T>  • 1 + <£> 


(G-l) 


Thus,  if  T„  ■ 5 x 10“*  and  A„,  - 20  ms,  we  would  insert  da  ss  * multipli- 
cative factor  in  Eq.  15,  where  da  lies  between  1 (for  ka  - 0)  and  0.975 
[for  ka  - 2.5  x 10“4  + ©(e*)  <*  2.5  x 10-*].  We  could  utilise  either  the 
mean  value  of  da,  calculated  by  assuming  a probability  density  function 
for  ke,  or  we  could  just  use  the  lower  limit  (0.975  in  this  case). 


Note  that  according  to  discussion  there  should  be  a factor  of  da 
in  Eq.  15  as  it  stands  in  Section  5 (with  Ta  - Tp,  the  shipboard  sampl- 
ing rate).  However,  this  factor  would  be  very  close  to  unity,  and  so 
has  been  deleted. 
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